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I.  INTRODUCTION  AND  OVERVIEW 

This  report  describes  research  performed  by  Berkeley 
Research  Associates  under  Contract  #N00014-81-C-2371 
with  the  Plasma  Physics  Division,  Naval  Research  Laboratory. 
-The- report  covers  the  period  June  1981  to  June  1983,  during 
which  time  investigations  were  conducted  into  several  major 
topics  concerning  the  behavior  of  an  intense  (multi-kilo- 
ampere)  electron  beam  in  the  modified  betatron  electron 
accelerator. -  The  work  described  here  was  performed  in  sup¬ 
port  of  and  in  close  association  with  the  staff  of  the  NRL 
special  focus  program,  "Advanced  Accelerators." 

- —  The  modified  betatron  has  been  selected  by  NRL  for 

experimental  evaluation  as  a  high-current  electron  accelera¬ 
tor.  Theoretical  support  has  been  directed  at  identifying 
those  phenomena  which  will  most  directly  affect  accelerator 
performance.  To  this  end,  research  has  been  carried  out  in 
the  following  areas:  (1)  Transverse  linear  beam  dynamics  in 
time  varying,  azimuthally  symmetric  fields,  (2)  Effects  of 
grad  Bq  induced  drifts,  (3)  Orbital  resonance  effects  due  to 
small  field  errors,  (4)  Nonlinear  effects,  especially  those 
due  to  non  constant  betatron  field  index,  (5)  Collective 
effects,  especially  the  negative  mass  instability,  and  (6) 
Strongly  focused  systems.  These  studies  were  undertaken 

f 

with  a  view  toward  assisting  in  the  choice  of  parameters  for 
an  experimental  device  to  be  constructed  at  NRL.  The  results 


of  this  work  have  served  as  intended,  to  focus  the  experi¬ 
mental  design  effort  in  that  region  of  parameter  space  most 
likely  to  produce  a  beam  of  the  required  characteristics. 

This  report  is  divided  into  eight  sections.  Follow¬ 
ing  the  introduction,  in  Section  II,  we  describe,  in  outline, 
the  assumptions,  analysis,  and  conclusions  of  our  study  of 
the  dynamics  of  an  intense,  azimuthally  symmetric  beam.  The 
beam  is  studied  in  the  paraxial  approximation  wherein  all 
fields,  both  self  and  externally  applied,  are  taken  to  vary 
linearly  with  displacement  from  the  nominal  design  orbit 
(assumed  to  be  planar) .  This  approximation  allows  us,  by 
performing  an  average  over  an  ensemble  of  initial  conditions, 
to  obtain  equations  for  the  motion  of  the  beam  centroid 
about  the  center  of  the  vacuum  chamber,  as  well  as  for  the 
motion  of  individual  particles  about  the  center  of  the  beam. 

A  WKB  solution  to  these  equations  allows  several  conclusions 
to  be  drawn  about  the  stability  of  the  betatron  oscillations 
and  the  adiabatic  behavior  of  these  oscillations,  both  during 
acceleration  and  the  subsequent  removal  of  the  toroidal 
magnetic  field. 

Section  III  outlines  our  analysis  of  the  effect  of 
the  radial  gradient  of  the  toroidal  field  on  beam  motion — a 
nonlinear  effect.  We  find  grad-B  drift  to  be  canceled  by 
the  weak  focusing  forces,  the  net  effect  being  a  slight  shift 
in  the  equilibrium  beam  position,  unless  it  happens  that  the 


focusing  forces  are  just  canceled  by  wall  image  forces. 

The  avoidance  of  this  transition  point,  where  betatron  focus¬ 
ing  and  wall  image  defocusing  forces  balance,  leads  to  an 
important  limit  on  beam  current  which  is  independent  of  the 
toroidal  field.  This  limit  is  given  in  Section  III. 

Sections  IV  and  V  deal  with  resonance  effects.  In 
Section  IV  we  describe  the  effect  on  the  beam  centroid  of 
small  azimuthal  variations  in  the  applied  fields.  We  obtain 
an  important  bound  on  the  allowable  azimuthal  field  perturba¬ 
tion  in  order  that  certain  so-called  integer  resonances  will 
not  disrupt  the  beam.  The  approach  again  employs  the  paraxial 
equations  of  motion  and  the  averaging  technique  developed  in 
our  study  of  the  dynamics  of  the  symmetric  beam.  The  analysis 
of  Section  V  presents  certain  results  on  nonlinear  resonance 
effects  caused  by  the  transverse  variation  of  the  field  index 
n. 

Next,  in  Section  VI,  some  results  on  collective  insta¬ 
bilities  are  given,  building  on  work  by  Sprangle  and  Vom- 
voridis.  A  novel  result- -a  double  valuedness  in  the  current 
vs  energy  spread  stability  curve — is  predicted  as  a  result 
of  competition  between  growth  and  stabilization  mechanisms. 

A  strongly  stabilizing  effect  of  the  toroidal  field  is 
evident. 

In  the  final  technical  section.  Section  VII,  a  dis¬ 
cussion  is  given  of  the  sensitivity  of  a  weakly-focused  system 
to  average  beam  momentum-vertical  field  mismatch. 
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Consideration  of  this  problem  has  led  to  a  proposal  for  the 
use  of  a  type  of  strong-focusing  coil  arrangement  in  which 
orbits  have  been  studied  in  the  linear  approximation. 

In  all  technical  sections,  calculations  are  only  sum¬ 
marized  or  outlined.  Details  are  relegated  to  Appendices. 
Additionally,  computer  codes  developed  in  the  course  of  this 
work  are  documented  in  Appendices. 

A  final  section,  Section  VIII,  briefly  summarizes  the 
work,  states  our  conclusions,  and  suggests  directions  for 
further  study. 


II.  TRANSVERSE  BEAM  DYNAMICS  IN  THE  MODIFIED  BETATRON 

The  modified  betatron  electron  accelerator  field  con¬ 
figuration  consists  of  a  conventional  weak  focusing  betatron 

field  upon  which  is  superimposed  a  toroidal  magnetic  field. 

1  2 

It  has  been  shown  '  that  this  toroidal  field  greatly  in¬ 
creases  the  amount  of  charge  that  can  be  confined  in  a  device 
of  given  size.  Acceleration  is  accomplished  as  in  a  conven¬ 
tional  betatron,  that  is,  by  changing  the  flux  through  the 
electrons'  orbit. 

The  orbit  of  any  particular  electron  depends,  however, 
not  only  on  the  external  fields  as  in  a  conventional  (low  v/y) 
device,  but  also  on  the  non  negligible  fields  produced  by 
all  other  particles  in  the  system.  These  fields  are  found 
by  solving  Maxwell's  equations  with  the  correct  sources. 

Since  we  do  not  attempt  here  to  calculate  the  particle  dynam¬ 
ics  self-consistently,  we  are  forced  to  make  some  approxima¬ 
tion  for  these  sources  which  we  do  by  taking  the  number  and 
current  densities  as  constants  across  the  beam  cross-section. 
This  approximation  appears  to  agree  fairly  well  with  the 
number  and  current  distributions  found  in  numerical  simula¬ 
tions  and  our  final  results  are  quite  insensitive  to  the 
exact  distributions  which  affect  only  a  certain  coefficient 
in  the  argument  of  a  logarithm.  Maxwell's  equations  are 
solved  through  first  order  in  the  inverse  aspect  ratio  of 
the  torus.  It  is  very  important  to  include  these  "toroidal 


corrections”  to  the  self-fields  for  intense  beams  (v/y^.02). 

They  affect  both  the  value  of  the  required  vertical  field 

needed  to  hold  the  beam  at  its  equilibrium  radius  as  well 

as  the  values  of  the  betatron  frequencies.  These  corrections 
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were  not  included  in  previous  treatments  '  . 

The  resulting  fields,  when  substituted  in  the  equa¬ 
tions  of  motion,  yield  a  paraxial  equation  description  of 
the  beam.  When  solved,  the  equations  illustrate  several 
interesting  phenomena  involving  the  beam.  First,  two  basic 
modes  of  oscillation  exist,  a  "fast"  mode  corresponding 
roughly  to  cyclotron  motion  about  the  toroidal  field,  and  a 
"slow"  mode  corresponding  to  an  IT  x  §  drift  motion  where 
here  the  force  ?  is  due  to  a  combination  of  the  ordinary 
weak  focusing  fields,  image  currents  and  charges  in  the  wall, 
and  hoop  stresses  on  the  ring.  It  may  be  arranged,  by  judi¬ 
cious  choice  of  parameters,  that  the  fast  mode  will  always 
be  stable,  throught  the  injection-acceleration-ejection  cycle 
The  slow  (drift)  mode,  however,  is  more  complicated.  Under 
certain  conditions  the  net  radial  focusing  force  (F)  may 
vanish,  leading  to  a  transition  to  unstable  behavior  fol¬ 
lowed  by  a  subsequent  reversal  in  sign  of  the  drift  motion 
that  is  described  in  detail  in  Appendix  A.  This  transition 
may  be  shown  to  occur  at  a  boundary  in  parameter  space  on 
one  side  of  which  the  toroidal  field  is  essential  for  beam 
stability,  while  on  the  other  side  the  toroidal  field  is 


superfluous  for  stability.  This  "instability  gap"  may  not 

be  so  serious  for  single  particle  motion  since  it  may  be 

shown  that  a  small  expansion  of  the  beam  restabilizes  the 

motion.  For  motion  of  the  beam  centroid,  however,  it  is 

much  more  serious.  In  fact,  we  have  concluded,  on  the  basis 

both  of  this  work  and  other  numerical  studies,  that  the  beam 

must  be  launched  and  accelerated  so  as  to  avoid  passage 

3 

through  the  instability  gap  for  beam  center  motion  .  For¬ 
tunately,  this  does  appear  possible  to  do,  though  it  does 
place  a  limit  on  the  current  that  may  be  accelerated  in  a 
device  of  given  aspect  ratio.  That  limit  is  given  by 

r2 

-f  ns  £  2J7<Va)  <1/2 
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where  rQ  and  a  are  the  major  and  minor  radii  of  the  toroidal 
chamber,  r.  is  the  beam  radius,  and  n  is  the  so-called  self 
field  index.  It  is  interesting  to  note  that  this  limit  is 
independent  of  the  value  of  the  toroidal  field  and  may  in 
fact  represent  a  more  stringent  requirement  than  the  basic 
stability  criterion  for  the  fast  mode,  which  can  always  in 
principle  be  satisfied  by  choosing  a  large  enough  toroidal 
field.  For  example,  according  to  the  above  constraint,  a 
10  kA  beam  must  be  injected  such  that  its  in  situ  energy 
(i.e.,  that  energy  retai  ad  as  ’  letic  energy  by  the  beam 
after  the  beam  has  given  up  some  fraction  of  its  energy  to 


the  fields  within  the  chamber,  after  leaving  the  diode) 
corresponds  to  a  value  of  y  in  excess  of  6.2  for  a  chamber 


with  rQ/a  =  10.  (Toroidal  corrections  modify  this  number 
somewhat;  these  are  described  in  Appendix  A.) 

This  constraint  (II-l)  also  affects  the  workability  of 
certain  injection  schemes  which  are  based  on  drifting  the 
beam  away  from  the  injector  structure  during  one  major 
transit  of  the  machine.  Since  the  drift  frequency  is  pro¬ 
portional  to  the  quantity  (1/2  -  r?n  /a2)  (B  /B  ) ,  one  does 
not  want  to  work  too  close  to  transition  since  a  minimum 
value  of  B  is  required  to  stabilize  the  fast  mode  and  one 
requires  that  the  drift  frequency  be  at  least  as  large  as 


“B>2^ 

A 


where  A  is  the  distance  from  the  injector  port  to  the  beam 
equilibrium  orbit  and  T  is  the  orbital  (major)  period.  (Con¬ 
straint  (II-2)  comes  from  requiring  the  beam  center  to  drift 
two  beam  radii  during  one  transit  around  the  machine — a  mini¬ 
mal  requirement  for  achieving  low  levels  of  loss  due  to  scat¬ 
tering  by  the  injector  structure.)  Both  constraints  (II-l) 
and  (II-2)  favor  a  large  value  for  y  at  injection  time. 

The  instability  gaps  for  particle  and  whole  beam  motion 
have  no  analog  in  a  conventional  betatron.  They  occur,  rough¬ 
ly  speaking,  when  self-fields  become  comparable  to  applied 
fields,  which  is  never  the  case  in  a  conventional,  low  current 


device.  The  modified  betatron,  if  it  operates  as  projected, 
will  be  the  first  cyclic  particle  accelerator  in  which  beam 
self-fields  play  a  significant  role  in  the  particle  dynamics. 

Another  phenomenon  occuring  in  the  modified  betatron 

does  have  an  analog  in  a  conventional  betatron.  This  is  the 

adiabatic  change  in  amplitude  of  the  betatron  oscillations, 

as  external  parameters  (Bz,  BQ,  field  index,  flux,  ...)  are 
4 

slowly  changed  .  The  solution  to  the  linearized  equations  of 
motion  allows  us  to  obtain  explicit  expressions  for  the  ratio 
of  the  beam  radius  at  the  end  of  the  acceleration  to  that  at 
the  beginning.  The  result  is  that  the  beam  undergoes  a 
slight  compression,  as  in  a  conventional  accelerator.  The 
beam  remains  well-behaved  during  all  slow  changes  in  param¬ 
eters  as  long  as  one  avoids  the  instability  gap,  the  bound¬ 
aries  of  which  appear  mathematically  as  two  turning  points  in 
the  WKB  solution. 

Unlike  the  case  in  a  conventional  betatron,  however, 
this  adiabatic  decrease  in  the  betatron  oscillation  ampli¬ 
tude  does  not  help  very  much  in  the  basic  injection  problem; 
that  is,  ensuring  that  the  beam  misses  the  injector  after  one 
turn.  In  a  conventional  betatron  one  has  at  least  a  few 
turns  to  accelerate  the  beam  before  a  particle  returns  to  the 
vicinity  of  the  injector,  since  the  betatron  wavelength  is 
somewhat  greater  than  the  machine  circumference.  With  a 
strong  toroidal  field  in  place,  however,  particles  in  the 


modified  betatron  tend  to  follow  the  toroidal  field  lines 


and  one  must  depend  on  the  slow  drift  to  carry  the  beam  away 
from  the  injection  port  in  only  one  turn.  The  beam,  if 
successful  in  missing  the  injector  on  the  first  turn,  may 
be  subsequently  trapped  by  changing  external  parameters. 

In  the  course  of  these  investigations  on  transverse 
beam  dynamics  in  azimuthally  symmetric  fields,  two  computer 
programs  were  developed  to  assist  in  our  understanding  of 
the  electron  orbits.  First,  a  code  LBE  (for  "Linearized 
Betatron  Equations")  was  written  which  integrates  the  linear¬ 
ized  equations  of  motion  for  either  an  individual  particle 
or  the  beam  centroid  in  arbitrarily  time-varying  external 
fields,  Eq (t) ,  B  (t) ,  B„(t).  The  code  includes  toroidal 
corrections  to  the  self-fields,  assuming  a  given  fixed  beam 
radius.  It  was  used  to  generate  Figs.  4-6  of  Appendix  A. 

The  code  itself  is  documented  by  I/O  description  in  Appendix 
AA,  where  a  listing  is  also  given. 

A  second  code,  a  single-particle  integrator,  was  also 
written  for  the  purpose  of  studying  nonlinear  dynamics  at 
high  energies,  where  self-field  effects  are  less  important. 
This  program,  named  BTRAK,  was  eventually  modified  to  include 
azimuthally  varying  fields  for  our  study  of  resonance  effects. 
(See  Sections  IV,  V,  and  VII.)  Its  documentation,  including 
listing,  appears  in  Appendix  BB. 
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III.  GRAD-Bg  drift  in  the  modified  betatron 

The  effect  of  nonlinearities  in  the  self-  and  applied 
fields  will  modify  the  results  of  Section  II;  in  general, 
however,  these  corrections  are  expected  to  be  small,  and  in 
any  event  are  hard  to  calculate.  If  the  toroidal  field  B^ 
is  strong,  however,  an  important  nonlinearity  to  consider 
is  the  radial  gradient  in  BQ. 

In  general,  an  electron  streaming  along  a  field  line 
executing  small  gyro  orbits  will  experience  a  drift  in  the 
direction  ^B  x  B  which,  in  the  modified  betatron,  is  vertical 
In  the  betatron  however,  the  situation  is  complicated  by  the 
presence  of  the  vertical  field  gradient  which  gives  a  verti¬ 
cal  restoring  force;  motion  in  the  vertical  field  clearly 
cannot  be  treated  in  the  drift  approximation,  since  the 
orbit  size  is  of  the  same  order  as  the  scale  length  of 
To  find  the  true  behavior  in  the  combination  of  vertical  and 
toroidal  fields,  we  must  solve  the  betatron  equations  of 
motion,  including  the  9BQ/9r  term.  Such  a  calculation  has 
been  carried  out;  the  details  appear  in  Appendix  B,  part  V. 

The  conclusion  reached  in  the  Appendix  is  that,  except 
for  the  exceptional  case  in  which  the  net  radial  focusing 
force  (due  to  the  betatron  field  and  image  fields)  vanishes, 
the  9Bg/9r  term  affects  the  motion  only  slightly,  giving  a 
radial  shift  in  the  position  of  the  equilibrium  orbit  and  a 
resulting  slight  change  in  the  betatron  frequencies  (which 


remain  real) .  The  radial  shift  in  equilibrium  may  be  inter 
preted  as  the  result  of  a  balance  between  the  outward  "dia¬ 
magnetic"  force,  which  tends  to  expel  the  beam  from  the 
high  field  region  and  the  inward  radial  restoring  force. 
When  the  radial  restoring  force  vanishes,  the  grad-B  drift 
is  free  to  operate.  Since  the  resulting  drift  is  extremely 
fast,  the  only  reasonable  experimental  alternative  is  to 
avoid  the  vanishing  point  for  the  restoring  force.  The  con 
dition  for  this  has  been  given  in  Section  II,  Eq.  ( II— 1 )  , 
which  provides  a  bound  on  the  accelerator  current  which  is 
independent  of  the  toroidal  field. 


IV.  INTEGER  RESONANCES  IN  THE  MODIFIED  BETATRON 


Throughout  the  analysis  of  the  previous  sections  it 
was  assumed  that  all  applied  fields  were  perfectly  azimuth- 
ally  symmetric.  In  practice,  of  course,  any  actual  accel¬ 
erator  will  have  slight  field  imperfections  in  the  toroidal 
direction,  which  will  be  encountered  periodically  by  each 
electron.  In  the  general  case  these  small  periodically 
applied  perturbations  to  the  electron's  orbit  cause  only  a 
small  response.  However,  if  it  happens  that  the  frequency 
of  the  betatron  oscillations  matches  the  circulation  fre¬ 
quency  (or  an  integer  multiple  thereof)  a  constant  phase 
relationship  is  maintained  over  many  circulation  times  be¬ 
tween  the  particle  motion  and  a  Fourier  component  of  the 
field  imperfection.  The  resulting  "integer  resonance"  can 
cause  an  enormous  buildup  of  betatron  oscillations  and  loss 
of  beam  confinement. 

In  a  conventional  betatron  integer  resonances  do  not 
occur  (neglecting  the  marginally  stable  cases  n  =  0  or  1) 
because  the  betatron  frequencies,  both  radial  and  axial,  are 
necessarily  always  lower  than  the  particle  circulation  fre¬ 
quency.  In  the  modified  betatron,  however,  the  fast  mode 
can  be  resonant.  For  motion  of  the  beam  center  the  8.-th 
resonance  occurs  when  the  ratio  of  the  toroidal  to  vertical 
fields  is 

=  t<*24  +  ~rn.  ) » 


( IV-1) 


where  l  is  the  Fourier  harmonic  number  of  the  field  error. 
Since  both  Bg  and  Bz  will  in  general  be  changing  in  time 
during  an  experiment,  some  of  these  resonances  may  have  to 
be  passed  through.  This  will  necessarily  be  the  case  if 
one  anticipates  removing  the  toroidal  field  prior  to  beam 
ejection.  The  question  then  arises  as  to  how  fast  the  £-th 
resonance  must  be  passed  through  in  order  to  avoid  beam 
disruption. 

To  answer  this  question  the  equations  of  motion  in 
the  presence  of  a  field  error  were  formulated  and  solved, 
assuming  that  the  fields  varied  slowly  over  a  circulation 
time.  The  result,  derived  in  Appendix  C,  gives  a  bound  on 
the  magnitude  of  the  field  error  that  may  be  tolerated. 

The  bound,  expressed  in  terms  of  the  acceleration  rate  (y)  , 
is  rather  restrictive,  in  a  practical  example  that  is  worked 
out  in  the  Appendix,  leading  us  to  speculate  on  ways  that  the 
resonant  effect  might  be  minimized. 

One  possibility  which  immediately  suggests  itself  is 
the  use  of  short  acceleration  times,  thereby  limiting  the 
time  during  which  the  resonance  effect  may  operate.  Very 
short  times  may  be  needed,  however,  since  the  required  accel¬ 
eration  time  for  a  given  final  oscillation  amplitude  scales 

_  2 

as  the  (field  error)  .  (See  Appendix  C,  Eq.  [14].) 

A  second  possibility  for  stabilization  investigated 
in  the  Appendix  is  thermal  spread.  Though  thermal  spread 


in  the  beam  introduces  a  spread  in  resonant  frequencies  and 
subsequent  reduction  of  the  response  of  the  beam  center 
motion  when  passing  through  resonance,  individual  particle 
motion  may  still  be  such  that  an  unacceptably  large  beam 
expansion  occurs. 

Yet  another  possibility,  which  may  be  practical  for 
certain  devices,  involves  maintaining  the  ratio  Bz/B0  con¬ 
stant  throughout  the  acceleration.  This  technique  will  keep 
the  tunes  constant  (save  for  the  tune  shift  due  to  space 
charge,  which  is  small  for  large  B0)  during  the  experiment. 

If  the  application  involves  use  of  the  beam  in  situ,  then 
the  presence  of  a  strong  toroidal  field  within  the  device 
at  the  end  of  the  acceleration  should  not  be  a  problem.  It 
would  probably  complicate  an  ejection  scheme,  however. 

A  final  possibility  that  was  investigated  for  sta¬ 
bilization  is  the  frequency  shifting  effect  of  nonlineari¬ 
ties.  Specifically,  both  the  toroidal  field  and  the  betatron 
field  index  will  generally  vary  with  radial  position.  Since 
the  betatron  frequencies  depend  on  the  values  of  these  quan¬ 
tities,  it  is  possible  that  the  frequencies  will  be  shifted 
sufficiently  by  a  small  (tolerable)  radial  displacement  so 
as  to  detune  the  resonance.  A  rather  strong  radial  gradient 
in  n  may  be  required  to  produce  the  desired  effect,  however. 
Results  of  an  investigation  into  this  question  are  presented 
in  the  next  section. 


In  order  to  be  able  to  predict  when  certain  resonances 
will  be  crossed  for  a  given  time  history  of  vertical  and 
toroidal  fields,  a  computer  program  named  TUNES  was  written. 
This  program  searches  for  a  specified  set  of  resonances  of 
the  form 

nfvf  +  nsvs  =  p  ( 

where  vf  and  vg  are  the  numbers  of  betatron  wavelengths  of 
the  fast  and  slow  oscillation  modes  within  the  machine  cir¬ 
cumference,  and  n.p,  n  ,  and  p  are  integers.  The  search  is 

restricted  to  [n^j  +  |n  |-3  and  p<p  where  p  is  speci- 

i  s  max  max 

fied  by  the  user.  Output  includes  the  resonance  label 
(nf ,  n  ,  p) ,  the  time  of  crossing,  and  the  values  of  various 
parameters  at  crossing.  TUNES  is  documented  in  Appendix  CC. 
Its  use  may  be  helpful  in  identifying  experimentally  observed 
resonance  effects. 


V.  EFFECT  OF  TRANSVERSELY  VARYING  FIELD  INDEX 
ON  SINGLE  PARTICLE  DYNAMICS 

As  discussed  in  the  last  section,  one  possibility  con¬ 
sidered  for  controlling  resonant  response  is  the  intentional 
introduction  of  strong  nonlinearities  (large  values  of  r3n/9r 
where  n  is  the  betatron  field  index)  which  would  result  in 
an  amplitude-dependent  betatron  frequency.  If  this  effect 
were  sufficiently  strong,  the  resonant  response  of  a  particle 
could  be  limited  to  small  values  if  its  finite  amplitude 
oscillations  lead  to  a  detuning  of  the  resonance.  We  have 
examined  the  effect  of  quadratic  nonlinearities,  limiting 
ourselves  to  single-particle  motion  for  simplicity.  Though 
this  does  not  allow  us  to  study  the  general  case  (since 
cubic  terms  also  contribute  to  frequency  shifts) ,  there  is 
a  special  circumstance  in  which  the  contribution  of  the  quad¬ 
ratic  terms  dominate.  This  occurs  when  a  coupling  resonance 
(defined  below)  coincides  with  an  integer  resonance.  In 
this  special  circumstance,  certain  progress  may  be  made 
analytically  in  studying  the  effect  of  resonance  detuning. 
Both  a  coupling  resonance  by  itself  and  this  "coincidence 
resonance"  have  been  examined.  Details  are  given  in  Appendix 
D.  Here  we  describe  the  results  of  this  investigation. 

The  equations  of  motion  of  a  particle  in  azimuthally 
symmetric  fields  are 


x"  +  (l-n)x  =  by'  +  (2n-l-^-)x2  -  (^^-)y2  +  i(x,2-y'2)  (V-la) 


y"  +  ny  =  -bx* 


(2n-n2)xy  +  x'y ' , 


I 
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correct  to  second  order  where  x  =  (r-rQ)/r  ,  rQ  is  the  major 

radius,  b  =  Bn  /B  ,  n  is  the  linear  field  index,  n„  is  the 
0o  zo  2 

second  order  field  index,  and  a  prime  indicates  3/90.  These 
equations  are  solved  perturbatively  to  second  order  in  Appen¬ 
dix  D.  To  linear  order  one  obtains  the  usual  betatron  oscil¬ 
lations  with  frequencies  (Appendix  D,  Eq.  [11].): 


Vf,s 


b2+l± [ (b2+l) 2-4n (1-n) ] 
2 


where  the  subscripts  f  and  s  refer  to  the  fast  (+)  and  slow 
(-)  modes  respectively.  These  single-particle  oscillations 
are  always  stable  if  h(l-n)>0. 

By  inserting  the  linear  solutions  into  the  equations 
of  motion,  we  find  that  the  second  order  correction  remains 
small  unless  it  happens  that  the  following  resonance  condi¬ 
tion  is  satisfied: 


vf  -  2v 


c 


This  condition  turns  out  to  be  a  generalization  of  the  so- 
called  Walkinshaw  resonance  condition  occuring  for  n=0.2  or 
0.8  in  conventional  accelerators  at  which  energy  is  exchanged 


between  radial  and  vertical  oscillation  modes.  This  phenom- 

1 


enon  has  been  observed  in  early  cyclotron  experiments  where, 
due  to  small  vertical  aperture  size,  it  has  led  to  loss  of 
the  beam.  In  the  modified  betatron  the  resonance  is  shown 


in  Appendix  D  to  lead  to  energy  exchange  between  fast  and 
slow  modes,  with  no  noticeable  growth  in  beam  size  under 
typical  experimental  conditions.  Thus  the  resonance  (V-3) 
is  fairly  inconsequential  in  the  modified  betatron. 

This  picture  changes  somewhat  when  departures  from 
azimuthal  symmetry  are  taken  into  account.  It  turns  out 
that  when  n=^  the  generalized  Walkinshaw  resonance  coin¬ 
cides  with  both  ordinary  integer  (t=l)  and  half-integer 
orbital  resonances.  By  including  field  error  terms  we  may 
derive  equations  governing  the  evolution  of  the  mode  ampli¬ 
tudes  for  this  "triple  coincidence"  resonance  and  use  their 
solution  to  study  the  effect  of  frequency  shifts  cn  resonance 
detuning  in  this  special  case.  This  program  is  described 
in  detail  and  carried  through  in  Appendix  D  where  particle 
orbits  under  resonance  conditions  are  illustrated  and  dis¬ 
cussed.  The  basic  result  from  this  analysis  is  that  even 
fairly  strong  gradients  in  n  (i.e.,  large  values  of  n2)  do 
not  adequately  control  the  resonant  response  of  a  single 
particle,  that  is,  the  frequency  shifting  effect  is  too 
small  to  be  helpful  in  the  case  we  have  studied. 

Our  conclusion  from  this  and  the  preceding  section 
seems  clear:  It  appears  to  be  important  to  avoid  machine 
operation  near  low  order  integer  resonances,  the  condition 
for  which  being  (IV-1) . 


VI.  BEAM  INSTABILITIES  IN  THE  MODIFIED  BETATRON 

The  question  of  what  limits  the  current  a  particu¬ 
lar  accelerator  can  carry  in  a  stable  manner  is  a  compli¬ 
cated  one.  In  the  case  of  the  modified  betatron,  the  first 
analyses  (Section  II,  Refs.  1  and  2)  suggested,  based  on 
examination  of  individual  particle  betatron  oscillations  in 
the  self-fields  of  the  beam,  that  particle  motion  could 
be  stable  for  large  currents  if  the  toroidal  field  were 
made  sufficiently  strong.  Later  work  (Section  II,  Ref.  3 
and  Appendix  B) ,  which  considered  motion  of  the  beam  centroid 
led  to  the  discovery  that  the  total  beam  current  must  satisfy 
the  constraint  given  by  (II-l) ,  which  is  independent  of  the 
strength  of  the  toroidal  field.  These  analyses,  however, 
treated  the  beam  as  smooth  and  azimuthally  symmetric.  It  is 
known  that  under  certain  conditions  small  azimuthally  varying 
density  perturbations  can  grow  exponentially  in  time  leading 
to  either  bunched  or  kinked  beams .  Such  longitudinal  and 
transverse  beam  instabilities  in  general  become  more  destruc¬ 
tive  (faster  growing)  the  higher  the  current  and  so  it  be¬ 
comes  important  to  consider  their  current  limiting  effect  in 
the  modified  betatron. 

A  dispersion  relation  for  longitudinal  and  transverse 

modes  for  a  beam  in  the  modified  betatron  has  been  derived 

by  Sprangle  and  Vomvoridis^,  where  a  stability  condition  is 
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also  given.  Building  on  this  work,  Sprangle  and  Chernin 


(Appendix  E)  have  considered  a  slightly  more  general  case, 
taking  into  account  short  wavelength  contributions  to  the 
wave  impedances  and  the  stabilizing  effect  of  finite  ampli¬ 
tude  betatron  oscillations.  In  Appendix  E  it  is  shown  that 
the  presence  of  the  toroidal  field  greatly  enhances  beam 
stability  to  both  longitudinal  and  transverse  modes,  thereby 
greatly  increasing  the  current  limit  over  that  of  a  conven¬ 
tional  betatron;  this  stabilizing  effect  was  noticed  in 
Reference  [1]  where  several  numerical  examples  are  worked 
out.  The  growth  rate  of  the  negative  mass/kink  mode,  for 
example,  scales  as  a  fact  attributed  to  the  inhibiting 

effect  the  toroidal  field  has  on  transverse  motion.  Stabili¬ 
zation  of  both  longitudinal  and  transverse  modes,  for  toroi¬ 
dal  mode  number  1/0,  is  due  to  energy  spread,  or,  more  pre¬ 
cisely,  angular  frequency  spread  in  the  beam;  if  two  parti¬ 
cles,  initially  traveling  together,  separate  by  a  wavelength 
or  more  in  a  (cold  beam)  growth  time,  clearly  the  coherence 
of  the  instability  will  be  lost  and  growth  will  stop.  If 
we  call  the  spread  in  angular  frequency  in  the  steady  state 
beam  Aft  and  the  growth  rate  in  the  absence  of  frequency  spread 
T,  then  we  expect,  on  the  above  grounds,  the  stability  condi¬ 
tion  to  be  given  by 

r  £  i  { Aft  | 

to  within  a  numerical  factor.  In  fact,  it  may  be  rigorously 
shown'*'  that  for  a  beam  with  a  Lorentzian  distribution  of 


canonical  angular  momentum  (VI-1) ,  is  the  exact  stability 
condition. 

The  frequency  spread  Aft  of  the  equilibrium  beam  is 
related  to  the  energy  spread  via  the  single  particle  momen¬ 
tum  compaction  factor;  the  relation  is 

| Aft |  =  |wc|a| (AE/E)  (VI 

1  -1  -2 

where  a  =  ( -  n  )  -  y  ,  w  is  the  cyclotron  frequency, 

4  O  v- 

and  AE/E  is  the  full  width  of  the  energy  distribution.  We 

note  the  importance  here  of  including  self-field  effects, 

represented  by  the  self-field  index  ng,  in  the  definition  of  a 

It  is  the  appearance  of  n  in  a,  which  leads  to  a  novel 

s 

effect,  predicted  in  Appendix  E  on  the  basis  of  (VI-1,2); 

For  low  currents  {ng<<^-) /  a  is  effectively  independent  of 
current  and  so,  since  the  cold  beam  growth  rate  increases 
with  current,  the  beam  energy  spread  required  for  stability 
also  increases.  As  one  continues  to  increase  the  current, 
however,  a  begins  to  increase  significantly,  eventually 
overcoming  the  increasing  growth  rate  beyond  which  point 
increasing  the  current  still  further  results  in  stabilization! 
In  fact,  for  ng=^  ,  virtually  no  energy  spread  is  required. 
There  results  therefore,  from  this  competition  between  growth 
and  stabilization  mechanisms,  a  double  valuedness  in  the 
current  vs  energy  spread  stability  curve,  illustrated  and 
discussed  in  Appendix  E.  The  prediction  of  a  second  stable 
operating  regime  for  accelerators  is  the  main  new  result  of 


this  investigation.  Our  conclusion  is  that  the  toroidal 
field  makes  possible  operation  at  high  currents  by  (1)  re¬ 
ducing  cold  beam  growth  rates  and,  (2)  giving  access  to 
stable,  high-current  sectors  of  the  current  vs  energy  spread 
diagram. 


VII.  STRONG-FOCUSING  SYSTEMS 

Conventional  betatrons  are  weak-f ocusing  accelerators, 
meaning  that  the  wavelengths  of  the  betatron  oscillations  are 
of  the  order  of  the  machine  circumference.  A  second  conse¬ 
quence  of  weak-focusing  is  that  the  momentum  compaction  fac¬ 
tor,  defined  as  the  fractional  radial  shift  in  a  mismatched 
beam  divided  by  its  fractional  momentum  mismatch: 

a  5  (Ar/ro)/(Ap/po)  (VII 

is  typically  of  order  unity.  For  a  conventional  betatron  it 
may  be  shown,  for  instance,  that  a  =  (1-n)  ^ .  As  a  result, 
for  weak-focusing  systems  one  can  typically  tolerate  only 
a  few  percent  momentum  mismatch  before  a  beam  is  lost  to  the 
chamber  wells.  Strong-focusing  systems,  on  the  other  hand, 
have  small  values  for  the  momentum  compaction  factor  and 
betatron  wavelengths  are  much  smaller  than  the  machine  cir¬ 
cumference.  A  strong-focusing  accelerator,  consequently,  can 
tolerate  a  relatively  large  momentum  mismatch.  It  was  the 
discovery  over  30  years  age  of  the  strong-focusing  principle 
which  has  allowed  the  construction  of  the  large  radius 
research  accelerators  in  use  today. 

A  modified  betatron  is  a  strong-focusing  system  with 
respect  to  particle  orbits  about  the  center  of  the  beam, 
but  is  a  weak-focusing  system  with  respect  to  motion  of  the 
beam  centroid  about  the  center  of  the  vacuum  chamber.  As  a 


result,  it  is  necessary  that  the  average  beam  energy  be 
matched  to  the  vertical  field  to  within  a  few  percent.  This 
may  be  done  with  careful  injector  and  magnetic  field  design 
but  design  tolerances  tend  to  be  tight. 

It  was  the  realization  of  this  tolerance  problem 
which  led  to  the  consideration  of  ways  to  implement  strong- 
focusing  in  a  modified  betatron  in  a  way  consistent  with  its 
design.  As  a  result  of  a  study,  it  was  proposed  to  add  so- 
called  £=2  stellarator  windings  to  the  betatron  in  order  to 
obtain  the  beneficial  effects  of  strong-focusing  on  the  beam 
centroid  motion.  This  extra  winding  is  basically  a  continu¬ 
ously  twisted  quadrupole,  the  limiting  case  of  conventional, 
closely-spaced  discrete  quadrupoles  which  we  used  in  alter¬ 
nating  gradient  focusing.  Beam  dynamics  in  the  resulting 
configuration,  consisting  of  a  conventional  weak-focusing 
betatron  field,  a  strong-toroidal  field,  and  an  £=2  stellara¬ 
tor  winding  have  been  analyzed  in  the  linear  approximation 
including  the  effects  of  self-fields.  The  results  are  given 
and  discussed  in  Appendix  F  where  expressions  for  the  betatron 
frequencies  and  momentum  compaction  factor  are  derived.  We 
find  the  results  to  be  encouraging  in  the  sense  that  the  addi¬ 
tion  of  the  stellarator  winding  allows  large  beam  currents 
to  be  confined  and  a  large  beam  momentum  mismatch  (~50%  is 
not  unreasonable)  to  be  tolerated. 


The  addition  of  the  stellarator  winding  is  not  with¬ 
out  some  drawbacks.  As  with  any  strong-focusing  system,  the 
effects  of  orbital  resonances  must  be  carefully  considered. 

As  discussed  above  in  Section  IV,  resonances  are  present  and 
may  be  a  problem  already  in  the  modified  betatron  without  the 
stellarator  field.  Addition  of  the  1=2  winding  introduces 
new  sets  of  resonances  which  must  be  examined.  This  work  is 
presently  in  progress  and  should  lead  to  important  guides  for 
design. 

From  a  practical  point  of  view,  the  stellarator  winding 
introduces  some  other  possible  complications.  Injection  may 
become  difficult  due  to  the  presence  of  the  separatrix, 
though  one  possible  way  to  avoid  this  problem  is  to  introduce 
straight  sections  along  which  to  inject.  Construction,  sup¬ 
port,  and  power  supply  questions  for  the  stellarator  winding 
also  need  to  be  examined.  Some  preliminary  study  of  the 
injection  and  coil  design  questions  in  fact  have  led  to  con¬ 
sideration  of  an  £.=  0  system  which  may  be  preferable  from  the 
point  of  view  of  some  of  these  problems.  The  £.=0  stellatron, 
or  "bumpy  torus  accelerator"  may  have  some  practical  advan¬ 
tages  over  the  1=2  system.  It  is  described  and  analyzed 
in  Appendix  G. 

Despite  possible  drawbacks,  we  conclude  that  strong- 
focusing  systems  show  significant  promise  as  high-current 
accelerators.  Basic  issues  in  orbital  stability  and 
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momentum  bandwidth  have  been  addressed.  Resonance  effects, 
beam  instabilities,  injection,  and  detailed  coil  design 
issues  remain  to  be  studied  further. 
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VIII.  CONCLUSIONS 

The  analysis  of  beam  behavior  in  the  modified  beta¬ 
tron  carried  out  by  Berkeley  Research  Associates  over  the 
past  two  years  has  led  to  a  detailed  understanding  of  how 
such  an  accelerator  should  operate,  its  advantages  compared 
to  conventional  betatrons,  and  its  limitations.  Among  the 
advantages  must  be  listed  its  ability  to  hold  large  currents 
stably  during  acceleration.  We  have  found  that  the  toroidal 
field  greatly  enhances  the  equilibrium  current  that  may  be 
carried  and  also  enhances  the  stability  of  the  beam  to  the 
longitudinal  (negative  mass)  and  transverse  collective  beam 
instabilities  which  may  affect  the  beam.  We  have  found  non¬ 
linear  effects,  including  grad-B  drifts  and  effects  of  trans¬ 
versely  varying  field  index  to  be  negligible  as  long  as  the 
net  radial  beam-focusing  forces  remain  finite  (II-l) .  Orbi¬ 
tal  resonances  may  be  a  problem  in  the  device  unless  thay  are 
passed  through  very  rapidly  or  avoided  altogether;  avoidance 
of  the  low  H-number  resonances,  at  least,  is  probably  essential 
and  possible  to  do  in  some  acceleration  scenarios  which  have 
been  discussed. 

Among  the  limiting  features  of  the  modified  betatron 
must  be  mentioned  the  sensitivity  of  the  position  of  the  beam 
equilibrium  orbit  with  respect  to  its  momentum  mismatch. 

This  feature,  a  consequence  of  the  weak-focusing  betatron 
fields,  may  be  overcome  by  the  addition  of  stellarator  fields. 


the  effect  of  which  is  to  allow  beams  with  a  large  momentum 
mismatch  to  be  confined. 

Studies  of  beam  dynamics  under  this  contract  have 
contributed  in  an  important  way  to  the  selection  of  parameters 
for  the  high-current  injection  experiment  now  being  constructed 
This  experiment  will  challenge  and  refine  our  understanding  of 
the  properties  of  high-current  electron  beams  in  toroidal 
devices.  Once  completed,  the  experiment  should  lead  to  an 
operational  high-current  accelerator . 

Theoretical  issues  which  remain  to  be  addressed  include: 
verification  of  the  double-valuedness  in  the  stability  curve 
discussed  in  Section  VI,  in  a  simple  but  rigorous  (Vlasov- 
Maxwell)  model;  analysis  of  the  expected  radiation  spectrum 
(for  diagnostic  purposes  or  for  radiation  source  development) ; 
and  investigation  of  resonance  effects,  injection  methods, 
and  coil  design  in  the  stellatron.  Work  in  these  areas  is 
presently  being  pursued  in  association  with  NRL  personnel. 
Combined  with  the  efforts  of  the  past  two  years,  this  con¬ 
tinuing  research  will  assist  in  a  significant  way  in  meeting 
the  goals  of  the  NRL  Special  Focus  Program  on  advanced  high- 
current  electron  accelerators. 
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The  linearized  equations  governing  the  motion  of  the  center  of  a  beam  about  its  equilibrium  position  in  a  modified 
betatron,  as  well  as  equations  governing  the  motion  of  an  individual  particle  about  the  beam  center,  are  presented 
and  solved.  Self  field  effects,  including  toroidal  hoop  stresses  and  wall  image  forces,  are  included  in  the  analysis. 

All  fields,  both  self  and  applied,  are  assumed  to  be  azimuthally  symmetric  but  are  allowed  to  have  arbitrary  time 
dependences.  The  solutions  to  the  equations  of  motion  are  analyzed  for  stability  and  conditions  for  stability  are 
obtained.  Further  study  of  the  solutions  illustrates  two  phenomena  of  experimental  interest:  (I)  the  unavoidable 
traversal  of  a  finite  “instability  gap”  in  parameter  space  during  acceleration  and  (2)  the  adiabatic  increase  in  the 
amplitude  of  the  betatron  oscillations  during  removal  of  the  toroidal  magnetic  field,  prior  to  beam  ejection.  By 
careful  design,  the  effects  of  these  phenomena  can  be  reduced  to  insignificant  levels  in  an  actual  accelerator. 


1.  INTRODUCTION 

It  has  been  suggested1-2  that  the  current  carrying 
capacity  of  a  conventional  betatron  accelerator 
might  be  improved  dramatically  by  the  addition 
of  a  strong  toroidal  magnetic  field.  Such  a  field 
acts  to  confine  the  beam  during  injection  and 
early  stages  of  acceleration  when  7,  the  usual 
relativistic  factor,  is  small  and  space  charge  ef¬ 
fects  which  tend  to  expand  the  beam  are  large. 
After  acceleration  is  complete,  7  is  large,  space 
charge  effects  are  small,  and  the  usual  weak  fo¬ 
cussing  betatron  fields  are  sufficient  to  confine 
the  beam;  the  toroidal  field  may  then  be  removed 
to  facilitate  beam  ejection.  In  general  both  ver¬ 
tical  and  toroidal  magnetic  fields  may  be  chang¬ 
ing  simultaneously  during  beam  injection  and 
ejection.  It  is  the  purpose  of  this  paper  to  examine 
the  behavior  of  the  beam  in  such  time-varying 
fields. 

Some  early,  though  unsuccessful  experiments 
using  this  modified  betatron  field  configuration 
were  carried  out  in  England  after  World  War  II;1 
subsequent  analysis4  attributed  the  poor  results 
to  the  injection  method  used  at  the  time  whereby 
significant  numbers  of  electrons  intersected  the 
back  of  the  injector  structure  after  a  few  trips 
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around  the  device.  Recently  other,  more  prom¬ 
ising  injection  schemes  have  been  proposed2  5  to 
take  full  advantage  of  the  focusing  action  of  the 
toroidal  field.  The  resulting  prospect  of  con¬ 
structing  a  very  high  current  (~  1-10  KA)  be¬ 
tatron  has  prompted  the  analysis  presented  here. 

We  shall  derive  and  solve  equations  governing 
the  motion  of  the  center  of  an  electron  beam  con¬ 
fined  in  a  modified  betatron  as  well  as  equations 
governing  the  motion  of  an  individual  particle 
within  the  beam.  Whole  beam  and  single  particle 
stability  criteria  will  be  presented;  the  stabilizing 
effect  of  the  toroidal  field  for  both  beam  and  sin¬ 
gle  particle  motions,  noted  earlier,'  2  will  be  ap¬ 
parent. 

When  the  fields  are  allowed  to  vary  in  time  two 
interesting  phenomena  occur.  The  first  phenom¬ 
enon,  which  occurs  during  acceleration,  has  no 
analogue  in  a  conventional  betatron:  As  the  beam 
accelerates  (7  increases)  the  betatron  makes  a 
transition  from  a  region  in  parameter  space  in 
which  the  toroidal  field  is  essential  to  stability 
(modified  betatron  regime)  to  a  region  in  which 
the  toroidal  field  is  superfluous  to  stability  (con¬ 
ventional  betetron  regime).  It  turns  out  that,  ex¬ 
cept  under  extraordinary  circumstances,  the  sys¬ 
tem  must  pass  through  an  “instability  gap" — a 
region  of  parameter  space,  separating  the  modi¬ 
fied  and  conventional  betatron  regimes,  in  which 
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single  particle  motion  is  unstable,  though  beam 
center  motion  may  not  be,  irrespective  of  the 
magnitude  of  the  toroidal  magnetic  field.  How¬ 
ever,  though  the  size  of  the  instability  gap  is  in¬ 
dependent  of  the  toroidal  field,  the  instability 
growth  rate  within  the  gap  is  inversely  propor¬ 
tional  to  this  field.  We  find  below  that  by  judi¬ 
cious  magnet  design  and  sufficiently  rapid  ac¬ 
celeration,  this  gap  may  be  successfully  traversed 
with  minimal  beam  disturbance. 

The  second  phenomenon  occuring  in  time 
varying  fields  does  have  an  analogue  in  a  con¬ 
ventional  betatron;  this  is  the  adiabatic  change 
in  the  amplitude  of  the  betatron  oscillations.6 
Since  the  frequency  of  these  oscillations  depends 
now  on  both  the  vertical  and  toroidal  fields  a  slow 
change  in  either  is  expected  to  alter  the  amplitude 
of  the  betatron  oscillations.  During  acceleration 
we  find,  as  in  a  conventional  accelerator,6  that 
the  oscillation  amplitude  decreases  as  the  vertical 
field  increases.  If  one  now  considers  removal  of 
the  toroidal  field  prior  to  beam  ejection,  we  find 
that,  as  long  as  the  toroidal  magnetic  field  is  much 
larger  than  the  vertical  field,  the  beam  motion 
will  describe  orbits  of  increasing  amplitude  as  the 
toroidal  field  is  decreased.  Once  the  toroidal  field 
becomes  comparable  to  the  vertical  field,  how¬ 
ever,  the  motion  becomes  more  complicated  and 
the  betatron  oscillations  no  longer  continue  to 
increase  in  amplitude.  We  find  that,  by  careful 
choice  of  field  strengths,  the  ratio  of  the  betatron 
oscillation  amplitude  before  acceleration  to  the 
amplitude  of  oscillation  following  complete  re¬ 
moval  of  the  toroidal  field  can  be  adjusted  to  be 
near  one. 

In  the  following  analysis  we  assume  “perfect,” 
i.e.,  azimuthally  symmetric  fields.  By  neglecting 
the  possibility  of  azimuthal  variation  in  the  self 
fields  (due  to  beam  bunching  or  kinking)  we  omit 
here  consideration  of  a  variety  of  beam  instabil¬ 
ities  that  may  occur;7  by  neglecting  similar  azi¬ 
muthal  variation  in  the  applied  fields  (“field  er¬ 
rors”)  we  neglect  the  effects  of  orbital  resonances. 
These  will  be  addressed  in  a  separate  report.8 

II.  EQUILIBRIUM  RADIAL  FORCE 
BALANCE 

The  geometry  of  the  modified  betatron  is  shown 
in  Fig.  1.  The  field  configuration  is  that  of  an 
ordinary  betatron  with  the  addition  of  a  toroidal 
magnetic  field,  B00,  here  taken  to  be  positive  and 
constant  across  the  minor  cross  section  of  the 


FIGURE  1  Cutaway  view  of  modified  betatron  geometry 


torus.  We  consider  an  electron  beam  of  circular 
cross  section,  as  shown  in  Fig.  2,  with  center 
located  at  ( rc ,  zr)  =  (r0  +  A r,  A z)  where  r0  is  the 
equilibrium  radius  for  the  center  of  the  beam  at 
which  the  electric,  magnetic,  and  centrifugal 
forces  on  a  particle  at  the  center  of  the  beam  are 
in  balance.  We  shall  take  r0  to  be  the  major  radius 
of  the  accelerator  chamber.  In  the  absence  of  self 
field  effects  radial  force  balance  requires  the 
electron  circulation  frequency  at  r  =  r0,  z  =  0 
to  be  given  by 

Bo  Dio 

s  eBi0/myoC  (no  self  field  effects), 


FIGURE  2  Coordinates  of  beam  and  particle  in  modified 
betatron.  Center  of  beam  is  at  (r,  i )  =  (r0  +  Ar,  A:).  Electron 
is  at  (r,  z)  =  (r0  +  Ar  4  8 r,  Az  4  8z). 
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where  Bz0  is  the  value  of  the  applied  vertical  be¬ 
tatron  field  at  the  location  of  the  orbit,  -y0  is  the 
usual  relativistic  factor,  e(>0)  is  the  magnitude 
of  the  electron  charge,  m  is  the  electron  rest 
mass,  and  c  is  the  speed  of  light. 

Self  field  effects  will  modify  Eq.  (1)  however.9 
A  nonneutral  current  ring  produces  both  a  zero 
order  vertical  magnetic  field  and  a  radial  electric 
field.  In  general,  for  a  reference  particle  at  r  = 
r„,  z  =  0,  radial  force  balance  requires 

-7 oroOo2  =  ~  [e/0)  +  ~  r0M^0)J  (2) 

„  where  Erm  and  Bzi0)  arc  the  zero  order  fields  at 
r  —  r0,  z  =  0.  From  Appendix  A,  Eqs.  (A-25c, 
26c,  26d) 

Bzi0)  =  Bz o  -  iT/JoePo  —  Ib  (3) 

TO 

Erm  =  -  irn0e  —  Ie  (4) 

r0 

where  the  notation  is  defined  in  Appendix  A. 

The  terms  proportional  to  lu  in  Eq.  (3)  and  lE 
in  Eq.  (4)  are  toroidal  corrections  to  the  self  fields 
of  a  cylindrical  beam.  They  represent  “hoop 
stresses” — self  forces  on  a  nonneutral  ring  of 
current  which  act  to  expand  the  ring.  Since  we 
do  not  attempt  here  to  construct  a  consistent 
equilibrium  for  the  beam10"12  we  leave  lB  and  lE 
arbitrary  in  the  analysis  below  since  their  precise 
values  depend  upon  the  particular  distributions 
of  charge  and  current  in  the  beam.  Still,  one  ex¬ 
pects  the  leading  order  logarithms  in  the  expres¬ 
sions  for  l„  and  lE,  Eqs.  (A-27,  28)  to  be  correct. 

Using  now  the  zero  order  fields,  Eqs.  (3,  4), 
in  Eq.  (2)  we  may  write  the  condition  for  radial 
force  balance  as 

j  1  +  —  Ib  1  Oo2  —  fi<o0o  +  j  In  =  0  (5) 

L  7o  J  7o  ro2 

where 


approaches  as  v/y0  — >  0  is,  to  first  order  in 
v/70 

0o  ~  Djo  £  1  —  —  It:  +  hi'j  j  ,  (7) 

where  a  =  flzor0/c.  Self  field  effects,  represented 
by  the  v/70  term,  are  seen  to  reduce  the  single 
particle  circulation  frequency  below  that  ex¬ 
pected  for. zero  density;  the  correction  term  can 
be  significant  (20-30%)  in  presently  contem¬ 
plated  devices.  The  general  result,  Eq.  (7),  will 
be  needed  below  in  the  derivation  of  the  first 
order  equations  of  motion. 


III.  FIRST  ORDER  EQUATIONS  OF 
MOTION 


In  this  section  the  equations  governing  the  mo¬ 
tion  of  a  beam  and  motion  of  an  electron  within 
the  beam  are  obtained  and  discussed.  We  shall 
consider  in  detail  only  motion  transverse  to  the 
toroidal  magnetic  field,  assuming  that  all  fields, 
both  self  and  applied,  are  independent  of  0. 

The  equations  of  motion  for  a  particle  in  the 
fields  of  (A-25,  26)  to  first  order  in  the  displace¬ 
ments  from  the  reference  orbit  (r„,  0),  are  derived 
in  Appendix  B.  They  are 
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and  where  to*  is  the  beam  plasma  frequency, 
(4irrt0e2/m7o)l/2.  Here  and  below  ns0  retains  the 
definition  assigned  to  it  in  Eq.  (1). 

Equation  (5)  is  a  quadratic  equation  for  the 
circulation  frequency,  0O.  The  solution  which 
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where 

n  -  r  -  r0  =  Ar  +  br 
Zi  =  z  =  Az  +  8z 

tij  =  Wfc2/(2-y02n.02) 
fteo  =  eBoolmyoC 

and  where  P6,  is  equal  to  the  canonical  angular 
momentum  of  the  particle  at  (r,  z)  minus  the  ca¬ 
nonical  angular  momentum  of  the  reference  par¬ 
ticle  at  (r0,  0),  to  first  order  in  small  quantities. 
It  may  be  shown,  using  the  definition  of  Pe ,  Pe 

=  r^rnyV*  -  ^Ae^,  that 

Pel  =  m-Yoro 

v  1  v  1<9) 

x  I  Veiyo2 - 5  nzoferi - £l;o/flAr  I, 

L  7o «  7o  J 

where  Ve,  =  Ve  -  Ve0. 

As  they  stand  Eqs.  (8a)  and  (8b)  are  not  easily 
solved  since,  before  they  can  be  solved  for  the 
coordinates  of  a  particle  (r  | ,  z,)  the  beam  position 
(Ar,  A?)  must  somehow  be  known  as  a  function 
of  time.  However,  a  set  of  consistent  equations 
for  beam  and  particle  motion  may  be  obtained  by 
performing  an  ensemble  average  of  Eqs.  (8a,  b) 


over  initial  particle  coordinates  and  velocities. 
Denoting  such  an  average  by  brackets  it  may  be 
shown  that,  as  long  as  the  beam  is  assumed  not 
to  kink  (Ar,  Az  independent  of  6),  we  will  have 

(ri)  =  Ar,  (8r)  =  (8r)  =  <8r)  =  0  (10a) 

(zi)  =  Az,  (8z)  =  (hz)  =  (8z>  =  0.  (10b) 

Upon  performing  this  averaging  procedure  on 
Eqs.  (8a,  b)  we  will  obtain  equations  governing 
the  motion  of  the  center  of  the  beam.  These  may 
subsequently  be  subtracted  from  the  original, 
unaveraged  Eqs.  (8a,  b)  to  obtain  equations  gov¬ 
erning  the  motion  of  a  single  particle  within  the 
beam.  Both  resulting  sets  of  equations  may  be 
summarized  by  the  following  single  set: 

x  +  to*2*  =  ne0y  +  ifW  +  P  (Ha) 

y  +  w>2y  =  —  fieo j  -  iHooJC  (lib) 

where  the  various  quantities  appearing  in  Eqs. 
(11a,  b)  are  defined  in  Table  I. 

Equations  (1  la,  b)  are  our  basic  starting  points 
for  the  analysis  to  be  presented  below.  In  the 
following  sections  we  will  derive  and  study  the 
WKB  solutions  to  Eqs.  (11a,  b).  First  we  make 
a  few  remarks  on  the  equations  themselves. 

The  term  proportional  to  x  on  the  left  hand  side 
of  Eq.  (1  la)  and  the  term  proportional  to  y  on  the 
left  hand  side  of  Eq.  (11b)  represent  radial  and 
vertical  focussing  forces  respectively.  In  general 
the  coefficients  of  x  and  y  in  these  terms  are  not 


TABLE  1 

Definition  of  Quantities  Appearing  in  Equations  of  Motion,  Eqs.  (11a, b) 
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equal  which  suggests  that  an  initially  circular 
beam  may  not  remain  circular.  The  value  of  n 
which  makes  these  terms  equal  (the  value  re¬ 
quired  to  maintain  a  circular  beam  cross  section) 
is 

+  02) 

which  depends  on  yu  and  therefore  on  time.  In 
what  follows  we  will  leave  n  arbitrary,  though  we 
shall  assume  implicitly  that  its  value  is  close  to 
ncir.  This  is  necessary  for  self  consistency  since 
we  obtained  the  beam  self  fields  Eqs.  (A-26)  as¬ 
suming  a  circular  beam  cross  section. 

In  the  case  of  constant  fields  Eqs.  (11a,  b)  are 
elementary.  For  this  case  we  have 


where  the  eigcnfrcquencies  (frequencies  of  be¬ 
tatron  oscillations)  are  given  by 

w,2  +  w,2  +  floo 
2 


vertical  magnetic  field.  It  therefore  moves  in  or 
out  slightly  depending  on  the  sign  of  the  energy 
mismatch.  If,  however,  the  radial  focussing 
forces,  represented  by  w2,  happen  to  vanish  the 
behavior  becomes  secular  (no  equilibrium  radius 
exists)  and  the  particle  moves  vertically,  up  or 
down  depending  on  the  sign  of  the  mismatch;  this 
secular  motion  is  just  the  so  called  “curvature” 
or  “centrifugal”  drift. 

The  solution  to  the  homogeneous  part  of  Eqs. 
(11a,  b)  also  becomes  secular  when  to,2  =  0.  In 
fact,  when  ux2  =  0  and  toy2  =£  w2  (n  A  ncjr),  the 
point  m2  =  0  corresponds  to  a  turning  point 
(transition  from  stable  to  unstable  behavior)  in 
the  WKB  solution  presented  in  the  next  section. 
Since  w2  for  particle  motion  will  pass  through 
zero  during  acceleration,  it  becomes  important 
to  examine  the  behavior  of  the  solutions  to  Eqs. 
(11a,  b)  for  time  dependent  fields.  In  general,  for 
slowly  time  varying  fields,  a  numerical  solution 
to  Eqs.  (1  la,  b)  over  the  entire  acceleration  cycle 
is  prohibitive  since  the  numerical  integration  time 
step  must  be  small  compared  to  fl<T0'  which  in 
turn  is  extremely  small  compared  to  typical  ac¬ 
celeration  times.  An  explicit  solution  for  this  case 
is  therefore  essential. 


IV.  MOTION  OF  BEAM  IN  SLOWLY 
VARYING  EXTERNAL  FIELDS 


±  [  ((Oj, 2  +  tOy2  +  flop)2  ~  4h)_t20)y2]l/2  j  1/2 

(14) 

and  where  the  Cj,j  =  1,  2,  3,  4  are  constants. 

Stability  conditions  result  in  the  usual  way  by 
requiring  w2/  >  0.  We  postpone  examination  of 
these  conditions,  however,  until  the  following 
section.  We  note  here  only  that  for  values  of  y 
above  a  value  dependent  on  geometry  ( rb ,  a,  r0> 
n )  but  not  on  beam  density,  the  self  field  contri¬ 
butions  to  ojj,2  and  wy2  fall  off  as  -y0_  1 .  rather  than 
y0-3-  For  whole  beam  motion  the  value  of  y  at 
which  the  v/y0  terms  become  comparable  to  the 
rb2nsla 2  term  can  be  modest  (y  ~  10)  for  typical 
laboratory  parameters  (rh  ~  1  cm,  a  =  10  cm, 
r0  =  100  cm,  n  =  0.5). 

The  particular  solution  in  Eq.  (13)  represents 
physically  for  particle  motion  a  first  order  radial 
shift  of  a  particle  which,  while  located  initially 
at  the  reference  orbit  (r0,  0)  does  not  have  the 
correct  energy  to  be  maintained  there  by  the  local 


A.  Stability  Considerations 

If  the  coefficients  of  the  derivatives  of  x  and  y 
in  Eqs.  (1  la,  b)  are  slowly  varying  during  a  period 
of  a  betatron  oscillation,  the  equations  may  be 
solved  by  the  WKB  method.  (See  Appendix  C.) 
To  leading  order  the  solution  is 

J_  V  A/  ( (wy2  ~  to/)1'2 

u)Aye,  <v/z  \(“/  -  ^2)l/2 

(15) 

X  cxpi/’vr  + /'rf, F(,-,. 

where  the  eigcnfrcquencies  are  those  given  in 
(14)  in  which  now  all  quantities  may  depend  on 
time, 

Wi  “  [k2  +  to  2  +  fliio)2  -  4wA2«i>y2l,M,  (16) 

and  where  the  kernels  Kx{t,  t')  and  Ks.(ty  /')  arc 
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given  in  Appendix  C.  The  Aj,j  =  1,  2,  3,  4  are 
constants  in  this  approximation. 

This  solution,  Eq.  (16),  is  valid  far  from  any 
turning  point,  i.e.  where  any  w,  vanishes.  Turning 
points  will  occur  if  wx2<«)y2  =  0  and  if  =£  wy2 . 
(See  below.)  Initially  we  shall  confine  attention 
to  a  cold  beam  (no  longitudinal  momentum 
spread)  for  which  the  particular  solution  in  (15) 
vanishes  identically.  Later  we  shall  comment  on 
the  effect  of  temperature. 

The  solution  is  unstable  (exponentially  grow¬ 
ing)  in  time  for  such  times  that  Imfioj)  <  0  for  any 
j.  Unstable  behavior  will  occur  therefore  when¬ 
ever  either  of  the  following  conditions  is  violated: 

u>,V2  >0  (17a) 

w*2  +  wy2  +  Aeo  >  2(co^2<oy2)l/2,  (17b) 

Forn  =  nc ir(w*2  =  «y2)  inequality  (17a)  is  trivial 
and  (17b)  gives  the  simplified  stability  condition 

fleo  >  max(0,  -  4to,2).  (18) 

If  n  ^  ncir  then  both  conditions  (17a,  b)  must  be 
simultaneously  satisfied  for  stability.  Condition 
(!7a)  in  particular  cannot  always  be  satisfied.  At 
injection  ns  is  typically  quite  large  and  both  wx2 
and  wy2  for  particle  motion  (and  perhaps  for  beam 
motion)  are  negative.  During  acceleration,  as  70 
increases  ns  decreases  ( ns  —  70  ~3)  and  w*2  and 
wy2  change  sign  (for  different  values  of  y0>  if  n 
nCir);  an  instability  “gap”  therefore  exists 
while  to*2  and  o>y2  have  opposite  signs. 

It  is  important  to  point  out  that  w*2  and  wy2  for 
beam  center  motion  (Re:  Table  I)  may  start  out 
and  remain  positive  throughout  the  injection-ac¬ 
celeration  cycle  while  wx2  and  wy2  for  particle 
motion  change  sign.  We  recall  from  Table  I  that 
the  small  quantity  ( rbla )2  multiplies  ns  in  the 
expressions  for  w*2  and  wy2  for  beam  center  mo¬ 
tion  but  not  for  single  particle  motion.  Therefore 
unless  ns  is  extremely  large  initially,  beam  center 
motion  will  remain  stable. 

The  inequalities  Eq.  (17a-b)  are  illustrated 
graphically  in  Fig.  3.  The  stable  regions  of  the 
(<i>y/fleo)2,  (ux/fleo)2  plane  are  those  shaded  re¬ 
gions  I  and  II  in  the  figure.  After  injection  but 
before  acceleration  both  (<i)y/fl„o)2  and  (w^/fieo)2 
for  particle  motion  are  negative  and  in  region  1. 
In  this  region  the  toroidal  magnetic  field  is  es¬ 
sential  for  stability  (modified  betatron  regime). 
Following  acceleration  both  (u>y/Aeo)2  and  (w,/ 
ft,o)2  are  positive,  i.e.,  in  region  II  in  which  the 


FIGURE  3  The  plane.  Shaded  regions  are 

stable.  Trajectors  a  and  c  pass  through  unstable  regions.  Only 
trajectories,  such  as  b .  avoid  all  unstable  behavior. 

toroidal  field  is  no  longer  required  for  stability 
(conventional  betatron  regime).  Only  by  passing 
precisely  through  the  origin  (e.g.,  trajectory  b  in 
Fig.  3)  can  instability  be  avoided  altogether. 
While  the  size  of  the  instability  gap  does  not  de¬ 
pend  on  the  magnitude  of  Be0  the  value  of  Im(wy) 
in  the  gap  does  and  is  inversely  proportional  to 
B0 o-  Therefore  by  choosing  a  sufficiently  large 
toroidal  field  it  should  be  possible  to  pass  through 
the  instability  gap  safely  (within  a  few  growth 
times,  or  less). 

We  may  be  quantitative  for  a  case  in  which 
toroidal  effects  may  be  neglected:  When  Eq. 
(17a)  is  violated  and  if  ft£0  >  |  w,2  |,  |  wy2  |  then 
for  the  unstable  mode,  from  Eq.  (14), 


which  has  a  peak  value,  assuming  only  70  and 
not  Be 0  is  changing  in  time,  of 


where  tx  and  t2  arc  the  times  at  which  the  insta- 
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bility  gap  is  entered  and  exited,  respectively,  then 
one  expects  that  the  transit  through  the  gap  will 
not  significantly  disrupt  the  beam;  Eq.  (21)  trans¬ 
lates  into  a  constraint  on  -yo 

-  })2.  (22) 

•yo  3  Beo 

If  the  acceleration  is  fast  enough  to  satisfy  Eq. 

(22)  particle  motion  will  be  essentially  unaffected 
by  passage  through  the  gap.  It  should  be  possible 
to  choose  a  machine  design  (i.e.,  a  sufficiently 
large  toroidal  field  and  a  field  index  close  to  i) 
so  that  Eq.  (22)  is  well  satisfied. 

The  instability  which  occurs  while  tox2Wy2  <  0 
has  an  interesting  dynamical  origin.  Let  us  con¬ 
sider  the  equations  of  motion,  Eqs.  (11a,  b),  tak¬ 
ing  F  =  0,  and  taking  the  external  fields  to  be 
constant  in  time: 

x  +  (i>*2jc  =  fteoy  (23a) 

y  +  toy2y  =  -fleoi.  (23b) 

These  equations  are  just  those  governing  the 
motion  of  a  particle  in  an  effective  electric  field 

£*cfr  = u>x2x  (24a) 

e 

Ey‘n  =  -  oiy2y  (24b) 

e 

and  a  magnetic  field  Be oho-  Converting  to  polar 
coordinates  p,  <t>  we  have 

£Pcfr  =  —  pfw*2  cos2  4>  +  <oy2  sin2  <}>]  (25a) 

e 

£<t>efr  -  —  p[o)y2  -  w*2]  sin<j>  cos<{>.  (25b) 

e 

The  particle  behavior  may  be  understood  as  fol¬ 
lows.  Let  us  assume  that  n  >  i,  from  which  it 
follows  that  <j>y2  >  to,2  always,  and  let  us  consider 
first  the  modified  betatron  regime  (to*2  <  0,  Wy2 
<  0).  £perr  in  this  regime  is  everywhere  negative 
thereby  giving  rise  to  a  clockwise  £  x  B  drift, 
assuming  B60  is  positive.  E+n ,  which  is  much 
smaller  in  magnitude  than  £pefr,  gives  a  radial 
drift  of  alternating  sign  as  the  particle  moves  from 
quadrant  to  quadrant,  thereby  producing  an  el¬ 
liptical  orbit.  Stable  motion  is  established  by  bal¬ 
ancing  the  outward  radial  electrostatic  +  out- 
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ward  centrifugal  forces  against  the  V  x  B 
confining  force. 

In  the  conventional  betatron  regime  w*2  >  0, 
toy2  >  0  and  the  sign  of  £pefr  is  reversed.  Azi¬ 
muthal  particle  drift  is  now  counter-clockwise 
and  the  major  axis  of  the  elliptical  orbit  is  rotated 
by  90°.  Stable  motion  is  achieved  by  balancing 
the  inward  radial  electrostatic  force  against  the 
centrifugal  force;  the  toroidal  field  is  no  longer 
needed. 

In  the  instability  gap  £pefr  has  zeroes  at  polar 
angles  given  by 


at  which  points  the  azimuthal  drift  velocity  van¬ 
ishes.  The  radial  drift  velocity,  cE^IBe,  cannot 
also  vanish  at  the  same  point.  Consequently  the 
particle  drifts  radially,  with  increasing  velocity, 
since  £ ~  p,  at  the  angle  4>0,  as  long  as  tox2wy2 
<  0.  Increasing  the  toroidal  B  field,  thereby  re¬ 
ducing  the  radial  drift  velocity,  reduces  the 
growth  rate  of  this  instability,  a  fact  reflected  in 
Eq.  (19). 

Typical  orbits  during  transit  of  the  instability 
gap  are  illustrated  for  a  simple  case  in  Figs.  4  and 
5  in  which  results  of  a  numerical  integration  of 


-1.0  8f  - -  1.0 


FIGURE  4  Particle  trajectory  (&z  vs.  6r)  in  the  modified 
betatron  during  transit  of  the  instability  gap.  y  varies  linearly 
in  time  from  7.0  to  16.1  in  2.4  ps.  /?»  =  600  gauss,  r0  =  100 
cm,  a  =  10  cm,  rk  =  I  cm,  n  -  0.53,  vly  =  8.4  x  10°  at 
t  =  0. 
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Eqs.  (1  la,  b)  are  plotted.  In  Fig.  4  condition  (21) 
is  not  well  satisfied.  The  dramatic  drift  direction 
reversal  and  instability  are  evident.  In  Fig.  5  con¬ 
dition  (21)  is  well  satisfied  ( n  is  near  i);  particle 
motion  is  virtually  unaffected,  except  for  the  rev¬ 
ersal  of  drift  direction,  by  passage  through  the 
gap.  The  two  graphs,  in  Figs.  4  and  5  differ  only 
by  the  value  of  n  used;  all  other  external  param¬ 
eters  and  total  integration  time  are  identical. 

So  far  no  mention  has  been  made  of  the  effect 
of  temperature,  the  inhomogeneous  term  in  Eqs. 
(11a,  b),  on  particle  orbit  behavior  in  or  near  the 
instability  gap.  Particles  having  an  energy  mis¬ 
match — either  too  little  or  too  much  energy  to 
be  maintained  at  the  reference  orbit  by  the  local 
vertical  field — will  seek  out  their  new  equilibrium 
orbits  about  which  they  will  execute  betatron 
oscillations.  Secular  behavior  is  expected,  as  dis¬ 
cussed  earlier,  when  vanishes. 

The  effect  of  energy  mismatch  on  a  particle 
orbit  is  illustrated  in  Fig.  6  where  the  particle  of 
Fig.  5  has  been  given  an  energy  mismatch  of 

Pei  -  (Pei)  ,  v 

- «7,  -  (-y,)  =  0.10. 

mr0c 

The  effect  is  twofold.  The  orbit  center  shifts 
slightly  outward  and  the  amplitude  of  betatron 
oscillations  following  passage  through  the  insta¬ 
bility  gap  has  increased  by  a  factor  of  —35  over 


-0.3  l- 
-0.3 


FIGURE  5  Particle  irajeclory  (Sz  vs.  hr)  in  the  modified 
betatron  during  transit  of  the  instability  gap.  All  parameters 
are  as  in  Fig.  4  except  n  =  0.51. 
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FIGURE  6  Particle  trajectory  (Sz.vs.  hr)  in  the  modified 
betatron  during  transit  of  the  instability  gap,  including  energy 
mismatch.  All  parameters  are  as  in  Fig.  5  except  an  energy 
mismatch  of  (Ft!  -  (Pei))tmr0c  =  0.10  has  been  introduced. 


the  zero  mismatch  case.  Such  a  large  expansion 
of  the  particle  orbits  cannot,  in  fact,  be  reliably 
computed  using  the  linearized  Eqs.  (1  la,  b)  used 
here.  One  non-linear  effect  in  particular,  namely 
the  reduction  of  beam  density  during  the  orbit 
expansion,  will  clearly  speed  the  passage  of  a 
particle  through  the  instability  gap.  (Recall  that 
n,  is  proportional  to  density.)  Due  to  this  density 
reduction  the  actual  degree  of  orbit  expansion  to 
be  anticipated  in  a  real  device  is  likely  to  be  sig¬ 
nificantly  less  than  that  seen  in  Fig.  6.  Still,  these 
calculations  suggest  that  a  fairly  cold  beam  will 
be  required  for  successful  acceleration  through 
the  instability  gap.  Poorly  “matched”  particles 
are  likely  to  be  lost  as  w*2  goes  through  zero.  It 
should  be  pointed  out  as  well  that  a  strong  to¬ 
roidal  field  greatly  reduces  the  effects  of  energy 
mismatch.  The  computer  runs  necessarily  em¬ 
ploy  a  very  modest  toroidal  field  (600  gauss  in 
the  case  of  Figs.  4-6)  due  to  time  step  consid¬ 
erations.  A  stronger  field,  by  further  restricting 
radial  motion,  is  expected  to  improve  the  con¬ 
finement  properties  of  a  warm  beam. 

B.  Adiabatic  Behavior 

Let  us  next  briefly  consider,  using  the  solutions 
to  the  equations  of  motion,  Eq.  (15),  the  effects 
on  the  particle  orbits  of  the  removal  of  the  to- 
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roidal  magnetic  field.  The  toroidal  field  may  need 
to  be  removed  in  order  to  facilitate  beam  extrac¬ 
tion  though  this  may  not  be  essential.  Let  us  as¬ 
sume  that  Eq.  (15)  is  valid  throughout  the  accel¬ 
eration  cycle,  i.c.,  that  or,2  and  to,2  pass  through 
zero  simultaneously  and  that  the  solution  to  the 
homogeneous  equation  (the  sum  in  Eq.  (15))  dom¬ 
inates  the  solution.  This  is  certainly  true  for 
matched  particles  (P6  =  ( Pe )  =  0)  when  n  =  $ 
and  when  toroidal  effects  may  be  neglected  {vly 
<  1).  One  may  show,  using  Eq.  (15)  for  such  a 
case,  that  for  beam  center  motion  in  either  the 
fast  or  slow  oscillation  mode 

l(Ar)2  +  (Az)2]/ 
l(Ar)2  +  (Az)2], 

(27) 

I 


while  for  particle  motion  about  the  beam  center13 

[(Sr)2  +  (8z)2]/ 
i(Sr)2  +  (5z)2], 

(28) 

[(*  ~  ns)B\ o  -f  j  flip],- 1 1,2 
l(i  -  ns)  Bio  +  1  Blo\f\  ' 

where  the  subscripts  /  and  /  correspond  to  any 
initial  and  final  states.  The  latter  expression,  Eq. 
(28),  may  be  interpreted  as  the  fractional  change 
in  beam  cross  sectional  area.  Note  that  for  large 
B9  the  area  of  the  orbits  —Bq~1,  as  expected. 

Expressions  for  these  ratios  in  the  case  that 
toroidal  effects  arc  not  negligible  and  n  J  may 
be  obtained  from  Eq.  (15).  The  expressions  arc 
complicated,  however,  and  will  not  be  cited  here. 

As  a  numerical  example  we  consider  a  1  kA 
beam  of  1  cm  initial  radius  in  an  initial  state  cor¬ 
responding  to  7,  =  7,  Bt 0>,  =  120  g,  B9o.i  -  L5 
kg  and  a  final  state  with  =  100,  Bz0,/  =  1.7 
kg,  and  B00j  =  0.  In  such  a  case  Eq.  (27)  gives 
for  the  orbital  area  ratio  a  value  of  0.63  while  Eq. 
(28)  gives  for  the  ratio  of  beam  cross  sectional 
areas  a  value  of  0.60. 

We  conclude  that  it  should  be  possible  both  to 
accelerate  the  beam  and  to  remove  the  toroidal 


field  to  facilitate  beam  ejection  without  causing 
either  the  beam  orbit  or  individual  particle  orbits 
to  expand  without  limit. 

V.  CONCLUSIONS 

The  beam  in  a  modified  betatron  can  be  stably 
confined  both  during  the  acceleration  phase  and 
during  the  subsequent  gradual  removal  of  the  to¬ 
roidal  magnetic  field  prior  to  beam  ejection.  As 
the  beam  is  accelerated,  however,  unless  very 
special  conditions  are  satisfied,  a  region  of  insta¬ 
bility  will  be  passed  through;  however  if  the  time 
of  transit  through  this  instability  gap  is  small  com¬ 
pared  to  the  time  specified  in  Eq.  (20)  the  net 
effect  should  be  small. 

As  the  toroidal  field  is  removed  to  facilitate 
beam  extraction  following  acceleration  no  further 
instability  gaps  occur  but  the  magnitude  of  the 
beam  betatron  oscillations  will  change  adiabati- 
cally.  By  arranging  that  the  ratios,  Eqs.  (27,  28), 
be  near  one,  one  expects  the  beam  to  be  well 
behaved  during  the  removal  of  the  toroidal  field. 

It  should  be  remarked  however  that  changing 
the  toroidal  field  changes  the  “tune"  of  the  be¬ 
tatron  which,  in  general,  will  necessitate  the  pas¬ 
sage  through  orbital  resonances  as  the  toroidal 
field  is  removed.  These  resonances,  due  to  the 
periodic  encounter  by  a  particle  of  a  field  error 
or  “bump”  are  currently  under  investigation.  It 
is  anticipated  that  a  condition  governing  the  min¬ 
imum  speed  with  which  B0  must  be  removed, 
expressed  as  a  function  of  the  magnitude  of  the 
field  error,  will  be  obtained. R 
Note  added  in  proof:  Due  to  a  quirk  in  the  pub¬ 
lication  process  the  work  of  reference  8,  while 
completed  and  submitted  for  publication  after  the 
present  work,  actually  appears  in  print  earlier  in 
this  volume  (Part.  Acc.  12,  329  (1982)) . 
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APPENDIX  A 

Fields  in  the  Modified  Betatron 

In  this  appendix  we  calculate  the  fields  seen 
by  a  particle  in  a  modified  betatron.  The  particle 
is  assumed  to  be  close  to  the  axis  of  the  torus, 
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that  is,  the  coordinates  of  the  particle  are  taken 
to  be  (refer  to  Figs.  1  and  2  in  the  text) 

(r,  z)  =  (r0  +  Ar  4-  hr,  A z  4-  8z) 

and  all  fields  will  be  calculated  to  first  order  in 
Ar,  8r,  A z,  and  8z.  Fields  will  be  given  in  the  (r, 
6,  z)  coordinate  system  of  Fig.  1  and  all  will  be 
assumed  to  be  independent  of  6.  Superscripts  a 
and  s  will  be  used  below  to  denote  applied  and 
self  fields,  respectively. 


Part  I  (Applied  Fields) 
Magnetic  Field 


Electric  Field 

All  applied  electric  fields  are  inductive.  The  to¬ 
roidal  electric  field  is  governed  by  the  changing 
central  flux  and  is  taken  to  be  a  specified  function 
of  time 

£e"  =  Eeo(t).  (A-4) 

Ee  is  negative  for  electron  acceleration  with  Bi0 
positive. 

Changing  the  toroidal  magnetic  field,  B00,  will 
induce  a  poloidal  electric  field,  the  r  and  z  com¬ 
ponents  of  which  are  easily  found 

Er°  =  -  ^-£»o(Az  +  8z)  (A-5) 


The  usual  weak  focussing  betatron  field  has  r  and 
z  components.  The  z  component  is  taken  to  be¬ 
have  near  r0  as 


B°  ~  Blo(,rofr)n  (A-l) 

„  Ar  4-  8r\ 

-"—7r~h 

where  Bl0  depends  only  on  time  and  n,  taken  as 
a  constant  to  this  order,  is  the  so-called  vacuum 
field  index.  The  radial  field  is  obtained  by  re¬ 
quiring  ($  x  i)e  =  0  and  Br(z  =  0)  =  0  (making 
the  z  =  0  plane  a  plane  of  symmetry).  The  result 
is 


•(A'2) 

The  applied  toroidal  field  generally  falls  off  as 
r-1  across  the  minor  cross  section  of  the  torus 

where  Bb0  depends  only  on  time.  However,  in 
the  equations  of  motion  Be  multiplies  only 
r  and  i  terms  which  are  already  first  order. 
Therefore  the  gradient  of  Bb  does  not  enter  the 
linearized  equations  of  motion  and  we  take  only 
the  zero  order  value, 


Eza  =  ^-£eo(Ar  +  hr),  (A-6) 

2c 

where  a  dot  indicates  a  time  derivative. 

Part  II  (Self  Fields) 

Since  we  neglect  beam  diamagnetism  and  the 
possibility  of  a  change  in  self  flux  due  to  time 
varying  beam  current  we  take  Bb  =  Ees  =  0.  It 
remains  to  calculate  the  r  and  z  components  of 
the  beam  self  electric  and  magnetic  fields. 

Consider  a  beam  circulating  inside  a  perfectly 
conducting  toroidal  chamber  of  circular  cross 
section  as  shown  in  Fig.  A-l.  (The  beam  dis- 


Bb“  —  B00 . 


(A-3) 


FIGURE  A-l  Geometry  for  self  field  calculation 
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placement  is  exaggerated  for  clarity;  we  will  as¬ 
sume  A  <?  a).  The  chamber  major  and  minor  radii 
are  r0  and  a  respectively.  The  beam  major  and 
minor  radii  are  Rb  and  rb  respectively.  We  must 
calculate  the  fields  inside  the  beam  (p  <  rb ),  as¬ 
suming  the  chamber  is  a  perfect  conductor.  To 
proceed  we  define  a  scalar  potential  4>(p,  «J>)  and 
a  magnetic  flux  or  stream  function  'kfp,  <J>)  =  rAa 
where  Ae  is  the  usual  vector  potential.  The  equa¬ 
tions  for  O  and  ¥  are 


1  £/  dd>\  1  a2<i> 

pdp  \p  dp)  +  p2a<j>2 


f  .  ad>  1  .  ,  ad)] 

cos  <t>  — - sin  4>  — 

L  dp  p  d<l>J 

R  b  +  p  cos  4> 


=  4ire/io(p,  4>)  - 


lj9_/  J_« 

pap  \p  ap/  +  p2, 


- - (Rb  +  pcos<j))7e(p,  4>) 

c 


1  • 

-  -sin4>  — 
p 


dp  p 
Rb  +  pcosc{> 


where  n0,  the  beam  number  density  and  Je,  the 
beam  current  density,  are  assumed  to  have  been 
specified.  Here  we  shall  take  both  n0  and  Je  con¬ 
stant,  independent  of  p  and  <{>. 

The  boundary  conditions  on  d>  and  ^  are  the 
same;  they  both  must  vanish  at  the  surface  of  the 
chamber,  specified  by 

p  =  a  -  A  cos  (i|»  -  4>).  (A-9) 

correct  to  first  order  in  A  la. 


Scalar  Potential  and  Electric  Field 

The  general  solution  for  d>,  including  the  first 
toroidal  correction,  is 


d>0  +  <7(1  -  p 2lrb7)  +  — — j  cos  <j> 

4  Rbrb 

P  P 

+  A  —  sin  4>  +  B  —  cos  <{>  p  <  rb 
rb  rb 


(A- 10) 


d>0  -  2 q  In  p lrb  +  ~  In  p !rb  cos  <j> 
Rb 


(a‘  —  +  C'  — ^ 
V  rb  p) 

(b'  -2-  +  D'  -) 
\  rb  p) 


cos  <}>  p  >  rb 


where  q  =  - en0-nrb 2  and  <!>„,  A,  B,  A',  B',  C', 
and  D'  are  constants. 

Applying  now  the  correct  boundary  conditions 
'  ’  both  at  the  beam  surface  and  the  wall  determines 


all  of  the  constants: 
d>0  =  2q\n  alrb 

A  -  A'  -  -  2q— —  sin 4» 

a  a 

*-'•- -*£'*** 

,  Ar* 

-2  q - cos  41 

a  a 

C’  =  0 

D'  =  . 

4  Rb 


(A- 11a) 
(A- lib) 


(A-l  lc) 


(A- lid) 
(A-l  le) 


Using  this  result  in  Eq.  (A-10)  we  may  calculate 
the  r  and  z  components  of  E 5  inside  the  beam, 
to  first  order: 

\ 

rr  j  d<J>  .  1  3«I>  .  , 

Er  =  -  —  cos  4>  +  -  —  sin  <J) 

ap  p  d<p 

-  %  [hr  +  ri ir] +  i ,n  z  (A-'2) 

rs  dd>  .  l  ad> 

Ei  -  ~  —  sin  4>  -  -  —  cos  <t> 

dp  p 


D.  CHERN1N  AND  P.  SPRANGLE 


where  (8r,  8z)  =  p  (cos  <j>,  sin  4>)  and  (Ar,  A z) 
=  A  (cos  sin  t}<). 


Magnetic  Flux  (or  Stream)  Function  and 
Magnetic  Field 

The  general  solution  for  ty,  including  the  first 
toroidal  correction,  is 

[  *0  +  0(1  -  P  i/r*2)  cos  $ 

I  4  Kb  T b 


-  P  .  .  -  p 

+  A  —  sin  '<})  +  B  —  cos  <j)  p  <  r* 
Tb  rb 


The  resulting  magnetic  field,  to  first  order,  is 


(A- 14) 


¥  =  np 

-  20  In  p/rb  -  — -  In  p/rb  cos  «{> 
Rb 


+  (k‘  —  +  C‘  — ^  sin  <j> 

V  rb  p/ 

+  ^2?'  ~  +  D'  cos  4>  p  >  rb 


where  0  =  'nrb2JBRb/c  =  -  Ttrb2en0$oR_b>  Po 
=  V'eo/c,  and  ^o.  A,  B,  A',  B’,  t" ,  and  £>'  are 
constants. 

Applying  the  boundary  conditions  gives 


¥o  =  20  In  - 
rb 

A  =  A'  =  -  20-  — sin 

a  a 


B  =  B'  +  0 


=  -  20  —  —  cos  ij>  +  0  —  In  — 
a  a  T  *  Rb  rb 


(A- 15a) 


(A- 15b) 


„ ,  1  fW  .  ,  ay  cos  4>\ 

BS  =  -  -  (  —  sin  <{>  +  — - ) 

r  \  Dp  3(J>  p  / 

=  -  27re/ioPo  ^Sz  +  Az^ 

„ ,  I  (dV  ,  3V  sin  cj>\ 

r  V  ap  a<j)  p  / 


(A- 16) 


=  2iten0po 


(A- 17) 


In  the  circulation  frequency,  0O,  rather  than 
the  current  itself,  is  taken  to  be  constant  across 
the  beam  (current  ~  r)  then  it  is  straightforward 
a 

to  show  that  In - Pi  in  Eq.  (A-17)  is  replaced 

by  In  —  +  2. 
rb 

If  the  magnetic  field  of  the  beam  has  diffused 
completely  through  the  wall  then  the  field  sur¬ 
rounding  the  beam  is  most  directly  calculated 
using  the  free  space  Green’s  function 

A(f)  =  -  \  dV  (A- 18) 

c  J  I  r  -  r  I 

If  J  -  Je  6  is  constant  across  the  beam  and  JB 
is  independent  of  0  then  A  =  A06  where 

Ae(r,  z)  =  —  I  r'dr’  \  dO'  dz' 

C  JRb-rk  J0  J-Zb(r) 


_ cos  9' _ 

l r 2  +  r'2  -  2rr'  cos  0'  +  (z  -  z')2]’ 


zb(r')  =  [rb2  -  (r’  -  K*)2)’ 


(A-19) 


The  integral  over  0'  may  be  expressed  in  terms 
of  the  complete  elliptic  integrals' 

4/e  +  ...  .  , 


4/b  r««>  +  r» 

Ae  = -  I  r'dr'  J  dz' 

C  J  Rb  -  rb  J-ZbW) 


_  1 J2  ^ 
Rb  4  Rb  a2 


C'  =  0 


(A- 15c)  *  ((r  _  r'Y  +  (Z  -  *»)*]•« 

(A-15d)  X  [(*  +  m)  ~  (A'20) 

(A-15e)  '  Handbook  of  Mathematical  Functions,  M.  Abramowitz 

and  1.  Stegun,  eds.  Dover  Publications,  ch  17. 
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where 

_ 4rr^ _ 

m  ~  (r  ~  r')2  +  (z-  z'f 

In  the  beam  interior  m  is  large.  Using  the 
asymptotic  expansions  for  K  and  E  one  may  show 
that 

(  1  +  —  )  K(  -  m)  -  —E(  -  m) 

\  mj  m 


(A-21) 

~  |m~,/2(ln  m  +  41n2  -  4). 


Using  Eq.  (A-21)  in  Eq.  (A-20)  the  resulting 
integrals  are  elementary.  The  result,  for  the  vec¬ 
tor  potential  inside  the  beam  is 

'*•' "  i  [2  ln  (^) + 1  “ pl/r,! 

(A-22) 

+  ^cos"’(H"?^43) 

]• 

where  I  ~  it rh2J0,  from  which  it  follows  that  the 
fields  inside  the  beam,  to  first  order  in  p,  are 

2/  p  .  t 

Br  = - 5  sin  4> 

c  rb 

(A-23) 

b,--  ' 

crb  IRb  rb  rb 

.  (A-24) 

We  may  summarize  all  of  the  foregoing  results 
as  follows: 

The  applied  fields  are. 

B,‘  -  -  + 

ro 

(A-25a) 

Be"  —  Be  o 

(A-25b) 

(A-25c) 

Er°  =  -  ^-5oo(Az  +  5 z) 

2c 

(A-25d) 

Eea  =  Ee  o 

(A-25e) 

Et°  =  ^-Beo(Ar  +  hr) 

2c 

(A-250 

where  Bz0,  BM,  and  £„0are  taken  to  be  prescribed 
functions  of  time. 

The  self  fields  are 

B*  =  -  2-iinoc3o^8z  ^-Az^ 

(A-26a) 

Be'  =  0 

(A-26b) 

B*  =  2TWioe(io^&r  +  ^  Ar  - 

(A-26c) 

E,'  =  -  2Ttn0e(br  +  r  + 

(A-26d) 

Eb'  =  0 

(A-26e) 

Ez'  =  -  2-nn0e  ^bz  +  ^  A 

(A-260 

where 

In  —  +  2  if  circulation  frequency,  0, 

Tb 

is  constant  across  the  beam 

Ib  = 

In  —  +  1  if  current  density  is  constant 
rb 

across  the  beam  (A-27) 

lE  =  in  —  if  density  is  constant  across  the  beam . 
rb 

(A-28) 

For  times  long  compared  to  the  time  it  takes 
the  magnetic  field  to  diffuse  through  the  chamber 
wall  the  result  (A-24)  shows  that  one  must  replace 
a  in  the  logarithm  in  the  definition  of  /„  by  (8  r0l 
e)  2.9r0.  This  suggests  that  fields  in  an  actual 
device  may  have  to  be  programmed  in  time  to 
compensate  for  this  extra  change  (reduction)  in 
Bt,  in  order  to  hold  the  beam  in  place. 
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APPENDIX  B 

Linearized  Equations  of  Motion  for  a 
Particle  in  the  Modified  Betatron 

In  this  appendix  the  equations  of  motion  for  a 
particle  in  the  fields  (A-25)  and  (A-26)  are  ob¬ 
tained,  correct  to  first  order  in  small  quantities. 
The  complete  equations  of  motion  are 

~-(yr)  -  y rO2  =  -  —  Fe,  +  -  (rQBz  -  zBa) 1 

at  m  L  c  J 

(B-l) 

-4(yr2e)  =  -  -  \e»  +  -  (zBr  -  rBz) 1 

rdt  m  L  c  J 


4(7 il- 

at  me 


We  consider  first  the  linearization  of  Eq.  (B- 
2).  This  equation  has  an  exact  first  integral,  as¬ 
suming  the  fields  do  not  depend  on  6;  it  is  the 
canonical  angular  momentum 


.  e 

■  rnyrH - Ae  • 

c 


We  now  write  all  quantities  Q  as  Q  =  Qo  + 
Q,  where  Q,  <S  £?0  and  Q0  refers  to  quantities 
evaluated  at  the  reference  orbit  (r,  z)  =  (r0,  0). 
Defining  Ve  =  r0  it  is  straightforward  to  show 
from  Eq.  (B-4)  using  (A- 14)  for  (rAeJ)(  that 


Vei  = 

my o 


r,  f. 

- _  o0 - - 

ro  yo  L  myoCj 


Ar  v  . 

-i - ;  llzo  —  IB, 

yo  yo 


where 


BfJ  =  Bto  -  -ir/ioePo  —  Ib 
fo 

v  1  (  ,  e2  \  1  0i2rb 2 

—  =  —  I  ttrb  n0  — 2 )  ~  1  _2 

•yo  70  \  me2/  4  c‘ 

on,2  =  4'trnoe2/myo 


and  where  we  have  used  7,  =  Ve07oiV'ei/ci. 
Now,  using  the  expression  for  0O  in  Eq.  (7),  one 
obtains 

Pe,  r,  v  I 

vei  -  — —5 —  +  5  ^2UzoIe 

my0  r0  70  7o  a 


Ar  v  , 

H - j  fljo  '  In 

7o  7o 

where  a  =  Clz0r0lc. 

The  expression  Eq.  (B-6)  will  be  needed  next 
in  the  linearization  of  the  radial  equation,  (B-l). 
Carrying  out  a  straightforward  linearization  of 
Eq.  (B-l),  using  the  zero  order  fields  from  Eqs. 
(A-25b,  c)  and  Eqs.  (A-26c,  d)  gives 

r.  =  _  |£r)  +  ^2^,1  +  fteoz, 

my0  L  c  J 

-  n  -  0o2ri 

70  (B-7) 

+  Veifljo-yo2 


where 


Heo  =  eBe0lmy0c. 


Using  now  Eq.  (B-6)  and  keeping  terms  only 
to  first  order  in  vly0  and  using  Eqs.  (A-25c,  d) 
and  Eqs.  (A-26c,  d)  to  write 

-  —  [Er,  +  —  Bzl  1 
m  70  L  c  J 

eBoo  ,  wh2  /  rh2  \ 

■ssii  +  vr?lr) 


+  nOofijof  1  (B-8) 


we  obtain  our  final  result  for  the  radial  equation: 

f,  +  —  n  +  fl2o  [l  -n  — — 

70  L  yo 

x  ^(3  -  n)~ih:  +  (2  -  n)/*)]r, 
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eBp o 
2nvyoc 


Zi  +  fleoZi  +  nsCllo 


X  f 8r  +  —  Ar^  +  —  fl20/«Ar 
V  «  /  Yo 

JL 

Yo  \  a2 


+  ftzO 


61 


YowroL 


(B-9) 


where  n,  =  <ot2/2"yo2fl?o  is  the  “self  field  index." 

The  analogous  linearization  of  the  z  equation 
[Eq.  (B-3)]  is  completely  straightforward,  using 
the  fields  (A-25a,  b,  f)  and  (A-26a,  b,  f).  The  result 
is 


Zi  +  ~  Zi  +  nfl|( 


Yo 


zO 


x[1-i(i'E  +  '")]1'  (B10) 


eB 
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2my0c 


fi  “  fW, 


n,fifo  Ar^  . 


APPENDIX  C 

WKB  Solution  of  Equations  of  Motion 

The  linearized  equations  of  motion  are  given  in 
the  text,  Eqs.  (11a,  b).  Below  we  shall  obtain 
first,  an  approximate  solution  to  the  homogene¬ 
ous  version  of  Eqs.  (Ua,  b),  assuming  that  all 
coefficients  are  slowly  varying.  We  will  then  give 
the  solution  to  the  full,  inhomogeneous  equa¬ 
tions. 

The  homogeneous  equations  are 

Jr  +  u>x2x  =  fieoy  +  ifieo  y  (C- 1  a) 

y  +  <o,2y  =  -ft8ox  -  Ifieo*  (C-lb) 

All  coefficients,  <ox2,  ioy2,  and  fi*,  will  Be  as¬ 
sumed  to  vary  significantly  only  over  a  slow  time 
scale.  To  carry  out  a  formal  asymptotic  expan¬ 
sion  then  we  define 


t  =  //A, 


where  A  is  a  large  dimensionless  parameter.  Dc- 
d 

noting  —  by  a  prime  ('),  Eqs.  (C-la,  b)  become 

OT 

x "  +  A2tox2Jt  =  A  fieoy’  +  ^  fieo'y  (C-2a) 

y"  +  A  2to,2y  =  —  Afteox  ’  —  —  fioo’x.  (C-2b) 

Now  writing 

x  =  a,(T;A)eft/“(T)rfT'  (C-3a) 

y  =  a2(T;  A)ea/“<T',jT'  (C-3b) 


we  proceed  to  express  a,  and  a2  in  formal  asymp¬ 
totic  series 

(C-4a) 

n  —  0  A 

a2( t;  X)  ~  2  (C-4b) 

n-0  A 

We  must  now  find  the  ain,  a2n ,  and  to. 

Substituting  Eqs.  (C-4a,  b)  in  Eqs.  (C-2a,  b) 
one  finds  the  leading  order  (A2)  result 


(tox2  -  a>2)fl|0  -  /wfleo02o  =  0  (C-5a) 
iwfleoo  to  +  (to,2  —  (o2)fl2o  =  0  (C-5b) 


from  which  it  follows  that  o>  must  be  one  of  the 
four  quantities 

<v*2  +  a >y2  +  Ooo  y 

to  = 


±  [(U)X2  -f  Uly2 


+  flip)2  -  4o>xW]1'2  j  1/2 

(C-6) 


The  next  order  (A1)  relation  may  be  written,  after 
some  manipulation,  as 


,  ,f  ,  1  <02  —  OJy2  fliol 

+  '  w  ~  X - 7T- 

L  2  to  ileoj 

=  J^-fleo  +  (w2  —  to^ljaio 

r  i  o,-  (a,2  -  a,,2)] 

+  [—  2^6°  +  -  — ---  — J«.0. 


(C-7) 


i 


i 


33.15 


D.  CHERN1N  AND  P.  SPRANGLE 


Using  Eq.  (C-5a)  or  Eq.  (C-5b)  and  Eq.  (C-7)  an 
equation  for  just  a,0  (or  just  a20 )  may  be  obtained. 
The  solutions  are 

«io  =  Aia  ~  l/2(o>,2  -  to2)— 

(to/  +  <0/  +  nio  -  2co2)"  —  (C-8a) 

a2 o  =  Aa>-  ,/2(u>2  -  to/)— 

(to/  +  to/  +  Hlo  -  2to2)",/2  (C-8b) 


mogeneous  equations  as 

jr'-'A 

yu>)j=  I,  2,  3,  4  (C- 1 1 ) 

one  finds  in  a  straightforward  way  that  a  partic¬ 
ular  solution  to  Eqs.  (11a,  b)  in  the  text  is  given 
by 

(0'/WwD'v)'  (c-,2) 


where  A  is  an  arbitrary  complex  constant.  where 

Using  Eq.  (C-8a,  b)  and  the  definition  of  to  we 
may  write  the  leading  order  WKB  solution  to  K  . 
Eqs.  (C-la,  b)  as  1 


~  tjum  x('\t)x<k>(t')yu,«')y(m\r) 


[(to/  +  to/  +  n2#0)2  —  4<o/io/]  1/4 


(C-9) 


W  =  tjktm  x<J,(l)xik\t)y(,)(t)y(m\t), 


x 


(to/  -  to/)-\ 
(«/  -  to/)-; 


where  the  sum  extends  over  the  four  values  of 
to  in  Eq.  (C-6)  and  where  the  Aj  are  constants. 

This  solution  is  expected  to  be  valid  as  long  as 
to*,  to, ,  and  fl80  are  slowly  varying  compared  to 
any  t ojt  i.e.. 


I 


—  In  to 
dt 


<»j  |.  J  =  1,2,  3,  4, 


where  ib  is  to*,  to,,  or  fteo-  The  solution  is  there¬ 
fore  expected  to  fail  when  to,  =  0,  that  is,  near 
a  turning  point.  From  Eq.  (C-6)  this  can  happen 
when 


to/ to/  =  0.  (C-10) 

Equation  (C-9)  is  bounded,  however,  if  in  addi¬ 
tion  to/  =  to/  =  0.  Breakdown  of  Eq.  (C-9)  (and 
a  transition  to  unstable  behavior)  occurs  only  if 
to*2 to,2  =  0  and  to*2  ^  to,2. 

Once  the  solution  to  the  homogeneous  equa¬ 
tions  have  been  found  the  solution  to  the  inhom¬ 
ogeneous  equations  follows  by  the  usual  varia¬ 
tion  of  parameters  or  some  similar  method. 
Writing  four  independent  solutions  to  the  ho- 


and  where  the  summation  convention  is  under¬ 
stood.  The  Wronskian  W  is  a  constant,  indepen¬ 
dent  of  time;  its  value  is  determined  once  a  choice 
is  made  for  the  x(J),  y<J). 
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The  dynamics  of  an  ultrahigh-current  electron  ring  in  a  modified-betatron  configuration  is 
addressed.  This  study  includes  analytical  and  numerical  results  for  both  “cold”  and  “hot”  rings. 

It  has  been  found  that  the  walls  surrounding  the  ring  and  toroidal  effects  play  a  very  important 
role  in  the  dynamics  of  the  ring,  even  when  the  wall  conductivity  is  infinite.  For  finite  wall 
conductivity,  the  diffusion  of  the  ring’s  self-magnetic  field  profoundly  impacts  its  dynamics  and  in 
general  the  equilibrium  could  be  lost  if  means  are  not  provided  to  compensate  for  this  effect.  In 
addition,  it  has  been  found  that  the  toroidal-field  grad  B  drift  is  not  important,  except  when  the 
bounce  frequency  is  very  small.  The  general  conclusion  of  these  studies  is  that  equilibrium  states 
of  ultrahigh-current  rings  in  modified-betatron  configuration  exist  over  a  wide  range  of 
parameters.  These  states  are  accessible  and  within  the  reach  of  existing  pulsed  power  technology. 


I.  INTRODUCTION 

Over  the  last  few  years,  several  laboratories1-7  have 
been  engaged  in  studies  that  are  intended  to  assess  the  feasi¬ 
bility  of  developing  ultrahigh-current  accelerators.  These 
studies  are  mainly  motivated  by  the  potential  application  of 
high-energy,  high-current  electron  beams  in  nuclear  physics 
research,  medical  radiography,  and  the  fusion  program. 

In  general,  the  effort  is  not  directed  toward  developing 
novel  accelerating  schemes  but  rather  toward  modifying  or 
improving  the  existing  accelerator  technology.  Among  the 
various  proposed  modifications,  the  addition  of  a  strong  tor¬ 
oidal  magnetic  field  to  a  conventional  betatron8-9  has  attract¬ 
ed  considerable  attention,4-7  and  the  configuration  has  been 
named  the  modified  betatron. 

In  this  paper  we  analyze  and  discuss  the  dynamics  of  a 
high-current  electron  ring  confined  in  a  modified-betatron 
configuration.  When  the  intense  electron  ring  is  surrounded 
by  a  finite  conductivity  wall,  its  dynamics  can  be  divided, 
rather  naturally,  into  three  distinct  phases:  The  pre-accel¬ 
eration  phase  follows  the  injection  and  trapping  of  the  elec¬ 
tron  beam.  The  duration  of  this  phase  is  short,  typically  a  few 
microseconds  and  the  conducting  wall  surrounding  the  ring 
can  be  treated  as  a  perfect  conductor.  The  pre-acceleration 
phase  is  followed  by  the  diffusion  phase,  which  lasts  for  a 
time  that  is  of  the  order  of  the  magnetic  field  diffusion  time. 
During  the  diffusion  phase  the  self-magnetic  field  of  the 
beam  diffuses  out  of  the  metal  torus.  As  a  result  the  electron 
energy  is  reduced  and  the  net  radial  force  on  the  ring  in¬ 
creases,  but  at  a  different  rate.  Therefore,  the  equilibrium 
will  be  lost  if  means  are  not  provided  to  balance  these  two 
effects.  The  last  phase,  i.e.,  the  main  acceleration,  starts  after 
the  self-magnetic  field  diffuses  out  of  the  torus  and  has  a 
duration  that  is  comparable  with  the  acceleration  time.  For 
most  of  the  third  phase  the  energy  of  the  ring  has  increased 
substantially  and  thus  the  effect  of  the  self-fields  is  apprecia¬ 
bly  diminished. 

•'Permanent  address:  Berkeley  Research  Associates,  Springfield,  Virginia 
22 1  SO. 

**  Permanent  address:  Sacns, .  ;  eeman  Associates,  Bowie,  Maryland  207 1 5. 


The  present  work  addresses  the  dynamics  of  the  high- 
current  ring  in  a  modified-betatron  configuration.  It  in¬ 
cludes  both  analytical  and  computational  studies  for  “cold” 
as  well  as  “hot”  electron  rings.  The  main  conclusions  are: 
First,  if  the  energy  of  the  injected  electrons  is  not  exactly 
equal  to  the  energy  corresponding  to  the  equilibrium  orbit, 
which  is  assumed  to  coincide  with  the  minor  axis  of  the 
torus,  the  center  of  the  orbit  is  displaced  from  the  center  of 
the  minor  cross  section  of  the  torus.  The  displacement  is 
proportional  to  the  energy  mismatch.  This  imposes  very 
stringent  constraints  on  the  injector.  Second,  during  the  dif¬ 
fusion  of  the  self-magnetic  field  of  the  ring  out  of  the 
chamber,  the  equilibrium  can  be  lost  if  means  are  not  pro¬ 
vided  to  balance  the  change  of  the  ring  equilibrium  radius. 
This  change  is  due  to  the  reduction  of  the  relativistic  factor 
y0  and  the  increased  forces  acting  on  the  center  of  the  ring. 
Third,  the  modification  in  the  flux  rule  for  high-current 
rings  can  be  ignored,  provided  that  v/y£  (v  —  Budker  pa¬ 
rameter)  is  very  small.  Fourth,  with  the  exception  of  the  ex¬ 
pansion  of  the  minor  cross  section  of  the  ring,  finite  emit- 
tance  does  not  have  any  other  noticeable  effect  on  the 
equilibrium,  and  fifth,  considerably  higher  axial  energy 
spread  can  be  tolerated  in  high-  than  in  low-current  rings. 

The  general  and  most  important  conclusion  of  these 
studies  is  that  equilibrium  states  of  high-current  rings  in  a 
modified-betatron  configuration  exist  over  a  wide  range  of 
parameters.  These  equilibria  are  realistic  and  accessible  with 
existing  technology. 


II.  ORBIT  STABILITY  OF  A  COLD  RING  WITHOUT 
SURROUNDING  WALLS  AND  TOROIDAL 
CORRECTIONS 

The  stability  of  a  single  particle  orbit,  when  thermal 
effects,  wall  effects,  and  toroidal  corrections  are  neglected, 
has  been  considered  previously.4  Using  the  coordinate  sys¬ 
tem  shown  in  Fig.  I  and  assuming  that  the  external-field 
components  vary  as 

)sxB0,{t  )[1  —  w(r  —  r0)/r0],  .  (la) 
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FIG.  1.  Schematic  of  the  modified -betatron  configuration  and  system  of 
coordinates  used  in  the  analysis. 


B,[r,t)&  - B0z(t)nz/r0, 

(lb) 

Ee(r,t )  —  -  -  fVr/r'-^MrV), 
rc  Jo  ot 

(lc) 

and 

B0  M  [  1  (r  r0)/r0], 

(Id) 

it  has  been  determined  that  for  Bg>Bz  and  v/y0<l  the  orbits 
are  stable,  provided  that  the  following  condition  is  satisfied: 

nt  <{Boo/EB0z  )2. 

(2) 

In  the  above  equations  BI(r,t )  is  the  axial  and  Br[r,t )  is  the 
radial  component  of  the  betatron  field,  B0[r)  is  the  toroidal 
magnetic  field,  E0{r,t )  is  the  induced  electric  field,  n  is  the 
external-field  index,  (ro,0)  are  the  coordinates  of  the  center  of 
the  electron  ring,  and  n,  is  the  self-field  index.  The  self-field 
index  is  defined  as 

n,  =  <ol/(2rof2 1 , )  =  2  (v/y0)[c/f20l  rb  )\  (3) 

where  <y*(  =  Ane^n^/m)  is  the  beam  plasma  frequency 
squared,  y0  is  the  relativistic  factor,  D0z(  =  eB0z/mc )  is  the 
cyclotron  frequency  corresponding  to  the  axial  component 
of  the  betatron  field,  rb  is  the  minor  radius  of  the  beam,  and  v 
is  Budker’s  parameter,  i.e.,  the  product  of  the  number  of 
electrons  per  unit  length  times  the  electron  classical  radius. 

The  electron-beam  current  that  can  be  stably  confined 
in  a  modified  betatron  can  be  obtained  by  substituting  Eq.  (3) 
into  Eq.  (2)  and  is 

Imb  <2.1  [rb/r0)Y0  Pl[B00/B0l  f  (kA)  (4) 

where P0  =  Vq/c.  If  y0>  1,  then  rjl0z/y0^c  and  Eq.  (4)  be¬ 
comes 

/mt<7.22XlO-7  ri YoBle  (kA), 

with  rb  in  cm  and  B00  in  G. 

In  addition,  it  has  been  shown  under  the  same  condi¬ 
tions,  but  with  Be  =  0  (conventional  betatron),  the  beam  is 
stable  provided 

n,  <  J,  conventional  betatron  [Bg  =  0).  (5) 

The  electron-beam  current  that  can  be  stably  confined  in  a 
conventional  betatron  can  be  obtained  by  substituting  Eq.  (3) 
into  Eq.  (3)  and  is 

/rt  <4.2(rb/r0)Vo  Pi  (kA),  (6) 

or 

/rt<  14.44  X  lO"1  rlyoB l,  (kA), 


when  y0>  1.  The  ratio  of 7„lfc  to/r,,  obtained  from  Eqs.  (4)and 
(6)  is 

I mb  /I cb  ~  J  {Boo/B0z)‘. 

The  above  relation  indicates  that  for  Bu0  >B(„ ,  the  electron- 
beam  current  that  can  be  stably  confined  in  a  modified  beta¬ 
tron  substantially  exceeds  the  current  that  can  be  confined  in 
a  conventional  betatron. 

The  modified-betatron  stability  condition  given  in  Eq. 
(2)  can  be  easily  obtained  from  the  well-known  confined 
equilibrium  condition10 

2(<uJ/72  7)<1  (nonrelativistic),  (7) 

where  a>b  is  the  beam  plasma  frequency  and  (lc  is  the  cyclo¬ 
tron  frequency.  For  relativistic  energies,  and  taking  into  ac¬ 
count  the  self-magnetic  field  of  the  beam,  (o7b  becomes 

*>\-+  K2/r0)(l/^),  (8) 

where  (l/j'o)  accounts  for  the  self-magnetic  field  of  the 
beam,  and 

n  l-+ni/ri.  (9) 

Substituting  Eqs.  (8)  and  (9)  into  Eq.  (7)  we  obtain  the  stabil¬ 
ity  condition  of  Eq.  (2).  The  equality  sign  in  Eq.  (2)  gives  the 
maximum  electron  density  that  can  be  supported  at  a  specif¬ 
ic  value  of  Boe/B0z  and  the  corresponding  equilibrium  is 
known  as  Brillouin  flow.10 

III.  WALL  EFFECTS  ON  THE  MACROSCOPIC  MOTION 
OF  A  COLD  BEAM 

In  this  section  we  analyze  the  effects  of  surrounding 
walls  on  the  motion  of  the  center  of  the  beam.  In  the  Sec. 
Ill  A,  it  is  assumed  that  the  perfect  conductor  that  sur¬ 
rounds  the  beam  is  a  straight  cyclindrical  pipe  of  circular 
cross  section  and  thus  toroidal  effects  (hoop  forces)  associat¬ 
ed  with  the  fields  are  omitted.  These  forces  act  to  expand  the 
ring  and  are  reduced  with  increasing  ring  major  radius.  The 
toroidal  effects  are  discussed  in  Sec.  Ill  B. 

A.  Without  toroidal  corrections 

Consider  a  pencil-like  electron  beam  inside  a  straight, 
perfectly  conducting  cylindrical  pipe  of  circular  cross  sec¬ 
tion  as  shown  in  Fig.  2.  The  center  of  the  beam  is  located  at  a 

z 


FIG.  2.  Watt  (images)  forces  acting  on  a  pencil-like  electron  beam,  situated 
inside  a  perfectly  conducting  cylindrical  pipe. 
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distance  Jr,  Az  from  the  center  of  the  minor  cross  section  of 
the  pipe.  As  a  result  of  the  induced  charges  on  the  wall,  the 
center  of  the  beam  will  experience  a  radial,  outward-directed 
force,  which  for  small  displacements,  i.e.,  Ar,  Az<,a  is  given 
by 

F£  =  2i Te2n0(rb/a)'2{Are,  +  Aze,),  (10) 

where  a  is  the  cylinder  radius  and  n0  the  uniform  beam  den¬ 
sity. 

Similarly,  as  a  result  of  the  induced  current  on  the  wall, 
the  center  of  the  beam  will  experience  a  radial  force  that  is 
directed  toward  the  opposite  direction  than  F£  and  is  given 
by 

TB=~PlFE.  (11) 

Using  the  external  fields  of  Eq.  (1)  and  the  induced  fields  of 
Eqs.  (10)  and  (11),  the  equations  describing  the  temporal  lin¬ 
ear  evolution  of  the  beam's  center,  for  time-independent  ap¬ 
plied  fields,  are 

Ar  4-  eb)Ar  =  ~^~Az  +  ^8Psl  ,  (12) 


(00,' 

j  (SPe) 

v  Yo  > 

'  Yo^ro 

Az  +  S7Az  =  ~{nog/y0]Ar,  (13) 

where 

5’  =  (^oz/ro)J(l  -  n  -  nSb/a7), 

S’  =  (rt0,/y0)J(n  -  n.d/a7), 

nog  =  eBog  /  me,  n0l  =  eB0z  /  me, 

SYtjYo  =  Po(SPg 

and  SPg  is  the  difference  between  the  canonical  angular  mo¬ 
mentum  of  an  electron  at  (r,z)  and  its  corresponding  value  at 
the  equilibrium  orbit  (ro,0).  The  average  is  over  initial  coordi¬ 
nates  and  velocities.  Equations  (12)  and  (13)  do  not  include 
the  self-electric  and  self-magnetic  fields,  because  both  these 
fields  are  zero  at  the  center  of  a  straight  beam.  In  addition, 
the  nonlinear  terms  (Ar/r0)(fl0g/y0)Az  and  [Ar/r0)[flog/ 
y0)Ar  have  been  omitted  from  Eqs.  (12)  and  (13).  These  two 
terms  have  their  origin  in  the  gradient  of  the  toroidal  mag¬ 
netic  field  and  are  considered  in  Sec.  V.  In  general  these 
terms  are  not  significant  except  in  the  limit  5,— *0. 

In  Eq.  (12),  Sy0=  P0(6Pg)/mr,/:  indicates  the  energy 
mismatch,  i.e.,  the  difference  between  the  energy  of  the  refer¬ 
ence  electron  (moving  along  the  axis  of  the  beam)  and  the 
energy  required  for  the  same  electron  to  move  on  the  equilib¬ 
rium  orbit  (r0,0).  The  solution  of  Eqs.  (12)  and  (13),  for  time- 
independent  fields,  is 

+  X  0(5*  - «?)' V*  (14) 

Y<i(0rrOm  J~  ' 

and 

jz=  (is) 

J-  i 

where  Cj  is  a  constant  and 

u]  -  +  S’  +  {Oo0/y0)2 

±  |[5’  +5’  +  [f2og/y0)7]7  -  45’SJ ) ,n).  (16) 


FIG.  3.  Bounce  frequency  squared  as  a  function  of  n,pt/a>,  with  the  field 
index  n  as  a  parameter.  The  orbits  are  open  when  (if],  <0  and  closed  when 

*i>o. 

The  first  term  on  the  right-hand  side  (RHS)  of  Eq.  (14) 
gives  the  displacement  of  the  center  of  the  orbit  from  the 
center  of  the  surrounding  cylindrical  pipe  and  can  be  written 
as 

Ar0=  - Wfj)/? ° - .  (17) 

o«  /Yo)i  1  -  n  -  n ,  rj  /a7)mr0 

Equation  (16)  has  four  roots.  Two  of  them  are  fast  (plus  sign) 
and  two  are  slow  (minus  sign).  When  B„  '>B1 ,  the  slow  modes 
become 

J  \O0B/y0) 

and  the  fast  modes  are  ~  ±  (Ahj/Vo)-  Equation  (18)  is  plot¬ 
ted  in  Fig.  3  for  two  values  of  the  external-field  index.  For 
n  >  J,  <o),  is  negative  when  1  —  n  <  n,  Pb/a7  <  n  and  for  n  <  J, 
a>l  is  negative  when/j<n,ri/a}<l  —  it.  Negative  values  of 
indicate  that  the  beam  motion  is  unstable  and  the  orbit  in 
the  rj  plane  is  open.  Since  the  parameter  n,  rj  /a7  scale  as 
Yo},  during  acceleration  n,ii/a7  decreases  rapidly.  There¬ 
fore,  in  order  to  avoid  the  instability,  it  is  necessary  that 
before  the  commencement  of  the  acceleration  the  parameter 
n,  rj,  /a7  <  1  —  n  when  n  >  $  and  <  n  when  n  <  J.  This  implies 
that  the  injected  beam  current  should  be  limited  to 
/<8.5(1  -  n]Ply-„a7/ri  (kA)  forn>J, 
and 

/<8.5n/?o7o0’/'o  (^A)  for  /><}. 

The  orbit  of  the  ring’s  center  is  described  by 
Aiii )  =  [dr(0)  -  dr0J  cos  coBt 

+  [Jr(0)/(uB]  sin  coBt  +  Ar0, 

Az{t )  =  Az(0)  cos (<y„r )  +  [Az{0)/cob]  sin(wfi/), 
where  the  initial  velocities  and  displacements  are  related  by 
M 0)  =  ~  [Z>l/[noa/y0)]Azl 0), 


Az[0)  =  — — Ar\  0)-(  A 

Woo/Yo)  ‘  \n09) 


(  ^op^Yo) 


’tto.) 

(SPe) 

■P^oe  ' 

Yo^i-o 

Ar0 ), 
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FIG.  4.  Snapshots  of  the  electron-ring  minor  cross  section  at  four  different 
times.  The  values  of  the  various  parameters  are  listed  in  Table  I.  The  center 
of  the  ring’s  minor  cross  section  describes  a  circle. 
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FIG.  5.  An  electron  ring  inside  a  conducting  toroidal  chamber.  Because  of 
the  toroidal  effects,  neither  the  self-electric  nor  the  self-magnetic  field  is 
zero  at  the  geometric  center  of  the  ring. 
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(SPg)  =  Y<fnr0{vg  -  rJ20,/y0). 

The  predictions  of  Eq.  (18)  are  in  good  agreement  with 
the  results  of  computer  simulation  shown  in  Fig.  4.  This 
figure  shows  four  snapshots  of  the  beam  in  the  plane.  At 
t  =  0  the  circular  electron  beam  is  injected  near  the  center  of 
the  pipe.  The  values  of  the  various  parameters  are  listed  in 
Table  I.  The  electron-beam  current  is  kept  low  (1  kA)  in 
order  to  minimize  the  toroidal  effects.  As  can  be  seen  from 
Fig.  4,  the  center  of  the  electron  beam  describes  a  circle  in  the 
rj  plane  with  a  period  of  about  188  nsec.  For  the  same  pa¬ 
rameters  Eq.  (18)  predicts  a  period  rB  =  2jr/<a^~180  nsec. 
These  numerical  results  are  discussed  further  in  the  Sec. 
IIIB. 

The  displacement  of  the  orbit’s  center  because  of  the 
energy  mismatch  [Eq.  (17)],  imposes  very  stringent  con¬ 
straints  on  the  injector.  This  becomes  apparent  when  we 
consider  some  limiting  cases.  For  example,  when  n  =  }  and 
ntr*k/a2<l,  Eq.  (17)  is  reduced  to 

Ar(/rQz;2(8y(/yQ).  (19) 

Equation  (19)  predicts  that  for  a  major  radius  r0  =  100  cm, 
the  ratio  8yo/y0  should  be  less  than  1  %  in  order  that  the 
displacement  of  the  orbit  will  be  less  than  2  cm.  The  condi¬ 
tion  8y,/y0<l%  requires  that  the  uncertainty  in  energy 
should  be  less  than  35  keV,  when  the  energy  of  the  injected 

TABLE  I.  Electron-ring  parameters  for  the  results  of  Fig.  4. 

Beam  energy  (MeV)  3 

Beam  current  (kA)  * 

Beam  minor  radius  (cm) 

Beam  major  radius  (cm)  50 

Torus  minor  radius  (cm)  6.4 

Vertical  magnetic  field  (G)  240 

Toroidal  magnetic  field  (kG)  2 

External-field  index  0-42 

Self-field  index 
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beam  is  3  MeV.  Although  such  a  small  energy  uncertainty 
can  be  attained  with  sophisticated  injectors,  there  are  other 
factors,  such  as  space  charge  and  inductive  effects,  which 
contribute  substantially  to  the  uncertainty  of  energy.  Al¬ 
though  the  displacement  cannot  be  eliminated,  its  negative 
impact  can  be  alleviated  by  elongating  the  minor  cross  sec¬ 
tion  of  the  torus  along  the  r  direction. 

B.  Toroidal  corrections5,11 

The  cause  of  these  effects  is  the  finite  curvature  of  the 
electron-beam  orbit.  For  relatively  large  aspect  ratio  (r,/ 
rb  >  1 )  beams,  the  toroidal  effects  become  important  when  v/ 
y0  exceeds  a  few  percent.  Previous  work12  on  toroidal  effects 
was  limited  to  "bare”  rings,  i.e.,  without  surrounding  con¬ 
ducting  walls.  In  this  subsection  we  briefly  analyze  the  toroi¬ 
dal  corrections  for  the  more  realistic  geometry  shown  in  Fig. 
5,  i.e.,  including  the  effect  of  conducting  wall  around  the 
beam.  Since  different  physics  issues  are  involved  in  the  case 
of  a  “bare”  and  a  "shielded”  ring,  it  is  appropriate  to  start 
our  discussion  with  a  bare  ring. 

Consider  an  intense  electron  ring  as  shown  in  Fig.  5,  but 
without  the  toroidal  chamber.  At  the  inner  edge  of  the  RHS 
cross  section  of  the  ring  the  self-magnetic  field  is  greater  than 
that  of  a  straight  beam  with  the  same  parameters,  because  of 
the  contribution  from  the  left-hand  side  (LHS)  of  the  ring.  At 
the  outer  edge  of  the  cross  section,  the  self-magnetic  field  is 
reduced  because  the  contribution  from  the  LHS  of  the  ring 
has  different  polarity  than  the  local  field.  Thus,  the  total  self- 
magnetic  field  is  shifted  upwards,  as  shown  schematically  in 
Fig.  5.  In  contrast,  the  self-electric  field  decreases  at  the  in¬ 
ner  radius  and  its  magnitude  increases  at  the  outer  radius  of 
the  ring.  Thus  the  self-electric  field  shifts  downwards.  As  a 
result  of  these  shifts  neither  the  magnetic  nor  the  electric 
field  are  zero  at  the  geometric  center  of  the  ring. 

When  the  electron  ring  is  surrounded  by  a  perfect  con¬ 
ductor,  the  shift  in  the  fields  discussed  above  is  reduced  but 
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additional  field  components  appear  as  a  result  of  induced 
charge  and  current  on  the  conducting  wall.  For  small  ring 
displacements  from  the  center  of  the  minor  cross  section  of 
the  torus,  the  induced  fields  vary  linearly  with  the  displace¬ 
ment  and  arc  identical  to  those  given  in  Eqs.  (10)  and  (11),  for 
a  cylindrical  pipe.  Quantitative  extension  of  the  previous 
qualitative  considerations  has  shown  that  the  fields  at  the 
center  of  a  uniform-density  electron  ring  inside  a  perfectly 
conducting  toroidal  chamber  of  circular  cross  section  are4,1 1 

EU--2,Wvo{[^  +  |iln(i)]#, 


d  Az  „ 

2  * 
a  r0 


Bind  =  —  2j7’|e|«0/?0r0|^-  e, 

-[?T-3(,+7h?)]4  ,211 

where  n0  is  the  ambient  density,  0O  =  vjc,  v0  is  the  azi¬ 
muthal  velocity  defined  by 

„  _  _ rtfloi/Yo _  Q21 

0  1  +  2(v/y0)[  1  +  \n(a/rb))  ' 

and  the  displacement  dr,  Az  of  the  ring  from  the  center  of  the 
torus  has  been  assumed  to  be  much  less  than  a.  In  addition, 
the  fields  given  by  Eqs.  (20)  and  (2 1 )  have  been  derived  under 
the  assumption  that  the  angular  frequency  of  the  electrons  is 
constant  and  therefore  the  electron  current  density  varies 
proportionally  to  r. 

Using  Eqs.  (20)-(22),  it  can  be  shown  that  the  center  of 
the  beam  is  described  by  Eqs.  (12)  and  (13)  with  to’  and  S’ 
replaced  by 

5j—  {fio,/ro)i<x  -n*-n, rj/o2),  (23) 

iJYokn*  -  n, d/a\  (24) 

and 

SPe—*iSPg, 

where 

«  =  £’=[  1  +  (2v/y0)[I  +  \n[a/rb\\\ (25) 


n*  =  n£. 

The  bounce  frequency  can  be  found  by  substituting  Eqs. 
(23)  and  (24)  into  Eq.  (18)  and  is 


f^A2  ( 

\B0J  \ 

<  Yo  /  V 

When  n *  =  a/2,  for  the  reasons  stated  in  the  paragraphs 
following  Eq.  (18),  the  orbits  are  closed  (stable)  as  y0  in¬ 
creases,  provided 

n,r3b/a1<a/2 

or, 


/<4.25(--)Yl  —  —  In -Vo  (kA)  for  r<>l.  (27) 

W  V  Yo  rj 

The  limiting  current  given  by  Eq.  (27)  is  considerably 
lower  than  that  of  Eq.  (4)  when  B00>B0l.  In  addition,  in 
contrast  to  the  current  given  by  Eq.  (4),  the  current  of  Eq. 
(27)  is  independent  of  the  toroidal  magnetic  field  B0$.  Toroi¬ 
dal  effects  shift  the  curves  of  Fig.  3  to  the  left  and  the  critical 
current  given  by  Eq.  (27)  is  lower  than  the  corresponding 
critical  current  when  toroidal  effects  are  omitted. 

At  this  point,  it  is  appropriate  to  return  and  discuss 
further  the  computer  simulation  results  of  Fig.  4.  For  the  1 
kA,  3  MeV  beam  the  v/y0  =  0.0084.  Even  at  this  small  value 
of  v/y0  toroidal  effects  are  noticeable.  The  value  of  £  com¬ 
puted  from  Eq.  (25)  using  the  above  value  of  v/y0  and  6.4  for 
the  ratio  a/rb  is  0.95.  According  to  Eqs.  (23)  and  (24)  the 
center  of  the  beam  will  describe  a  circle  when  a  —  n*  =  n* 
or  for  n  =  J[1  —  (v/y0)  ln(a/rb  )]£  =  0.467,  which  is  in  excel¬ 
lent  agreement  with  the  simulation.  In  addition,  Eq.  (26)  pre¬ 
dicts  a  period  of  1 88  nsec,  which  also  is  in  excellent  agree¬ 
ment  with  the  simulation. 

The  most  striking  manifestation  of  toroidal  effects  is  in 
the  value  of  the  betatron  magnetic  field  required  to  confine 
the  rotating  beam  at  a  specific  radius.  When  the  axis  of  the 
beam  lies  along  the  axis  of  the  torus,  i.e.,  when  Ar  —  Az  —  0, 
it  can  be  shown  from  Eqs.  (20)  and  (21)  that  the  external 
magnetic  field  required  for  the  beam  to  rotate  with  a  radius 
r0is 

K  =B0\\+[2V/y0)[\+\n[a/rb))\, 

where  B0  is  the  magnetic  field  necessary  for  a  single  particle 
of  the  same  energy  to  rotate  with  a  radius  r0. 

The  above  expression  for  the  magnetic  field  is  based  on 
the  assumption  that  all  the  electrons  rotate  with  a  constant 
angular  frequency,  i.e.,  the  current  density  increases  linearly 
across  the  beam.  If  the  current  density  is  constant  across  the 
beam,  the  above  expression  is  slightly  modified  and  becomes 

B0,  =  30\  l  +  (2v/yo)[0.i  +  ln(a/rb)]). 

For  a  10  kA,  2  MeV  uniform  current  density  beam  with 
a  ratio  a/rb  =  6.4,  the  correction  is  55%,  i.e.,  B0z/ 
B0  =  1.55.  This  effect  is  demonstrated  clearly  by  the  results 
of  computer  simulation  shown  in  Fig.  6.  The  three  snapshots 
of  the  electron-ring  minor  cross  section  in  a  modified-beta- 
tron  field  correspond  to  f~0, 20,  and  40  nsec.  For  all  practi¬ 
cal  purposes,  the  minor  radius  of  the  beam  remains  constant. 
The  external  betatron  magnetic  field  is  127  G,  i.e.,  approxi- 


R(cm>  RUm)  R(cm) 


FIG.  6.  Three  snapshots  of  the  electron-ring  minor  cross  section  in  a  modi- 
fied-betatron  field.  The  external  magnetic  field  required  to  confine  the  ring 
is  50%  higher  than  the  single-particle  field.  The  values  of  the  various  pa¬ 
rameters  in  this  run  are  listed  in  Table  11. 
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TABLE  II.  Parameters  of  the  uniform  current  density  ring  for  the  results 
shown  in  Fig.  6. 


Beam  energy  (MeV)  2 

Beam  current  (kA)  10 

Beam  minor  radius  (cm)  I 

Beam  major  radius  (cm)  100 

Torus  minor  radius  (cm)  6.4 

Vertical  magnetic  field  (G)  127 

Toroidal  magnetic  field  (kG)  1.7 

External-field  index  0.25 


mately  50%  higher  than  the  single  particle  field.  The  re¬ 
maining  parameters  are  summarized  in  Table  II. 

Another  consequence  of  the  toroidal  effects  is  the  in¬ 
crease  in  the  radial  displacement  of  the  orbit  for  fixed  energy 
mismatch  and  major  radius.  For  n,rj/a7<  1,  the  orbit  dis¬ 
placement  becomes 

Ar0/r0~r0[5r0/(y20  -  1)]/(1  -«/£). 

In  addition,  th  e<u*  vs  n,p-b/a2  curves  of  Fig.  3  are  shifted  to 
the  left  and  thus  the  maximum  permissible  current  that  can 
be  accelerated  is  reduced,  as  may  be  seen  from  Eq.  (27). 

IV.  ACCELERATION 

After  injection  and  trapping  the  beam  is  accelerated  to 
high  energy  by  the  inductive  electric  field  generated  by  the 
time  varying  betatron  magnetic  field.  As  shown  in  Fig.  7,  the 
acceleration  may  be  divided  into  three  phases:(l)  pre-accel¬ 
eration,  (2)  diffusion,  and  (3)  main  acceleration. 

The  pre-acceleration  phase  occurs  for  times  much 
shorter  than  the  field-diffusion  time,  i.e.,  l<tD  =  (An/ 
c*)a8a  In (tg/a),  where  a  is  the  conductivity,  S  is  the  thick¬ 
ness,  a  is  the  minor  radius,  and  r0  is  the  major  radius  of  the 
torus.  During  the  pre-acceleration  phase  the  self-magnetic 
field  of  the  beam  does  not  have  time  to  diffuse  out  of  the 
conducting  torus.  In  the  example  shown  in  Fig.  7  the  total 
acceleration  time  was  chosen  to  be  1  msec  and  the  diffusion 
time  10/rsec.  The  ratio  of  the  temporal  extent  of  each  phase 
to  the  total  acceleration  time  t  /t„  is  also  given  in  the  figure.  It 
is  apparent  that  the  beam  spends  most  of  its  time  in  the  main 
acceleration  phase. 

During  the  pre-acceleration  phase  the  metal  wall  sur¬ 
rounding  the  beam  can  be  treated  as  a  perfect  conductor.  As 
a  result,  for  small  displacement  of  the  beam,  the  self-flux 


fiCMtofBtloil  PlWlH 


FIG.  7.  The  three  phases  of  acceleration.  In  the  example  shown  the  accel¬ 
eration  time  is  I  msec  and  the  diffusion  time  is  10/iiec. 
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linking  the  axis  of  the  electron  ring  remains  constant  to  lo 
est  order.  This  point  is  discussed  further  later. 

The  condition  for  the  major  radius  of  the  electron  ring 
to  remain  constant  can  be  derived  from  the  conservation  of 
canonical  angular  momentum 

Pe  =  my,/vog  -  (e/c)r(A  -f  A  ■„),  (28) 

where  voe  is  the  toroidal  velocity,  A  is  the  vector  potential 
that  describes  the  external-magnetic  (betatron)  field,  A’e  is 
the  vector  potential  that  describes  the  self-magnetic  field  of 
the  beam  and  r  is  the  instantaneous  radius  of  the  beam. 

The  self-flux  linking  the  axis  of  the  ring  <P*  is 

4>‘  =  2  nrA’0,  (29) 

and  to  lowest  order  is  independent  of  the  beam  displacement 
and  thus  of  time  [see  Eq.  (35)].  Substituting  Eq.  (29)  into  Eq. 
(28)  and  using  the  expression  for 

v  _  _ rnQt/y0 _ 

06  1  +  2(v/ y0) [  1  +  \n(a/rb)]  ’ 

Eq.  (28)  predicts  that  r  remains  constant,  provided 

—  f - - )  =  — —  <2?"*).  (30) 

Ml+2(v/r0)[l+ln(o//*)J/  2  dl 

In  Eq.  (30),  B0l  (t )  is  the  local  magnetic  field,  v  is  the 
Budker’s  parameter,  a  is  the  minor  radius  of  the  torus,  rb  is 
the  beam  radius,  and  (/T*' )  is  the  average  applied  magnetic 
field  within  the  electron  ring. 

Using  the  equation 

mc2^  =  -evE,  (31a) 

dl 

together  with 

£*=-  <*"’>.  (31b) 

2c  at 

Eq.  (30)  becomes 

^J  =  l[,  +  f(,+'^)]l<*'“>-  1321 

A  similar  condition  has  been  derived  previously.13 

Equation  (32)  is  the  condition  that  must  be  satisfied  in 
order  for  the  radius  of  the  accelerated  electron  ring  to  remain 
constant.  For  low-v/y0  beams,  Eq.  (32)  is  reduced  lo  the 
well-known  flux  rule,  i.e., 

dB0r  ( t  )  _  d_  /  n  CIV 

dl  2  dl 

The  correction  term  (2v/yj])[l  -f  \n(a/rh)]  in  Eq.  (32)  is 
very  sensitive  to  the  beam  energy.  For  a  10  kA,  3  MeV  beam 
injected  into  a  10cm  minor-radius  torus  with  rb  =  1  cm,  the 
correction  is  only  1  %  and  therefore  can  be  neglected.  How¬ 
ever,  when  the  energy  of  the  same  beam  is  reduced  to  0.5 
MeV  the  correction  is  48%,  i.e.,  very  substantia).  The  design 
of  the  accelerator  is  simplified  considerably  by  choosing  the 
beam  parameters  such  that  the  correction  term  is  negligible. 

The  instantaneous  value  of  y^t )  may  be  determined 
from  Eqs.  (30)  and  (3 1 )  and  for  y£(r  )>  1  is 

yjjt)  _  ( 1  +  2[ v/y0(Q)l  [  1  +  ln(fl/r„)] )  B0,  (/ ) 
yo(0)  [  1  +  2[v/y0(t ))  [  1  +  In (a/rb)) )  ^(0) ' 
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FIG.  8.  System  of  coordinates:  (a)  depicts  the  various  parameters  appearing  in  Eq.  (34);  and  (b)  depicts  a  displaced  straight  electron  beam  inside  a  perfectly 
conducting  pipe.  The  image  current  is  located  at  a  distance  a’/l^r5  +  Az*)'n  from  the  center  of  the  pipe. 


In  contrast  to  the  low-current  beams,  in  which  y0(t  }/B0l  (t ) 
remains  approximately  constant  during  acceleration,  for 
high-current  beams  with  toroidal  corrections  the  quantity 
that  remains  constant  during  acceleration  is 


Jjslil  i  4.  2  — —  [  1  -f  ln^— }]]  =  const. 

5or(Ol  7o(')l  Vfc/JJ 


As  y0(t )  increases  with  lime,  the  correction  term  decreases 
and  Eq.  (33)  is  reduced  to  y0(t)/B0z[t )  =  const. 

Now,  we  return  to  discuss  the  self-flux  linking  the  axis 
of  the  ring.  In  the  notation  of  the  toroidal  geometry  shown  in 
Fig.  8(a),  the  magnetic  stream  function  ¥  =  rA‘g  for/?  <  rb  is 
given,  to  lowest  order,  by" 


(34) 


where  Q  —  —  \e\n0irr],OR  l/c,  12  —  v9/r,  A  2  =  Ar3  -f  Az3, 
rt0  is  the  uniform  density,  and  it  has  been  assumed  that  o<r0, 
rb  <a,  and  Ar^Az^a. 

The  self-flux  through  the  axis  of  the  ring  ( p  —  0)  is 


<t>'  =  2jt!P.  (35) 

Substituting  Eq.  (34)  into  Eq.  (35),  we  get 

=27r0[l  +  21n(a/ri,)],  (36) 

which  does  not  depend  on  the  displacement  Ar,  Az  of  the 
beam.  This  is  a  rather  unusual  result  and  deserves  further 
discussion. 

To  gain  some  insight  into  the  problem,  we  have  com¬ 
puted  the  flux  linking  a  horizontal  surface  s  extended  from 
the  axis  of  a  straight  beam  to  the  inner  wall  of  a  perfectly 
conducting  cylinder  of  circular  cross  section  as  shown  in  Fig. 
8(b).  The  beam  displacement  is  arbitrary  but  the  ratio  rb/ 
a<  1 .  It  is  straightforward  to  show  that  the  flux  through  such 
a  surface  of  length  /  is 


^)] 


(37) 


where /is  the  beam  current 

For  (dr2  +  Az2)in/a<l,  Eq.  (37j  indicates  that  the  flux 
has  a  quadratic  dependence  on  the  beam  displacement. 
Therefore,  since  only  the  linear  terms  on  beam  displacement 
were  kept  in  the  derivation  of  Eq.  (34),  it  is  not  surprising  that 
the  flux  given  by  Eq.  (36)  is  independent  of  the  displacement 
of  the  beam. 

When  the  time  approaches  the  magnetic-field-penetra¬ 
tion  time  tD,  the  self-magnetic  field  of  the  electron  beam 
starts  to  diffuse  out  of  the  finite  conductivity  metal  torus. 
Using  the  geometry  shown  in  Fig.  9,  it  is  shown  in  Appendix 
A,  that  for  a  very  thin  conductor,  i.e.,  o~b,  to  lowest  order, 
the  self-magnetic  field  of  the  beam  at  time  t  is 


Bb[r,t)~[2I/cr)(l  —e~,/'D),  r>b, 
B*  [r,t ) a:  21  /cr,  rb  <r<a, 


FIG,  9.  System  of  coordinates  used  to  compute  the  magnetic  field  during 
diffusion. 
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where 


lD  ^{4iro5a/c')  lnj r^/a). 


These  equations  are  based  on  the  assumptions  that  S  / 
<i<1,  [(5/6)ln(r0/6)]l,2<!>  and  the  inductive  electric  field 
varies  as  Et  (r,t )  =  {b  /c)[dB4  (b,t  )/dt  ]  ln(r/r0),  i.e.,  goes  to 
zero  at  r  —  r0  (pseudotoroidal  geometry). 

The  self-magnetic  flux  in  the  three  regions  is  given  by 


f  (4nro/ c]I\n[a/rb ), 


<P'  ~ 


0. 


\{Avrf/c)I[\  -e-,/'D)ln{r0/b), 


r<a, 

a<r<b, 

r>b. 


and  it  diffuses  out  at  the  rate 


iY  (39) 

dt  tDc  \  b  / 


The  inductive  electric  field  generated  by  the  changing 
flux  given  in  Eq.  (39)  acts  to  slow  down  the  beam.  In  addi¬ 
tion,  for  a  constant  current  density  ring,  the  hoop  forces 
increase  by  the  ratio  [1  +  \n[ir(/rb  )]/[l  +  In (a/rb )],  and  the 
induced  magnetic-field  components  go  to  zero  at  the  end  of 
the  diffusion.  However,  the  electric-field  components  re¬ 
main  the  same.  As  a  result,  for  an  electron  ring  that  is  situat¬ 
ed  along  the  minor  axis  of  the  torus,  the  decrease  in  its  equi¬ 
librium  radius  associated  with  the  reduction  of  y0  is  greater 
than  the  corresponding  increase  of  the  equilibrium  radius 
associated  with  the  enhanced  hoop  forces  and  thus  equilibri¬ 
um  can  be  lost.  This  difficulty  can  be  avoided  by  placing  a  set 
of  external  conductors  along  the  minor  cross  section  of  the 
torus  having  a  poloidal  distribution  that  closely  resembles 
the  distribution  of  wall  currents  in  a  perfect  conductor,  i.e., 

Iw  -  —  (/„/27ra)[l  -  (a/2r0\  cos  <p  ], 


where  Iw  is  the  wall  current  per  unit  length  and  lb  the  ring 
current. 

This  compensation  is  satisfactory  even  when  the  beam 
is  displaced  off  center  and  rotates  around  the  equilibrium 
position  with  <oB ,  provided  that  cjbtd>1.  The  reason  is  that 
the  correction  terra  in  the  fields  at  the  center  of  the  ring,14 
when  the  skin  depth  is  much  greater  than  the  thickness  of  the 
conductor,  is 


(  —  2e  ,,Td/o)btd)  sin  aBt, 

and  therefore,  it  can  be  neglected.  Similarly,  when  the  skin 
depth  is  much  smaller  than  the  thickness,  the  correction 
term  is  also  small  and  is  given  by 

(2 [b-a)/a\n  .  2e-',TD  . 

(-5 - — )  +  - sin  wBt. 

In  addition,  we  have  shown  that  the  components  of  the 
magnetic  field  that  are  proportional  to  the  displacement  of 
the  beam  do  not  diffuse  out  of  the  chamber  when  aB  td  >  yj . 
As  a  consequence  n,  rj /a2  does  not  increase  during  diffusion 
and  o)B  does  not  change  polarity.  Therefore,  the  drag  insta¬ 
bility 14  can  be  avoided  by  choosing  the  various  parameters  to 
give  (oB  >  0  at  the  commencement  of  the  diffusion  process. 

During  the  main  acceleration  phase  the  significance  of 
toroidal  effects  is  reduced  because  v/y0-* 0.  When  Be  re¬ 
mains  approximately  constant  (as  shown  in  Fig.  10)  or  in¬ 
creases  with  time,  the  accelerated  beam  moves  closer  to  the 


t 

Injection 


FIG.  10.  This  figure  illustrates  a  possible  time  dependence  of  the  toroidal 
and  betatron  fields.  From  the  stability  point  of  view  other  lime  variations 
are  more  appropriate. 


center  of  the  minor  cross  section  of  the  torus.  This  may  be 
seen  from  the  beam  equations,  which,  for  Boe  =  const,  are 

Ar  +  i-Ar  +  a>2,Ar  =  Q!£-Az+  |4o) 

Y  Y  Y  ymr0 

Az  +  [y/y]Az  +  cb)Az  =  ~{f200/y]Ar.  (41) 

With  the  exception  of  the  second  term  ( y/y]Ar ,  Eqs.  (40)  and 
(41)  are  identical  to  Eqs.  (12)  and  (13).  For  n=\,  Eqs.  (40) 
and  (41)  can  be  readily  solved.  Introducing  a  new  variable 
ip  =  Ar  +  iAz,  these  two  equations  can  be  combined  into  a 
single  equation 

i>+-{y+  inoe)j  +  =  ^  ,  (42) 

y  Y  Ymro 

where  o>\  =  a)  —  The  general  solution  of  Eq.  (42)  is 

if  =  fa  +  fa+fa. 

where 

=  a  1  cxp!  —  *  ^  rf/') 

'  Ym[Woe/Y?  +  4<olYIA' 

_  a 2  exp(  -  /  Jq  a/"1  dt ') 

7'l/W,/?f +  4^]l/4’ 

and 

j,  =j,  f  t  ~  (^oe/r)((SPg)/ymr0\i/>2]  dt’ 

’  Jo  w[fa,fa) 

,  .  p  ( (A>s /Y)( (SPg )/ymr0 ) ]fadt' 

2  Jo  w(p„fa) 

In  the  above  equations  ru1®’  =  (fl00/2y)  ±  [(/20ft/ 
2y)?  +  col  ] 1/2  and  w[ip„p2)  is  the  Wronskian  of  the  two  inde¬ 
pendent  solutions  i/>,  and  fa  of  the  homogeneous  equation. 

Since  the  denominators  of  ip,  and  ip2  increase  when  the 
betatron  field  increases  in  time,  the  center  of  the  beam  moves 
toward  the  center  of  the  minor  cross  section  of  the  torus 
during  the  acceleration. 
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FIG.  1 1 .  System  of  coordinates  used  to  analyze  the  drift  resulting  from  the 
gradient  in  the  toroidal  magnetic  field. 


V.  GRAD  B  DRIFT  IN  THE  MODIFIED  BETATRON 

Until  now  we  have  neglected  in  our  analysis  the  radial 
gradient  in  the  toroidal  magnetic  fields .  In  this  section  we 
consider  carefully  the  effect  of  this  nonlinearity  on  particle 
motion,  assuming  that  it  is  the  dominant  nonlinear  effect. 

The  linearized  equations  of  motion  of  a  particle  located 
at  r  =  r0  -f  Ar  -f  Sr,  z  =  Az  4-  Sz,  where  Ar  and  Az  define 
the  beam  position  with  respect  to  the  center  of  the  torus  (see 
Fig.  11),  are 

r,  +  (1  -  n]n2l0rl-n,nla(sr+ 

'  r0/  Yomro 

=  "M1  ~ 


where 

r,  s=  Ar  +  Sr,  z,  =  Az  +  Sz, 

(43) 

n,0  —  fto,^Yo<  ft  6  O  =  ft  OB  /Vo- 

Choosing  n  —  J  simplifies  the  subsequent  analysis. 
Making  this  choice  and  defining  ip  =  r  +  iz  gives 

+  \ft\a^\  ~  nfiloW  +  (>l/ai)A'l>] 

=  -  ift c  1 1  -  (Re  4-  ft,ol&f>e/Yomro)- 

(44) 

We  proceed  to  solve  (44)  perturbatively,  assuming  Re  ip,  <r0. 

The  zero-order  equation,  neglecting  the  nonlinearity,  is 
easily  solved.  First  an  averaged  is  performed  over  initial  po¬ 
sitions  and  velocities  of  the  particles  to  obtain  a  single  equa¬ 
tion  governing  the  motion  of  the  beam  center  A  ip.  Denoting 
this  average  by  brackets  we  will  have  (Sip)  =  0  by  definition 
and,  as  long  as  no  kinks  develop  in  the  beam,  it  may  be  shown 
that  (Sip)  =0  =  (Sip).  Once  the  average  equation  is  ob¬ 


tained  it  may  be  subtracted  from  (44)  to  obtain  an  equation 
for  Sip.  Carrying  out  this  program,  we  find  for  the  zero-order 
solution: 


Ai/>'°>  =  -= - \°L±>  - +  Ae~  *“*'  +  Be~  '"v, 

ft,  oYomro  m -id  /°K  ]  (45  j 

Sip «»  =  -I-  Ce  ~  -f  De~ 

ft,oY</»ro([  -  ns) 

where  A,  B,  C,  and  D  are  complex  constants  dependent  on 
initial  coordinates  and  where  the  frequencies  &> ,— <u4  are  giv¬ 
en  by 

o>i.z  =#*o±  l-Geo  +4/2^  [J -('»/*>*  ))WS). 

(46) 

^3.4  =  if^SO  ±  [ft  00  +  4/2  Jj(J  —  «,)]  ,/J). 

We  shall  take  subscripts  1  and  3  to  correspond  to  the  -f 
signs. 

In  writing  (45)  we  have  assumed  that  neither 
[i  -  ( d/a 2  )n,  ]  nor  (J  —  n, )  is  zero.  If  either  of  these  quanti¬ 
ties  does  vanish  (corresponding  physically  to  the  vanishing 
of  net  radial  restoring  forces)  the  corresponding  solutions  to 
(44)  grow  secularly,  indicating  a  curvature  or  centrifugal 
drift  in  the  vertical  direction.  Below  we  assume  that  the  radi¬ 
al  restoring  forces  do  not  vanish  for  either  beam  or  particle 
motion.  In  addition  we  shall  make  the  assumption  that 
(SPe)  =  0.  This  is  the  same  as  the  requirement  that  the  equi¬ 
librium  position  of  the  beam  be  at  the  center  of  the  minor 
cross  section  of  the  torus.  This  assumption  has  no  effect 
whatsoever  on  the  basic  physical  results  and  conclusions  but 
does  simplify  the  mathematics  somewhat. 

We  may  return  to  (44)  and  calculate  the  first-order  cor¬ 
rection  to  the  beam  position  A  ipn).  The  equation  to  be  solved 
is 

Aty"  4-  +  -  (rj/a1)/!,]^’" 

=  MjReiOtf*).  (47) 

r0 

Substituting  from  (45)  we  find  the  right-hand  side  of  (47)  is 
<(Re #•)#■>  =  -I  i  i  <ojPjpke-'“f~°> 

J~ 1  A«* I 


+  Ce~^+De 


—  <1 «■>/  -  a} 


X cosK /  -  a* )  -  -j-  V  Vj(p7j),  (48) 

2  /-i 

where  we  have  defined 

A,B,Cfi  ^pf',  y=l, 2,3,4, 
and  where  we  have  assumed  that 
<C>  =  <D>=  0, 

with  the  average  taken  over  the  initial  positions  and  veloc¬ 
ities  of  a  particle. 

From  (47)  and  (48)  we  can  see  that  apart  from  the  oscil¬ 
lating  terms  the  net  effect  of  the  radial  gradient  in  Be  is  to 
cause  an  outward  shift  in  the  equilibrium  position  of  the 
beam: 

,  /in  ,  <Oj  pj) 

Aip,'  =  — =— - - — -  +  oscillating  terms. 

n*>[\-(d/a KJ 

(49) 
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This  result  is  simply  understood  as  the  result  of  a  balance 
between  the  outward,  "diamagnetic”  force  (which  tends  to 
expel  the  beam  frc.  -  the  high  field  region)  and  the  inward 
radial  restoring  force. 

Since  the/>y  values  in  (49)  depend  on  the  details  of  the 
injection  process,  it  is  difficult  to  draw  practical,  quantita¬ 
tive  conclusions  from  this  result.  However  it  is  probably  safe 
•o  conclude  quite  generally  that  any  device  should  be  de¬ 
signed  so  that  ^(rj/a2)/*,  for  all  times,  or 

(50) 

for  the  case  of  a  perfectly  conducting  wall.  The  Yo  should  be 
omitted  in  the  parentheses  on  the  right-hand  side  of  (50)  in 
the  case  of  a  poorly  conducting  wall  (diffusion  time  short 
compared  to  a  beam  oscillation  period).  We  note  that  the 
constraint  in  (50)  is  independent  of  the  strength  of  the  toroi¬ 
dal  magnetic  field. 

Two-dimensional  computer  simulations  bear  out  our 
claim  that  self-consistent  beam  equilibria  exist  in  the  pres¬ 
ence  of  a  gradient  in  Be,  as  long  as  the  net  radial  focusing 
forces  [proportional,  basically,  to  the  denominator  in  (49) 
but  generalized  to  include  the  case  n  ^  $,  and  to  include  toroi¬ 
dal  corrections  to  the  self-fields]  do  not  vanish.  In  Fig.  12  we 
show  a  succession  of  “snapshots”  of  a  beam  cross  section, 
which  remains  in  its  equilibrium  position  for  significant 
times  compared  to  rb/VD,  where  VD  is  the  single-particle 
drift  velocity, 

V D  —  fi0op*/2ro, 

and  where p  is  a  particle  gy radius.  No  drift  is  observed.  (This 
is  not  just  a  visual  observation  but  is  obtained  from  a  plot  of 
average  particle  position  versus  time.) 

Ib  -  10  kA  Eb  •  3  MeV 

|  r0  •  lOOcm  o  •  6.4  cm  Fj,  •  I  cm 

Bz  •  160  G  B  ■  1.4  kG  n  •  0.3 


FIG.  12.  Four  snapshots  of  the  electron-ring  minor  cross  section.  The  val¬ 
ues  of  the  various  parameters  used  in  this  computer  simulation  run  are 
shown  at  the  top  of  the  figure. 


r0  •  100  cm  o  •  6.4  cm  rb  *  I  cm 

B2  •  150  G  B  •  1.4  kG  n  •  0  3 


FIG.  13.  Reducing  the  betatron  field  by  10  G,  but  keeping  the  remaining 
parameters  fixed,  as  in  the  run  of  Fig.  12,  the  electron  ring  drifts  rapidly  in 
the  vertical  direction. 


In  Fig.  1 3,  on  the  other  hand,  we  present  a  case  in  which 
net  radial  focusing  does  nearly  vanish.  Loss  of  confinement 
is  extremely  rapid  under  such  conditions.  The  ring  drifts 
vertically  with  an  average  speed  of  0.25  cm/nsec. 


VI.  THE  EFFECT  OF  EMITTANCE 


Up  to  this  point,  we  have  dealt  with  the  equations  de¬ 
scribing  the  motion  of  the  center  of  an  electron  ring  in  a 
modified-betatron  configuration.  In  this  section,  we  discuss 
the  effect  of  the  finite  emittance  on  the  equilibrium  of  the 
gyrating  electron  ring  using  the  beam  envelope  equation. 
When  the  major  beam  radius  r0  is  large,  n  =  ],  y0  =  const,  v/ 
Y<  1,  the  thermal  energy  spread  Ayo  —  0,  and  the  effect  of 
surrounding  walls  is  neglected,  the  beam  envelope  equation 
in  the  paraxial  approximation  for  Bog  >B0x  becomes, 13  in  the 
Larmor  frame  of  reference, 


*w  + 


2v/ro 

Plrb(s) 


(51) 


where  e  is  the  beam  emittance  (unnormalized),  s  =  0ro  is  the 
length  along  the  minor  axis  of  the  torus,  and  r'b(s)=drb/ds. 
For  a  zero-emittance  beam  the  motion  of  the  particles  is 
laminar  and  the  equilibrium  is  called  either  laminar  flow  or 
Brillouin  flow.  The  equilibrium  radius  in  this  case  is  ob¬ 
tained  from  Eq.  (51)  by  setting 


r;  =  f  =  0, 


and  is 


ffceo  te(2c/f20fl)(2v/j/0)'n,  (52) 

when  y0>\.  For  a  finite  emittance  beam,  the  equilibrium 
radius  can  be  determined  from  Eq.  (51)  by  taking  r£(s)  =  0 
and  is 
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RADIAL  ENVELOPE. 

8S ■  1 4 i 5  G  ;  f  *0 


RING  CROSS  SECTION 

«•  50  m rod 'Cm 


0.0  1.0  2.0  3.0  4.0  3.0 

TIME  (nsec) 

FIG.  1 4.  Root-mean  distance  of  the  electrons  from  the  center  of  the  ring  as  a 
function  of  time,  when  the  emittance  ( is  zero. 


For  small  envelope  oscillations  rb  =  rfceq  +  5,  with 
6  //»_„,  <  1,  and  Eq.  (51)  gives 

6  '  +  2 [ (noe/2ydcf  +  ]S  =  0, 
which  has  a  period 

T  =  - - .  (54) 

2[(nM/2r„c)2  +  fVr^q]l/2 

The  effect  of  emittance  on  the  equilibrium  of  the  ring 
has  been  studied  extensively  using  a  computer-simulation 
code.  Numerical  results  from  the  computer  simulation  are 
given  in  Figs.  14-17.  For  a  beam  with  y0—  7, 1  —  \0  kA, 
r„—  100  cm  ,B0l  =  160  G.and  Bog  =  1415G,  Eq.  (52)  pre¬ 
dicts  that  the  equilibrium  radius  for  the  Brillouin  flow  is  1 
cm.  The  numerical  results  of  Fig.  14  also  give  a  radius  of  1 
cm,  that  for  all  practical  purposes  remains  constant  in  time. 
In  this  run  the  electron  beam  is  injected  into  the  torus  with  a 
rotational  frequency  that  is  half  of  the  local  cyclotron  fre¬ 
quency.  Figure  1 5  shows  the  envelope  in  time  of  a  nonrotat¬ 
ing  beam.  The  various  parameters  in  this  run  have  the  same 
values  as  those  in  Fig.  14,  except  now  the  emittance  is  non 
zero  in  the  Larmor  frame.  For  the  equivalent  emittance  of 
€  =  50  mrad  cm,  Eq.  (53)  gives  /•4cq  =  1.21  cm  and  the  nu- 

o  RADIAL  ENVELOPE 

B#*I4I5G;  «  *50  mrod-cm 


0.0  1.0  Z.O  3.0  4.0  9.0 

TIME  (nsec) 

FIG.  I }.  Root-mean  distance  of  the  electron  from  the  center  of  the  ring  as  a 
function  of  time.  In  this  run  t  =  50  mrad  cm. 
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FIG.  16.  Four  snapshots  of  the  electron-ring  minor  cross  section  when 
€  =  50  mrad  cm. 

merical  results  give  an  equilibrium  radius  approximately 
equal  to  1.24  cm.  Similarly,  Eq.  (54)  gives  a  period  of  2.13 
nsec  and  the  numerical  results  give  2.28  nsec.  Snapshots  of 
the  beam  cross  section  from  the  same  run  are  given  in  Fig. 
14.  The  oscillations  observed  in  the  run  of  Fig.  15  can  be 
avoided  by  “matching”  the  beam,  i.e.,  by  raising  the  Boe 
magnetic  field  to  1830  G.  This  value  of  the  magnetic  field 
gives  an  rbcq  ~  1  cm,  which  is  the  radius  of  the  beam  for 
Brillouin  flow.  Numerical  results  from  this  run  are  shown  in 
Fig.  17. 

The  electron  beams  discussed  so  far  in  this  section  were 
monoenergetic  with  finite  emittance.  Such  beams  have  an 
axial  velocity  spread  equivalent  to  that  of  a  cold  beam  with 
energy  spread  Ay  that  is  given  by 15 

Wr/r),  =  M: rPe/rb)7. 

Hot  beams  have  an  additional  energy  spread.  This  thermal 
energy  spread  in  the  direction  of  beam  propagation  has  an 
important  effect  on  the  dynamics  of  electrons  as  may  be  seen 
as  follows.  The  equations  describing  the  motion  of  individual 
electrons  in  cylindrical  geometry  are  identical  with  those 

RADIAL  ENVELOPE 

O 

9-i  B*’I830G;  «’50mrod'tm 


0.0  1.0  Z.O  3.0  4.0  3.0 

TIME  (nsec) 

FIG.  17.  Root-mean  distance  of  the  electrons  from  the  center  of  the  ring. 
The  envelope  oscillations  are  drastically  reduced  when  the  toroidal  magnet¬ 
ic  field  is  increased  from  1415  G  to  1830  G. 
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describing  the  motion  of  the  center  of  the  beam,  provided 
that  5*  and  <0]  are  replaced  by 

s;  =  (/2oz/ro)J(i 
=  (^o ,/rJri*  -  «.). 

and 

(SPe)  =  c2  Ay 
Wo  y0 

The  above  set  of  equations  are  based  on  the  assumption  that 
toroidal  effects  can  be  neglected.  In  addition,  it  should  be 
emphasized  thatdyis  the  thermal  energy  spread  and  not  the 
energy  mismatch  6y0  discussed  in  Sec.  III. 

For  n  =  $,  the  solution  of  the  individual-particle  equa¬ 
tions  are 

r  —  r0  —  8r0  cos  cogl  —  6z0  sin  eoet  +  r0 

Yo 

(1  —  COS  C0g  t  ) 

*  ~T^r~ '  <55' 

_  _  .  4ysin<u»/ 

oz  =  6z0  cos  aBt  -f  <5r0  sin  coBt  —  r0 — - ,  (56) 

Yo  \~n, 

where 

*>»  =  (^ot/Yo)(B,/B0)(\  -  n,). 

According  to  Eq.  (55),  when  n,  <1,  i.e.,  for  low-current 
beams,  thermal  efTects  substantially  increase  the  minor  radi¬ 
us  of  the  beam.  In  such  beams  the  minor  radius  varies  as 
Irtfiy/Yo- 

In  contrast,  when  n,>l,  i.e.,  for  high-current  beams, 
thermal  effects  do  not  change  significantly  the  minor  radius 
of  the  beam,  which  varies  as  \rjns\hyjy^. 


LOW  CURRENT  (100  A) 
At/t  =  10% 


TIME  -  0.00  me 


HIGH  CURRENT  (10  kA| 


tu/i  ■=  10% 
TIME  ■  0.00  n*ec 


TIME  -  10.00  ns#c 


FIG.  1 8.  Snapshots  of  the  electron-ring  minor  cross  section  for  low  and  high 
current.  The  values  of  the  various  parameters  for  this  run  are  listed  in  Table 
II.  The  energy  spread  in  both  cases  is  10%. 


TABLE  111  Electron-ring  parameters  for  the  two  computer  runs  shown  in 
Fig.  18. 


Parameters 

Low  current 

High  current 

Beam  current  (k A.) 

0.1 

10 

Beam  energy  (MeV) 

3 

3 

Energy  spread  (%) 

10 

10 

Major  radius  (cm) 

100 

100 

Beam  minor  radius  (cm) 

1 

1 

Initial  betatron  field  (G) 

116 

146 

Toroidal  field  (kG) 

1.4 

1.4 

External-field  index 

0.5 

0.35 

Self-field  index  (n,) 

0.37 

23.4 

Torus  minor  radius  (cm) 

3.2 

3.2 

The  effect  of  the  axial  energy  spread  on  the  minor  cross 
section  of  the  beam  in  a  modified-betatron  geometry  has 
been  studied  numerically.  Results  from  both  the  high-  and 
low-current  beams  are  given  in  Fig.  1 8.  In  these  two  runs  (he 
various  parameters  have  the  values  listed  in  Table  III.  In  the 
high-current-beam  case  the  minor  radius  expands  by  ap¬ 
proximately  a  factor  of  two.  However,  in  the  low-current- 
bfun  run  the  beam  expands  significantly  and  strikes  the 
wall.  Therefore,  a  substantial  energy  spread  can  be  tolerated 
in  high-current  beams  without  a  catastrophic  expansion  of 
the  minor  radius  of  the  beam.  Such  energy  spread  may  be 
required  to  stabilize  the  various  disruptive  instabilities.16"18 

VII.  SUMMARY 

In  this  paper  we  review  the  dynamics  of  ultrahigh-cur- 
rent  electron  rings  in  the  modified-betatron  configuration. 
Our  discussion  mainly  addresses  the  evolution  of  the  elec¬ 
tron  ring  after  injection.  The  formation  of  the  ring  during 
injection  has  been  analyzed  and  reported  previously. 19 

Our  work  includes  both  analytical  and  numerical  re¬ 
sults  for  "cold”  and  “hot”  rings.  The  conclusion  of  these 
studies  is  that  equilibrium  states  of  ultrahigh-current  rings  in 
a  modified  betatron  exist  over  a  wide  range  of  parameters. 
These  equilibra  are  realistic  and  accessible  with  state-of-the- 
art  injectors. 

The  results  presented  in  this  paper  are  based  on  several 
simplified  assumptions.  Among  them,  we  have  assumed  that 
the  various  fields  are  free  of  errors,  the  conducting  wall  that 
surrounds  the  electron  ring  was  assumed  to  be  smooth  (i.e., 
without  ports  and  gaps),  and  the  external-field  index  was 
assumed  to  be  constant  in  time  and  space.  The  consequences 
of  these  assumptions  are  presently  under  investigation. 
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APPENDIX  A:  DIFFUSION  OF  THE  SELF-MAGNETIC 
FIELD  OF  THE  BEAM  THROUGH  A  CONDUCTING  LINER 

The  purpose  of  this  appendix  is  to  briefly  outline  (he 
calculation  of  the  diffusion  of  the  self-magnetic  field  of  the 
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beam  through  a  conducting  liner.  To  simplify  the  analysis,  it 
is  assumed  that  the  electron  beam  is  located  along  the  axis  of 
a  straight,  circular  cylinder  of  inner  radius  a,  outer  radius  b, 
and  thickness  S  =  b  —  a,  as  shown  in  Fig.  9.  In  addition,  it  is 
assumed  that  the  current  of  the  electron  beam  is  a  step  func¬ 
tion  that  is  turned  on  at  t  =  0.  Since  the  problem  of  interest  is 
that  of  an  electron  ring  inside  a  torus,  it  is  further  assumed 
that  the  axial  inductive  electric  field  goes  to  zero  at  r  =  r0, 
where  r0  is  the  major  radius  of  the  torus. 

Neglecting  the  displacement  current  {quasistatic  ap¬ 
proximation),  the  fields  inside  the  cylindrical  conducting 
shell  (o<r<6 )  are  given  in  mks  units  by 


VXH  =  J, 

(Al) 

VB  =  0. 

(A2) 

VXE=  - 

dt  ’ 

(A3) 

J  =  crE, 

(A4) 

where  J  is  the  current  density,  H  is  the  magnetic  field,  B  is 
the  magnetic  induction,  and  a  is  the  conductivity  of  the  con¬ 
ductor. 

Assuming  that  all  quantities  are  independent  of  z  and  <f>, 
Eqs.  (A  1 )  to  (A4)  can  be  combined  to  give  a  diffusion  equa¬ 
tion  for  the  azimuthal  component  of  the  magnetic  field 
d  {  1  d  \  dlh 

!(t .  'A5> 

For  r  <a  and  r>  b,  the  conductivity  is  equal  to  zero  and  Eq. 


(A  5)  becomes 

(A6) 

with  the  only  acceptable  solutions 

)  =  (/ /2irr)0(t ),  r*  <r<a, 

(A7) 

Ht[r,t)  =  H4[b,t  )b/r,  r>b. 

(A  8) 

where  O  [l )  is  the  step  function  and  rb  is  the  beam  radius. 

To  complete  the  specification  of  the  problem  we  have  to 
introduce  boundary  conditions.  From  Eqs.  (A3)  and  ( A8),  we 
get 

BE,  _  ..  dH,[b,t)  b 
dr  dt  r 

which  after  integration  yields 

E,  ( r,t )  =  p  dH*£i'±  b  ln(-^  ,  (A9) 

where  we  have  assumed  that  E,(r,t )  is  zero  at  r  =  r0. 

Combining  Eqs.  (Al),  (A4),  and  (A9)  at  r  =  b,  we  obtain 
the  first  boundary  condition,  namely 


(t  ,M01 


The  second  boundary  condition  is  furnished  from  the  con¬ 
tinuity  of  the  magnetic  field  at  r  =  a,  i.e., 

{I/2ira)0[i)  =  N^a,t).  (All) 

Since  the  magnetic  field  is  zero  at  f  =  0,  Eqs.  (A5),  (A9),  and 
(A  10)  take  the  following  forms  in  their  transformed  state 

_  A  * 

d'H* 


dr2 


L  +  T^T  ~ {?  +  afip)fl*  =  °*  (A12) 


( j UrH4 1),  ^  =  appH.  (b,t  )b 1  ln(A)  , 


and 


(/  /2-rtap)  =  [H4{r,t  )],_„, 
where,  the  Laplace  transform  of  H ^  is  defined  by 


(A  13) 

(A  14) 


tf*{r,p)=  f 
Jo 


e-'”H+(r,t)dt.  Re(  p)  >  r. 


Equation  (A12)  is  the  modified  Bessel  equation  and  its  solu¬ 
tion  is 


H+(r,p)=AIx(Ar)  +  BKl[Ar), 


where  A  and  B  are  constants  and  A  2  =  app. 

The  two  constant  coefficients  .4  and  B,  determined  from 
Eqs.  (A  13)  and  (A  14),  are 

A  =  (I  /2irap)[K0{Ab )  +  Ab  ln(6  /r0)Kx(Ab  ))/A, 

B  =  (I/2nap)[I0(Ab )  -  Ab  In (b/rJI^Ab  )]/A, 
where 

A  =  K0[Ab  ]I,[Aa)  +  I0(Ab  )K,(Aa )  +  Ab  \n[b  /r0) 

X[Kt{Ab)r,(Aa)  ~  XifibVMb)).  (A  15) 

The  magnetic  field  in  the  region  a  <r  <b  as  a  function  of  space 
and  time  can  be  obtained  by  inverting  the  Laplace  transform, 
i.e., 


J_  r4'"  dpt" 

2777  Jr  _  ,  „  2irap 

X  j  ( [x-0(/lb )  +  Ab  ln^*  ,(Ab  )j/,(/lr) 

+  [l0\Ab )  -  Ab  In )]AT,(/lr)J,  (A16) 


where  the  path  of  integration  is  a  vertical  line  in  the  complex 
P  plane  to  the  right  of  all  singularities  of  the  integrand. 

Equation  (A  16)  has  a  simple  pole  atp  =  Oof  residue  a/r 
and  an  infinite  set  of  simple  poles  at  A  =0. 

Contour  integration  ofEq.  (A  16)  gives 


^(r,t)=y-l-+  £exp(— ^-) 

2ira  I  r  , «  l  \  a/i0  / 

xjjyor.ft  */,(«,/■)-  Y^a.r^a.b)  +  a,b  ln(^)[  Yx[a,b  )/,{a,r)  -  Y^r^a.b)]  ]  J  ]~‘  (A17) 

where  p  =  /r0,  A  =  iff,  Fand  J  are  the  Bessel  functions,  a,  are  the  roots  of  equation 
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Y0{a,b\Ji(a,a)~  Yilfl.aVJfit.b  )  +  a,6  ln(6/r0)[  Y^a.b  y,(a,a)  -  ^(a.ay.la.h )]  =0.  (A18) 

and 

a\~)a  o  =  a,a[Y0(a,b)J0(a,a)  -  Y0(a,a)J0(a,b }]  -  a,b  [  Y^a.b  \Jt(a,a)  -  Y^a^^a.b)] 

+  a,b\n(b/r0)\ba,[Y0{a,b)Jt{a1a)  -  Y^a^^a.b)]  +  a,a[70(aIa)}'j(a,h)--  Y0[a,a\J^a,b)\  J. 
Fora,a>  1,  Eq.  (A  17)  is  reduced  to 

X  |ar,<5  sin  (a,6)  —  a]b ln^— j  cos  a,S  j  |, 


(A19) 


where  the  a,  are  now  the  roots  of 

cos(ar,<5)  +asb  ln(6  /r0)  sin(a,<5 )  =  0.  (A20) 

Equations  (A  19)  and  (A20)  are  further  simplified  when 
a,<5<l.  In  this  case 

"'|r',,=  2£'l7_(7)’1“s“-(>-'’1 

+  a,b  ln^-j  sin  a,(b  —  r)j  exp(  — — jj, 

which  at  r  =  b,  becomes 

(b,t )  =  (/  /2vb )  [  1  —  exp(  —  a)t  / ofi0 )  ] , 

and 

a]=  [8b\n[ro/b)}-'.  (A21) 

Therefore,  to  lowest  order,  the  magnetic  field  in  the  region 
r>b  is 

H* ^ ) =  2wr  [ 1  ~  CXP(  r>b  *mkS*’ 

cp0  -  (a.22) 

or 

H4  ( r,t )  =  (21  /rc)  [  1  —  exp(  —  a2,tc7/4va)  ]  (cgs). 

(A23) 

Under  the  same  approximations,  the  electric  field  of  the 
inner  edge  (r  —  a)  of  the  conducting  shell  is 

EAa,t)=  -  ~~a\  ln(y)  exp('~g  ~  )  <mks)> 

(A24) 

Since  the  electric  field  is  uniform  in  the  region  r<,a,  Eq.  ( A24) 
also  gives  the  electric  field  that  acts  on  the  beam.  Substitut¬ 
ing  Eq.  (A24)  into  the  energy-rate  equation 

me2  —  =  —  evE, 
dt 

we  obtain  for  highly  relativistic  beams 

Ay/y~  -  2(v/y)  ln(ro /b ).  (A25) 

For  tf/b  —  10 and  v/y  =  0.1,  Eq.  (A25)  gives  Ay/y  =  0.46, 
i.e.,  a  substantia]  reduction  in  the  energy  of  the  beam.  How¬ 
ever,  for  highly  relativistic  beams  v,  ~c  and  thus  the  current 
of  the  beam  remains  approximately  constant. 
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APPENDIX  B:  DESCRIPTION  OF  THE  PARTICLE-IN¬ 
CELL  COMPUTER  CODE 

The  Naval  Research  Laboratory  (NRL)  modified-beta- 
tron  accelerator  is  designed  for  a  maximum  v/y0~ 0.1,  i.e., 
the  current  is  high  enough  so  that  the  self-fields  of  the  ring 
can  exceed  the  externally  applied  fields.  A  realistic  theoreti¬ 
cal  description  must  therefore  self-consistently  include  the 
beam’s  self-fields  as  well  as  the  effect  of  surrounding  walls. 
Since  this  is  difficult  analytically,  particularly  if  the  ring  is 
displaced  substantially  from  the  center  of  the  minor  cross 
section  of  the  torus,  numerical  simulations  are  useful  both  in 
gaining  insight  into  the  important  physical  processes  as  well 
as  to  provide  a  method  to  check  the  applicability  of  specific 
assumptions  in  an  analytic  model. 

The  dynamics  of  the  accelerated  electron  ring  are  deter¬ 
mined  by  forces  that  vary  on  a  number  of  different  time 
scales  which  range  from  the  electron  cyclotron  period,  i.e.,  a 
few  nanoseconds,  to  the  beam  acceleration  time,  which  is  of 
order  of  a  millisecond.  The  code  described  here  is  tailored  to 
simulate  efficiently  the  various  phenomena  on  the  interme¬ 
diate  time  scale.  This  time  scale  is  characteristic  of  the  drift 
(bounce)  motion  of  the  ring  after  equilibrium  has  been  estab¬ 
lished,  rather  than  the  rapid  evolution  occurring  at  injection. 
Simulation  of  a  single  turn  around  the  major  axis  that  lasts 
about  20  nsec  using  4000  particles  on  a  64  X  64  grid  typically 
takes  about  one  minute  on  the  NRLTexas  Instruments  ASC 
(Advanced  Scientific  Computer). 

The  simulation  co  J->  is  r  —  z,  spatially  two-dimension¬ 
al,  i.e.,  d/dO  —  0,  but  with  three  velocity  components.  Al¬ 
though  Be  is  used  in  calculating  the  particle  trajectories,  it  is 
not  solved  self-consistently,  i.e.,  it  is  assumed  to  be  generated 
from  external  coils  only.  This  assumption  is  valid  to  first 
order  in  v/y.  The  radiative  term  (displacement  current)  is 
also  ignored,  i.e.,  the  code  uses  the  Darwin  model  for  Max¬ 
well’s  equations. 

The  electrostatic  potential  is  computed  from  Poisson’s 
equation 

V20=p/eo,  (Bl) 

and  the  magnetic  vector  potential  from 

V2Ae-A0t?  =  -nJe,  (B2) 

with  the  boundary  condition  <P  —  A0  =  0  at  the  conducting 
wall. 

Equations  (B 1 )  and  (B2)  are  solved  by  Fourier  decompo¬ 
sition  in  thez  direction  and  then  by  Gaussian  elimination  of 
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the  resultant  tridiagonal  matrix  of  equations  obtained  from  a 
three-point  differencing  scheme  for  Y2.  The  inverse  Fourier 
transform  yields  A0  and  <P  on  the  grid.  Note  that  the  8- 
particlc  velocities  are  advanced  using  the  conservation  of 
canonical  momentum  in  the  6  direction,  the  equation  for  A  0 
is  therefore  not  properly  time-centered  since  the  velocities 
from  the  previous  time  step  are  used  to  calculate  the  currents 
from  the  canonical  momenta.  This  method  was  chosen  pri¬ 
marily  for  its  speed  and  simplicity  but  care  must  be  taken  in 
applying  the  code  when  the  inductive  acceleration  of  parti¬ 
cles  in  the  6  direction  is  significant. 

If  boundary  conditions  other  than  Ae  or  <t>  equaling 
zero  on  a  rectangular  grid  are  desired,  it  is  possible  to  obtain 
relatively  arbitrary  boundary  conditions  using  the  capaci- 
tive-matrix(Buneman)  technique.20  In  this  method  a  matrix 
is  generated  which  is  the  Green’s  function  for  discrete  “wall” 
points  within  the  system.  Then  at  each  time  step  the  field 
solver  described  above  is  used.  The  potential  at  the  discrete 
(wall)  points  is  obtained.  By  multiplying  this  vector  potential 
by  the  invense  matrix  obtained  previously  a  set  of  (wall)  cur¬ 
rents  or  changes  is  generated,  which  is  used  to  specify  the 
wall  potential.  These  (wall)  sources  are  then  added  to  the 
original  beam  source  and  the  field  solver  is  used  again.  The 
result  by  the  principle  of  superposition  is  correct  inside  the 
system  and  has  the  correct  boundary  condition  on  the 
“wall." 

The  motion  of  the  electrons  is  governed  by  the  Lorentz 
force 

=  J-IE+VXB), 
at  m0 

where  E  and  B  arc  the  total  electric  and  magnetic  field,  re¬ 
spectively. 

In  component  form  the  equations  used  to  update  the 
velocities  and  positions  at  each  time  step  are 

yVg  =  Pg/mj  -  (q/m0]Ae,  (B3) 


U„xn,  (  U,xH0 


=  3-  (Er  +  VgxB ,  -  V'XBoi  +  — ,  (B4) 

dt  m0  r 

iMA  =  J-(Et  -  VexBr  +  VrxBe ),  (B5) 

dt  m0 

(B6) 

dt 

—  =  VI.  (B7) 

dt 

Equation  (B3)  is  used  at  each  timestep  to  compute  yV0. 
Equations  (B4)  and  (B5)  are  coupled.  To  advance  these  veloc¬ 
ities  a  leap-frog  scheme  is  employed.  At  time  /  all  fields, 
positions,  and  Ve  are  known  exactly.  Vr  and  V,  however  are 
known  at  t  ■+  1  and  the  velocities  will  be  advanced  to  t  + 
Before  the  equations  are  differenced,  it  is  convenient  to  re¬ 
write  the  equations  in  terms  of  the  relativistic  momenta.  Let¬ 
ting 

u  =  rV, 

so  Eqs.  (B3),  (B4),  and  (B5)  can  be  written  as 

Ug  =  Pg  /m,f  -  ( q/tn0]A  g,  (B3'j 

=  JL-  {rE,  +  UexB ,  -  U,xB0)  +  ,  (B4') 

dt  ym0  yr 


This  formulation  is  then  differenced  by  substituting 


U' +  1/2  +  U'  - 
2 


for  LJ. 


After  making  these  substitutions  in  Eq.  (B4')  and  (B5')  it  is 
straightforward,  if  somewhat  tedious,  to  solve  the  two  cou¬ 
pled  equations  for  t/'+,/2  and  U'r+I/ 2.  Since  t/,  and  V ,  are 
calculated  at  the  half-timesteps  while  U0  is  known  at  the  full 
timesteps,  y  is  not  known  at  time  t.  This  difficulty  is  over¬ 
come  using  iteration. 

After  advancing  the  velocities,  the  particle  positions  are 
advanced  using  a  simple  centered  difference 
i*"  —  1*  +AtUr/y, 
and 

z'+l  =  z'  +AtU'/y. 

Once  all  the  velocities  and  positions  are  found  the  new  cur¬ 
rent  and  charge  densities  are  used  to  update  the  fields  and  a 
new  time  step  begins. 
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The  integer  resonances  affecting  beam  motion  in  the  presence  of  externa)  field  imperfections  in  the  modified  betatron 
are  studied.  An  upper  bound  is  obtained  on  the  magnitude  of  field  error  that  may  be  tolerated.  A  numerical  example 
shows  that  for  practical  parameters  the  resulting  bound  is  very  restrictive.  The  effect  of  longitudinal  temperature 
and  other  possible  stabilizing  effects  are  discussed. 


I.  INTRODUCTION 

In  a  conventional  betatron,  low-order  resonances 
between  particle  motion  and  field  imperfections 
can  be  avoided  by  restricting  the  beam  current 
so  that  the  tune  shift'  remains  sufficiently  small. 
In  a  modified  betatron,  the  addition  of  a  strong 
toroidal  magnetic  field  may  allow  large  currents 
to  be  accelerated,2-3  but  resonances  become 
much  more  difficult  to  avoid,  especially  if  one 
contemplates  removing  the  toroidal  field  before 
the  beam  is  extracted.  This  paper  examines  the 
problem  of  integer  resonances  in  the  modified 
betatron4  and  obtains  a  condition  bounding  the 
rate  of  change  of  the  fields;  when  the  condition 
is  satisfied  the  resonances  are  passed  through 
with  sufficient  speed  so  that  the  beam  is  not  sig¬ 
nificantly  disturbed.  We  consider  here  only  er¬ 
rors  in  the  fields  themselves,  not  in  field  gra¬ 
dients,  so  we  discuss  only  integer,  not  half- 
integer  resonances. 

In  what  follows,  we  will  first  consider  a  “cold” 
beam,  that  is,  one  in  which  there  is  no  spread  in 
longitudinal  energy  or,  therefore,  in  circulation 
frequency  about  the  machine.  The  effect  of  or¬ 
bital  resonances  on  such  a  cold  beam  will  be  seen 
to  place  rather  severe  limits  on  the  magnitude  of 
the  tolerable  field  imperfections.  When  the  ef¬ 
fects  of  temperature  are  taken  into  account,  how¬ 
ever,  a  numerical  example  below  will  illustrate 
a  reduction  of  the  effect  of  the  resonance  on  the 

*  Supported  by  the  Office  of  Naval  Research. 


motion  of  the  beam  center  of  mass.  An  expla¬ 
nation  of  this  temperature  effect  will  be  given. 

II.  ORBITAL  RESONANCES  FOR  A  COLD 
BEAM 

We  consider  a  beam  of  circular  cross  section  and 
uniform  density  and  current  profiles,  as  shown 
in  Fig.  1.  The  torus  has  a  major  radius  r0  and 
minor  radius  a;  the  chamber  is  assumed  to  be 
perfectly  conducting  as  far  as  the  rapidly  varying 
part  of  the  self-fields  is  concerned.  The  beam 
radius  is  rb  with  center  located  at  r  =  r0  +  A r, 
z  =  Az,  as  shown  in  the  figure.  If  we  define  the 
displacement  of  a  particle  from  the  design  orbit 
r  =  r0,  z  =  0  as  rx  =  A r  +  8r,  z,  =  Az  +  8z 
then  the  equations  of  motion  for  ri  and  z,  are,  to 
first  order  in  the  displacement  from  the  design 
orbit 

rj  +  -^ri  +  0*o0  -  n  -  ns)rt 


2ntyoC 


Zi  +  HeoZi 


- I  Er  +  Po Bl 

m  y<>  l 

+  J'  dt'E'W)  j 
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FIGURE  1  Minor  cross  section  of  modified  betatron  show¬ 
ing  beam  center  and  particle  coordinates.  The  major  radius 
of  the  device  is  r0. 


Zi  +  —  Zt  +  n|n(n  -  n,)zi 
7o 


where  70  is  the  particle  energy  on  the  "ideal” 
design  orbit  (in  units  of  me2),  f lz0  =  eB^lmyoC, 
fteo  =  eBaohn'1oC,  Bt 0  is  the  vertical  magnetic 
field  at  the  design  orbit,  Be 0  is  the  toroidal  mag¬ 
netic  field,  n  is  the  betatron  field  index  (assumed 
constant  here),  n,  =  wfc2/(2-y02f^)  is  the  self-field 
index,  w*  =  ( 4Ttn0e7/myoy' 2  is  the  beam  plasma 
frequency  and  e,  m  are  the  magnitude  of  electron 
charge  and  rest  mass,  and  where  field  compo¬ 
nents  with  a  “wiggle”  on  the  right  hand  sides  of 
(la,  b)  denote  the  value  of  the  field  imperfection 
which  in  general  will  depend  on  the  value  of  0, 
the  azimuthal  position  of  the  particle.  In  deriving 
(la,  b)  we  have  allowed  all  fields  to  depend  on 
time;  we  have  therefore  included  the  inductive 
poloidal  electric  field  (Bw  terms).  Also  included 
are  the  effect  of  wall  image  charges  and  currents 


( rb2la 2  terms),  present  when  the  beam  is  dis¬ 
placed  from  the  center  of  the  chamber. 

We  desire  to  have  equations  which  describe 
only  the  motion  of  the  beam  center,  Ar(0,  t), 
Az(0,  t).  To  this  end  we  define  a  distribution  func¬ 
tion  /  as 


f(r,  0,  z,  v„  ve,  vt,  /)  =  £  g(r(0),  v<0>) 

(40)^0) 

X  8(r  -  f)  -6  ~  —  8 (z  -  £)8(3)(v  -  *),  (2) 


where  r(0)  and  v(0)  are  particle  initial  conditions, 
f,  0,  i,  and  $  are  the  solutions  for  the  particle 
trajectories  as  functions  of  initial  position,  ve¬ 
locity,  and  time,  and  where  g(r<0>,  v(0>)  is  a  weight¬ 
ing  function.  We  then  have  that 


A r(0,  l ) 


s  r  dr  dz  d\{r  -  r0)f 
f  r  dr  dz  dvf 


S  g(r(0),  ~  ro)8(0  -  0) 

2  £(r(0).  vlO))8(0  -  0) 


(3) 


A  similar  expression  holds  for  Az(0,  t).  It  may 
be  similarly  shown  that,  for  a  cold  beam, 

<h>  -  (i, +  n"s)'4r  '4) 

<i,'>'(s  +  n»s) ir-  <5) 

where,  of  course,  analogous  expressions  hold  for 
(zi),  (51)  and  (tj).  In  Eqs.  (4)  and  (5)  we  have 
assumed  that  all  particles  circulate  the  machine 

with  0  =  f I#.  This  assumption  will  be  relaxed 
in  the  next  section  where  the  effects  of  finite 
longitudinal  temperature  are  considered. 

Using  this  averaging  procedure  on  Eqs.  (Ia,  b), 
one  obtains  equations  for  the  beam-center  mo¬ 
tion.  Though  these  may  be  solved  in  general,  the 
special  choice  n  =  i  (which  is  consistent  with 
our  assumption  of  a  circular  beam)  simplifies  the 
analysis.  With  n  —  {  and  defining 

A«Ji  =  Ar  +  i'Az  =  2  Ail ucm,  (6) 
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the  equation  for  At},,  is 


resonance).  Evaluating  Eq.  (8)  then  gives 


Ai|>,  +  f—  +  ,fieo  +  2//nJ 
dl  |_f°  J  dt 

*  [n;"G '  ^  +  "n“w 


+  i  - - —  IClooClio  +  i/ftjo  I  Ai|», 

2my0c  J 


where  F,  is  the  /-th  Fourier  component  of 


-  —  \'Er 
my  o 


+  PoFc  +  i(Ez  —  PoBr) 


ft*o  J/'£e(/')^j  . 


Equation  (7)  may  be  solved,  assuming  the  func¬ 
tions  multiplying  derivatives  of  A«|i,  are  slowly 
varying  over  the  period  of  a  betatron  oscillation: 

Ai{//  ~  (-yo<D0)~ in  J  dt'e~iS‘drli^eo  +  lilto1  ^ 

where 

wo(0  33  Jo  n^j  +  4^®°j  •  (9) 

For  long  times  (many  betatron  periods),  the  in¬ 
tegral  in  (8)  may  be  evaluated  by  the  method  of 
stationary  phase.  The  points  of  stationary  phase 
(resonance  points)  occur  when 

ft,~  =  —  Jfteo  —  /ft*o  ±  wo  =  0  (10) 

for  a  given  /.  This  is  just  the  condition  that  the 
betatron  frequency  be  /  times  the  fundamental 
cyclotron  frequency,  (1*0.  Condition  (10)  may 
also  be  written 

Be o  =  “  7  n>  ~  2)  ^zo  ■  (11) 

For  positive  Bw  and  Bt0,  Eq.  (1 1)  may  be  satisfied 
only  by  negative  /  and  for  such  /,  Eq.  (10)  may 
be  satisfied  only  for  the  lower  sign  (fast  mode 


/  \  1/2  r 
'■(!)  [ 


yo(t-) 

yo(  t )  Wo(  1  )o>o(  1 -)_ 


y  o, -dr- itt/4 

|fir(r-)|,/2 

where  /_  is  the  time  at  which  ft,-  =  Oand  where 
the  +  or  —  sign  is  used  in  the  exponent  according 
as  ft,-(/_)  >  Oor  <  0  respectively.  If  we  neglect 
the  possibility  of  cancellation  due  to  different 
phases  as  we  pass  through  different  resonances 
and  if  we  interpret  F,  generically  as  [  -  (e/my0)bfi) 
where  8/,  is  the  /-th  Fourier  component  of  any 
field  error,  we  may  obtain  a  lower  bound 
on  |  ft,-  |  by  requiring 


|  AiJj,|  «  a, 


which  gives 


Iftrl-fMA-T.  (i4, 

2  L/nyoWoaJ 

which  is  our  basic  result.  For  y0  large  enough 
that  we  may  neglect  ft*0  compared  with  fteo  and 
rb2n,la2  compared  with  1/2,  this  constraint  may 
be  rewritten,  using  the  relations 

=  ~FJnh(leo  (,5) 


l2  +  1/2 


I  ft  00  I  ^ 


l2  +  1/2 


[-£T 


As  an  example,  we  consider  the  problem  of 
passing  through  the  /  =  -  1  resonance.  We  con¬ 
sider  a  hypothetical  experiment  (r0  =  1  m,  a  = 
10  cm,  rb  =  1  cm)  in  which  y0  is  increased  linearly 
in  time  from  an  initial  (injection)  value  of  7  to  a 
final  value  (rr,naI  =  1  millisecond)  of  100,  while 
simultaneously  Be „  is  decreased  from  1.5  kG  to 
0.  The  /  =  -  1  resonance  will  occur  at  /  =  627 
p.scc,  at  which  time  Bl0  =  1 120  G,  Bb0  =  560  G 
and  y0  -  65.3.  At  resonance,  ft0o  =  -6.2  x 
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10"  sec-2.  Substituting  in  the  expression  (17)  we 
obtain  an  upper  bound  on  the  allowable  field 
error 


— 1.3  x  I0-4. 

Bzo 

a  rather  severe  requirement. 

We  conclude  that,  at  least  for  the  case  of  a  cold 
beam,  it  may  not  be  desirable  to  remove  the  to¬ 
roidal  field  and  pass  through  these  resonances. 
Perhaps  the  toroidal  field  may  be  reduced  some¬ 
what  from  its  initial  value,  assuming  the  high-/ 
resonances  are  not  too  important  and  can  be 
passed  through  easily.  It  may  then  be  possible, 
by  the  use  of  an  intentionally  introduced  field 
perturbation,  to  use  a  low-/  resonance  in  a  con¬ 
trolled  way  to  extract  the  beam  before  Z?B  is  com¬ 
pletely  removed. 

It  should  be  noted  that  it  is  possible,  at  least 
in  principle,  to  avoid  the  integer  resonances  al¬ 
together  by  raising  both  Bi0  and  Be0  proportion¬ 
ately  and  in  such  a  way  that  condition  (11)  is 
never  satisfied  for  any  /.  At  the  end  of  such  an 
acceleration  cycle,  however,  one  will  have  a  very 
large  toroidal  magnetic  field  in  the  device,  pos¬ 
sibly  complicating  the  extraction  process. 

The  above  results  apply  to  a  beam  all  of  whose 
particles  are  traveling  to  lowest  order  at  the  same 
azimuthal  angular  velocity.  All  particles  are  then 
in  resonance  at  precisely  the  same  moment  and 
receive  the  same  periodic  perturbations  to  their 
orbits.  In  the  next  section  we  relax  this  assump¬ 
tion  and  examine  the  behavior  of  a  beam,  the 
particles  of  which  possess  a  spread  in  energy. 


III.  EFFECT  OF  FINITE  BEAM 

TEMPERATURE  ON  RESONANCES 

To  calculate  the  effect  of  beam  temperature  on 
beam  behavior  near  a  resonance,  we  consider  an 
ensemble  of  beams,  each  cold  and  each  consist¬ 
ing  of  particles  traveling  Avith  a  zero-order  an¬ 
gular  frequency  0  0  given  by 

0o  =  n,0  -  kP,  (18) 

where  P  is  the  canonical  angular  momentum  of 
a  particle,  which  is  related  to  the  difference  in 
energy  between  the  particle  under  consideration 
and  the  (reference)  particle  maintained  at  the  de¬ 


sign  orbit  r  =  r0l  z  =  0  by 
=  Ay  me2 


(19) 


and  where,  in  Eq.  (18), 

(20) 

For  each  cold  beam,  the  relations  (4)  and  (5)  are 
then  modified  by  the  replacement 

0,o-+0  z0-kP,  (21) 

and  therefore  we  may  obtain  the  solution  for  each 
cold  beam  by  making  the  replacement,  in  Eq. 
(8), 


/0*o-»/(Ozo  -  kP).  (22) 

The  behavior  of  the  actual  warm  beam  will  then 
be  given  by 

=  ^(yowo )'m  f  dt'e-M™ + 

(23) 

where  the  average  is  defined  over  some  normal¬ 
ized  distribution  function  in  P,  i.e. 

<•••>„=  f  dPG(P).  .  .  .  (24) 

In  Eq.  (23)  we  can  immediately  anticipate  the 
effect  of  temperature  on  the  behavior  of  the 
beam;  the  entire  effect  is  included  in  the  phase 
factor,  in  the  term  kP.  Such  a  term,  when  aver¬ 
aged  over  any  reasonable  momentum  distribu¬ 
tion,  will  give  a  reduction  in  amplitude  of  the 
average  as  the  “width”  of  G(P)  is  increased. 
Physically  this  means  that  the  various  particles 
of  different  energies  within  the  beam  receive, 
when  passing  through  resonance,  displacements 
in  slightly  different  directions.  The  net  effect  on 
the  motion  of  the  beam  center  is  therefore  re¬ 
duced.  (Though  our  linearized  treatment  here 
necessarily  includes  a  fixed  beam  size,  it  may  in 
fact  be  the  case  that  a  warm  beam  will  just  expand 
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FIGURE  2  Rc  (Ai(i_  i)  vs.  time  for  7"/.  =  0.  FIGURE  4  Re  (Aiji- 1)  vs.  lime  for  Tj.  =  1.0. 


slightly  while  passing  through  resonance  while 
the  motion  of  the  beam  center  remains  relatively 
undisturbed.) 

As  an  example,  we  consider  a  beam  made  un 
of  particles  having  the  energy  distribution 

G, Av)  -  \ ]ITl  (25) 

GdAy)  ~  |  0  |  Ay  |  >  TJ2  ’ 

where  TL  is  a  measure  of  the  longitudinal  tem¬ 


perature  and  where  A-y  is  related  to  P  by  Eq. 
(19).  We  consider  again  Ihe  hypothetical  exper¬ 
iment  described  in  the  preceding  section.  For 
hf-\/Bc0  =  5  x  I0~3  the  results  of  a  numerical 
evaluation  of  Eq.  (23)  are  shown  in  Figs.  2-5, 
which  correspond  to  TL  =  (0.,  0.5,  1 .0,  2.0).  In 
each  figure,  the  real  part  of  Ai|;_ ,  in  centimeters 
is  plotted  versus  time  in  seconds.  The  resonance 
condition,  Eq.  (II),  is  satisfied  at  the  center  of 
the  time  axis.  Total  elapsed  time  is  2.1  p.sec.  The 
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chamber  diameter,  2a  =  20  cm  is  indicated  by 
solid  horizontal  lines  on  each  plot.  We  observe 
that  for  this  example,  TL  -  1.0,  or  a  0.5-McV 
energy  spread,  is  adequate  to  smooth  out  the  ef¬ 
fect  of  the  resonance.  This  is  the  same  order  of 
magnitude  of  spread  needed  to  damp  the  nega¬ 
tive-mass/kink  instability  in  this  device3. 


IV.  CONCLUSIONS 

We  have  obtained  a  bound  on  the  magnitude  of 
field  errors  that  can  be  tolerated  in  a  modified 
betatron  in  order  that  certain  integer  resonances 
may  be  safely  passed  through.  We  have  found 
that  for  practical  parameters  the  bound  is  ex¬ 
tremely  restrictive.  The  basic  difficulty  stems 
from  the  fact  that  unless  the  externa!  parameters 
of  the  system  are  changed  very  quickly,  the  orbits 
remain  in  or  near  resonance  for  many  betatron 
oscillations,  allowing  the  displacements  to  grow 
to  large  levels.  Such  a  result  suggests  that  non¬ 
linear  effects  may  play  an  important  role  in  beam 
behavior  near  a  resonance.  For  example,  one 
may  ask  whether  the  radial  dependence  of 
would  be  sufficient  to  “detune”  the  resonance 
as  the  beam  moves  a  finite  but  small  distance 
from  its  equilibrium  position.  This  possibility  is 
receiving  further  study. 

We  have  also  shown  that  a  finite  longitudinal 
beam  temperature  acts  to  reduce  the  effect  of  the 


resonance  on  the  motion  of  the  center  of  the 
beam.  The  temperature  spread  required  appears 
to  be  comparable,  in  a  specific  example,  to  that 
needed  to  stabilize  certain  microinstabilities. 5  It 
remains  unresolved  in  this  analysis  whether  the 
beam  expands  when  passing  through  a  reso¬ 
nance.  Such  behavior,  of  course,  if  severe,  could 
be  as  unacceptable  as  large  whole-beam  displace¬ 
ment. 

Should  it  be  possible  to  achieve  significantly 
lower  field  errors  than  those  used  in  our  example 
(0.5%),  or  if  it  is  possible  experimentally  to  detect 
and  correct  by  some  feedback  mechanism  the 
sudden,  resonant  displacement  of  the  beam,  then 
perhaps  lower  toroidal  fields  may  be  employed 
initially  and  be  removed  either  during  or  follow¬ 
ing  acceleration.  The  effects  of  passage  through 
the  low-/  resonances  may  thereby  be  reduced  to 
a  tolerably  small  level. 
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MODE  COUPLING  IN  THE  MODIFIED  BETATRON 

I.  Introduction 

For  small  displacements  about  a  planar  reference  orbit,  particle  motion 
transverse  to  the  toroidal  field  in  the  modified  betatron  may  be  represented  as  a 
linear  superposition  of  two  eigenmodes  of  motion:  a  "fast”  mode,  corresponding 
to  gyration  about  the  toroidal  field  lines  and  a  "slow"  mode,  corresponding  to  an 
drift  motion,  where  here  the  force  ?  is  due  to  the  weak  focusing  betatron 

■U 

fields,  space  charge  forces,  and  induced  (wall  image)  fields.  The  linear  theory 
of  orbits  in  the  modified  betatron  has  been  worked  out  in  some  detail^-'*  and  will 
only  be  reviewed  as  needed  here.  In  the  present  paper  we  will  mainly  discuss  the 
effect  of  quadratic  non-linearities  on  the  motion. 

Non-linear  terms  In  the  equations  of  motion  become  Important  to  consider 
if:  (1)  displacements  from  the  reference  orbit  become  large,  due  say  to  the 

method  of  injection  used  or  due  to  the  operation  of  an  Instability  of  some  kind, 
(2)  strong  non-linearities  (e.g.  large  values  of  9n/3r)  are  present  in  the  magnet 
design,  or  (3)  the  non-linear  term  itself  contains  a  resonant  part.  In  the 
following  we  will  illustrate  two  effects  of  quadratic  non-linearities  on  single 
particle  motion,  viz.  the  amplitude  dependence  of  the  betatron  frequencies  and 
the  exchange  of  energy  between  the  oscillation  modes  under  certain  conditions. 
These  conditions  turn  out  to  be  analogous  to  the  so-called  Walkinshaw  resonance^ 
in  accelerators  without  a  toroidal  magnetic  field^-^.  VJe  will  limit  ourselves 
here  to  consideration  of  single  particle  motion  only,  neglecting  the  effects  of 
self  fields;  the  treatment  here  then  will  only  be  valid  for  fairly  large  values 
of  y  In  high  current  devices,  such  that  v/y  «  1>  where  v  is  Budker's  parameter. 

Four  sections  follow.  In  the  first  we  introduce  some  notation  and  sketch 
the  derivation  of  the  equations  of  motion  to  second  order  in  displacements  from 
and  transverse  velocities  about  the  reference  orbit,  taken  to  be  a  circle  in  the 

symmetry  plane.  In  the  second  section  the  equations  of  motion  are  solved 
Manuscript  approved  February  15, 1983. 
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perturbatively  and  a  condition  for  the  generalized  Valkinshaw  resonance  is 
obtained.  Under  this  condition  we  study  the  behavior  of  the  betatron 
oscillations,  giving  a  numerical  example  as  illustration.  As  an  interesting 
result  we  find  a  particular  value  of  field  index  gradient  for  which  the  resonance 
ceases  to  have  any  effect  on  particle  motion. 

In  these  first  two  sections  the  discussion  will  assume  that  all  applied 
fields  are  azimuthally  symmetric.  In  the  presence  of  field  perturbations  other 
orbital  resonances  may  occur  and  it  is  interesting  to  ask  whether  the  amplitude 
dependence  of  the  betatron  frequencies  induced  by  the  non-linearities  are 
sufficient  to  keep  the  oscillation  amplitudes  at  finite,  but  tolerably  small 
values.  Though  In  general  this  is  a  difficult  question  we  will  discuss  a 
special,  simple  case  in  the  third  section  of  this  paper  in  which  the  Walkinshaw 
resonance  coincides  with  both  an  integer  and  half-integer  orbital  resonance. 

A  final  section  summarizes  these  results  and  states  some  conclusions  and 
conjectures. 
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II.  The  Equations  of  Motion 

The  geometry  of  the  modified  betatron  is  shown  in  Fig.  1.  We  employ 
standard  (r,0,z)  cylindrical  coordinates.  The  exact  equations  of  motion,  using 
0  in  favor  of  time  for  our  independent  variable,  are  written 


|g-  [r'(r'2  +  r2  +  z'2)"1/2]  =  r(r'2  +  r2  +  z'2)“1/2  +  A_1  (rBz  -  z'BQ)(l) 


d  ,  x,  .2  .  2  .2-1/2 

-rr  fz  (r  +  r  +  z'  ) 


)  J  =  A  (r'B0  -  rBr) 


2 

where  $  is  any  suitable  function  of  position  (r,0,z),  A  =  -me  By/e,  e  and  m  are 
the  magnitudes  of  the  electron  charge  (e  >  o)  and  mass,  0  and  y  are  the  usual 
relativistic  factors,  c  Is  the  speed  of  light  and  a  prime  {' )  denotes  d/d0. 

We  shall  assume  that  Br  vanishes  on  the  plane  z  =  0  and  take  all  fields  to 
be  independent  of  0.  (The  assumption  of  azimuthal  symmetry  will  be  relaxed  in 
Section  IV,  below.)  We  take  the  equilibrium  orbit  of  a  particle  of  relativistic 
factor  yQ  at  r  =  rQ,  z  =  0,  so 


A  =  -  r  B  (r  ,o). 
o  z  o 


Let  us  now  define  the  normalized  coordinates  x  =  (r  -  r  )/r  and  y  =  z/r  . 

o  o  o 

The  vector  potential  is  given  correctly  to  third  order  by 


.  v  l""11  2  n  2  2  2  1  .  ,  0  v  3. 

A0  *  roBzo11  +  —  x  +  2  y  “  —  xy  +  6<n  +  n2“3)xl 


where  Bzo»  n>  and  are  constants.  The  corresponding  fields  are 


B  “  -  B  [n  -  n„x]y 
r  zo  1  2 
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(6) 


~  n2  2  n  n2  2 
Bz  =  Bzo  f1  -  nx  +  ir  xz  +  —  y  ) 


from  which  n.  is  identified  as  the  second  radial  logarithmic  derivative  of  B  . 

^  z 


The  toroidal  field  is  assumed  to  be  given  by 


\  '  V(1+X)  "  B8o(1-X  +  '  '  •> 


where  Bqq  is  the  value  of  the  toroidal  field  at  the  reference  orbit,  x  =  y  =  0. 

Using  the  fields  (5,6,7)  in  the  equations  of  motion  (1,2)  and  keeping  terms 
only  of  quadratic  order  gives  the  coupled  equations: 


^  ^  9  T  9  9 

x~  +  (l-n)x  =  by"  +  (2n-l-  -  (— ~)y^  +  j  (x'  -  y'Z) 


y"  +  ny  =  -bx-  -  (Zn-n^xy  +  x' y' 


where  b  =  BQc/Bzo*  These  equations,  (8)  and  (9),  are  our  starting  points.  In 
the  following  section  we  examine  the  behavior  of  an  approximate  solution  to  (8) 
and  (9)  for  various  values  of  n,  n^,  and  b. 
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III.  Perturbative  Solution  of  Equations  of  Motion 


In  general  the  quadratic  terms  in  (8,9)  will  be  small  so  we  attempt  to  treat 
the  equations  perturbatively .  Neglecting  the  non-linear  terms  altogether  one  has 
the  solution  to  the  linear  equations: 


where  and  Ag  are  complex  numbers  depending  on  particle  initial  conditions  and 
the  frequencies  are  given  by 


.2.-  +  ra2An2  /  v,1/2  l/2 

yj  =  [ b  +1  ~  Kb  +1>  ~  ^(l-n)]  j 

s  ^ 


(ID 


The  subscripts  ?  and  s  are  used  here  and  below  to  label  the  amplitudes  and 
frequencies  of  the  fast  and  slow  oscillation  modes.  We  will  assume  that  the 
linear  motion  is  stable,  that  is  n(l-n)  >  0. 

We  may  calculate  the  correction  to  (10)  due  to  the  non-linear  terms  by 
inserting  (10)  in  (8,9)  and  resolving.  The  resulting  equations  will  be 
Inhomogeneous  with  various  "driving”  terms  at  the  frequencies 
2v_,  2v  ,  0,  and  v,  +  v  .  Consequently,  the  non-linear  correction  to  (10)  will 

t  S  JL  ~  S 

remain  small  unless  it  happens  that 


V 


f 


v  =  v 


s  s 


(12) 


the  condition  for  which,  from  (11),  being 


f  (n(l-n) ]  -1 . 


(13) 


In  the  absence  of  a  toroidal  field,  (13)  is  satisfied  for  n  =■  0.2  or  0.8  which  we 

identify  as  the  Walkinshaw  resonance^,  the  consequences  of  which  were  first 

12 

observed  in  cyclotrons  .  We  proceed  to  examine  particle  behavior  on  this 
resonance  in  the  modified  betatron. 

On  resonance  conventional  perturbation  theory  fails  and  one  must  resort  to 

some  other  method.  A  multiple  "time"  scale  analysis  of  the  problem  gives  a 

solution  of  the  form  (10)  in  which  the  complex  amplitudes  kc  and  A  are  no  longer 

Is 

strictly  constant  but  vary  slowly  with  0;  they  are  found  to  obey  the  equations 

dA  * 

3eT  '  lr1<n-Ws  <lA> 

dAf  2 

W  =  ir2(n’n2)As  (15> 


where  1*^  2  are  two  rea*  valued  functions  of  the  field  indices  n  and  n^. 


2  2 
v  -  n  „  2v  +  n 

T,  »  — — —  1 2n  +  6v  +  n9( — ^5 -  )] 


1  6v  3 
s 


2  2 
v  -  n 
s 


v  -  n 

-  -4  (-^ - )  r2 

vf  -  n 


(16) 


and  where,  on  resonance,  we  have 


2  «  2  =  ..1/2 

vf  “  4vg  =  2(n(l-n)) 


(17) 


(An  asterisk  denotes  complex  conjugate  in  (14).) 

The  question  of  orbital  stability  on  resonance  is  thus  reduced  to  the 
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question  of  the  behavior  of  the  mode  amplitudes  in  (14)  and  (15). 

The  equations  (14)  and  (15)  may  be  completely  solved  in  a  straightforward 
manner;  the  solution  is  obtained  and  discussed  in  the  Appendix.  To  settle  the 
stability  issue,  however,  it  is  sufficient  to  note  that  there  is  a  simple 
integral  of  motion 

J  |Af |2  +  1  r2  1asI2  =  o/r1,  (18) 

and  consequently  the  motion  is  necessarily  bounded  if  r  =  T ^  >_  0  which,  in 

fact  is  true  for  all  n  and  n^,  as  follows  from  (16).  On  this  resonance  energy  is 

simply  exchanged  back  and  forth  between  the  fast  and  slow  modes  of  motion. 

Though  we  have  argued  that  particle  motion  is  bounded  on  this  resonance  we 

have  not  in  fact  specified  a  bound  or  showed  that  the  bound  is  acceptable,  in 

terms  of  some  machine  aperture.  One  might  conjecture,  from  (18),  that  if  one  of 

2  were  significantly  larger  than  the  other  then  transfer  of  energy  from  the 

more  "stiff"  (larger  r  coefficient)  mode  to  the  less  "stiff”  mode  would  result  in 

increasing  particle  oscillation  amplitude.  Hence  one  would  be  concerned  if,  from 
2  2 

(16),  either  v  ■  n  or  vf  =  n.  It  follows  from  (17)  and  (18),  however,  that 

S  JL 

this  can  occur  only  for  b  =  0,  n  a  .2  or  .8.  If  n  is  chosen  so  that  b  is 
0(1)  when  the  resonance  is  crossed  then  and  r2  are  of  the  same  order  of 
magnitude  and  one  expects  this  resonance  to  be  quite  harmless.  For  specific 
initial  conditions  ft  is  possible  to  find  a  bound  by  calculating  the  turning 
point  of  a  certain  particle-in-a-well  problem,  as  shown  in  the  Appendix. 

Curiously,  one  can  render  this  resonance  completely  inoperative  for  any 
particular  n  by  choosing  n2  so  that  and  ^  both  vanish.  From  (16)  and  (17) 
this  value  is  found  to  be 
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n2  E 


r7n  -  3  +  4[(l-n)n]1/2 
1  +  [(l-n)/n]1/2 


] 


(19) 


Choosing  this  value  ensures  that  the  mode  amplitudes  remain  constant  when  passing 
through  the  "exchange"  resonance. 

We  proceed  to  illustrate  some  of  these  results  using  a  simple  single 
particle  numerical  orbit  integration.  The  algorithm  includes  the  fields  (5-7) 
but  does  not  use  an  expansion  of  the  force  or  acceleration.  Figures  2  and  3  show 
the  solutions  to  equations  (14-15)  and  (1-2)  respectively  for  the  case  n  **  0.5, 

A 

n2  =  n2  “  0.625,  b  51  0.5.  The  mode  amplitudes  are  strictly  constant,  no  exchange 
occurs,  and  the  particle  orbit  projection  retraces  itself  in  a  stable  manner  over 
and  over  again.  We  contrast  this  case  with  that  illustrated  in  Figures  4  and  5 
for  which  the  parameters  are  again  n  =  0.5  and  b  =*  0.5  but  now  with 
Uj  =  0.  Now  the  mode  amplitudes  oscillate;  one  rises  while  the  other  falls  in 
order  to  conserve  D  (Eq  (18)).  The  oscillation  period,  from  Equation  (16)  in  the 
Appendix,  is  46.6  major  periods  for  the  particular  initial  conditions  chosen. 

The  particle  orbit  projection  now  simply  (and  harmlessly)  rotates  slowly  counter¬ 
clockwise. 

Our  conclusion  is  that  in  the  case  of  azimuthally  symmetric  fields  the 
generalized  Walkinshaw  resonance  is  quite  harmless  in  the  modified  betatron.  The 
exchange  of  energy  between  fast  and  slow  modes  is  expected  to  cause  no  major 
changes  in  the  beam  dimensions.  When  azimuthal  field  variations  are  present, 
however,  the  situation  changes  dramatically  due  to  a  coincidence  described  in  the 
next  section. 


J 

J 


■i 


tv  i 


1 


.1 


►•J 


36.11 


kv 


,;>■>  m  %  kl  m  I  mJt 


IV.  A  Triple  Coincidence  Resonance 

When  n  =  1/2  (the  case  illustrated  in  Figures  2-5)  the  values  of  the  tunes 

V-  and  v  at  the  exchange  resonance  (13)  are  v,  =  1,  v  =  1/2;  therefore  in  the 
X  S  X  s 

presence  of  field  and  focusing  errors  the  generalized  Walkinshaw  resonance 
coincides  with  an  integer  and  a  half  integer  orbital  resonance.  This  triple 
coincidence  allows  us  to  study  in  detail  in  this  special  case  the  effect  of  mode 
coupling  (and  the  amplitude  dependence  of  the  betatron  frequencies)  on  the  orbit 
at  an  integer  and  half  integer  resonance.  Though  the  restriction  to  n  =  1/2  is 
necessary  for  there  to  be  a  true  coincidence,  if  n  is  near  but  not  exactly  1/2, 
the  three  resonances  will  be  ‘’nearby"  and  will  occur  nearly  simultaneously  and 
the  analysis  below  should  still  hold  in  an  approximate  way. 

At  the  triple  coincidence  resonance  the  mode  evolution  equations  become 

^As  *  * 

ar  =  lriAfAs  +  eias  <20> 

de”  =  lr2As  +  e2  +  e3Af  *  (21) 

The  e's  in  (20)  and  (21)  are  complex  constants  proportional  to  certain 
Fourier  coefficients  in  expansions  of  the  fields  and  their  gradients; 
specifically,  is  due  to  an  l  =  1  term  in  the  field  gradient,  leading  to  a  half 
integer  resonance,  is  due  to  an  %  -  1  term  in  the  field,  leading  to  an  integer 
resonance,  and  is  due  to  an  l  =  2  term  in  the  field  gradient,  leading  to  a  ”2- 
halves"  integer  resonance.  Were  they  present  alone  (i.e.  with  no  non-linearity) 
in  (20)  and  (21)  the  field  imperfection  terms  are  observed  to  lead  to  the  usual 
linear  and/or  exponential  growth  characteristic  of  integer  and/or  half  integer 
resonances.  In  the  presence  of  mode  coupling  the  situation  is  much  less  clear. 
Since,  as  they  stand,  equations  (20)  and  (21)  cannot  be  solved  analytically  we 


must  in  general  resort  to  numerical  integration.  First  let  us  comment  on  what  we 
might  expect  to  see  in  the  solution. 

The  field  imperfection  terms  in  (20)  and  (21)  act  roughly  speaking  as  source 
terms,  pumping  energy  from  longitudinal  to  transverse  motion.  If  this  energy 
flow  continues,  and  if  there  is  no  mechanism  to  return  this  energy  to 
longitudinal  motion  the  result  is  disastrous  —  a  linear  orbital  resonance.  Non- 
linearities,  however,  can  shift  the  betatron  frequency  of  the  resonant  particle 
off  resonance  for  some  finite  amplitude  of  betatron  oscillation.  (The  quantity 

A 

(n^  -  n^)  is  presumably  a  measure  of  the  frequency  shift  induced  by  a  given 
amplitude  oscillation.)  In  practice  though  one  can  not  say  a  priori  how  much 
frequency  shift  will  be  sufficient  to  terminate  the  growth  of  the  resonant  mode. 

A  numerical  study  seems  to  be  essential. 

For  a  numerical  example  we  will  examine  the  effect  of  the  non-linear  terms 
in  (20)  and  (21)  on  an  integer  resonance,  that  is,  in  the  following  we  shall  take 
=  0.  Cases  have  been  examined  numerically  for  various  other  combinations 
of  values  for  the  e's  with  no  major  differences  appearing  in  the  results. 

In  Figures  6  and  7  we  illustrate  the  mode  amplitudes  and  orbit  projection 

A 

for  a  pure  integer  resonance  with  no  mode  coupling  (n^  =  n^  =  0.625, 
e 2  "  .005).  The  (resonant)  fast  mode  amplitude  grows  linearly  without  limit;  the 
(decoupled)  slow  mode  stays  at  a  fixed,  small  value.  The  particle  orbit  size 
(Figure  7)  consequently  grows  continuously. 

Turning  on  the  mode  coupling  changes  the  behavior  of  the  mode  amplitudes 
dramatically  but  has  little  apparent  effect  on  the  particle  orbit,  which  still 
appears  to  grow  to  intolerable  size.  This  is  illustrated  in  Figures  8  and  9 
where  we  see  that  the  mode  amplitudes  grow  to  a  certain  size  and  then  turn  over  - 
-  presumably  a  reflection  of  the  detuning  of  the  resonance  due  to  the  frequency 
shift.  The  "turnover",  however,  is  at  extremely  large  amplitudes  (Recall  that 
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the  mode  amplitudes  are  normalized  to  the  radius  of  the  device,  rQ.),  therefore 
the  particle  motion  appears  to  be  relatively  unaffected,  practically  speaking,  by 
changing  nj  from  5/8  to  0.  This  result  does  suggest  though  that  by 

a 

increasing  In^  -  02!  we  might  reduce  the  resonant  response  to  a  tolerable  value. 

Figures  10-11,  12-13,  and  14-15  show  our  results  for  n2  =  -1.,  -4.,  and 

A 

-10.  respectively.  We  see  that  as  ^2  ”  n2^  *s  increased  the  mode  turnover 
amplitude  is  reduced  and  the  particle  orbit  becomes  somewhat  more  compact, 
staying  within  the  plot  boundaries  for  significantly  longer  times.  (Even  so,  the 

A 

transverse  orbit  size  is  rather  large,  even  for  the  largest  values  of  |n2  -  n2l 
we  have  tried.) 

These  results  suggest  that  to  stabilize  the  l  m  1  resonance  a  significant 
non-linearity  (large  value  of  02)  could  be  intentionally  introduced  in  the 
betatron  field.  One  must  be  careful  in  drawing  this  conclusion,  however,  because 
such  a  non-linearity  has  well  known  adverse  effects,  among  them  a  sensitivity  of 
the  behavior  of  the  orbit  to  initial  conditions;  that  is,  only  some  special  class 
of  particles  may  be  confined  while  others  are  lost.  Also,  if  which  is 
effectively  the  radial  derivative  of  n,  is  very  large,  it  then  becomes  difficult 
to  keep  n  itself  within  the  stable  range  0  n  1  everywhere  within  the 
aperture.  Consequently  we  conclude  that,  as  a  practical  matter,  it  is  best  not 
to  rely  on  non-linearities  to  stabilize  the  l  =  1  resonance  in  the  modified 
betatron  and  to  design  the  machine  with  a  flat  radial  index  profile.  Avoidance 
of  this  resonance  as  well  as  other  low  order  resonances  -  -  which  then  becomes 
the  only  reasonable  experimental  alternative  -  -  is  possible  in  principle  by 
accelerating  with  constant  b  (i.e.  B^o  «  B^),  thereby  keeping  the  tunes  fixed** 

-  -  except  for  the  tune  shift  due  to  space  charge,  which  affects  the  fast  mode 
tune  only  very  slightly  for  large  B^;  the  slow  mode  tune  can  generally  be  chosen 
to  be  very  small  (~  .2  -  .3)  for  all  time. 
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V.  Summary  and  Conclusions 

We  have  examined  Che  effects  of  mode  coupling  on  single  particle  orbits  in 
the  modified  betatron.  We  find  that  a  generalization  of  the  Walkinshaw 
(exchange)  resonance  can  occur  for  any  value  of  field  index  in  the  range 
.2  n  .8  but  that  its  effect  on  particle  orbits  in  general  is  quite  modest 
(Figures  2—5)  and  may  be  rendered  completely  ineffective  by  a  special  choice  of 
field  index  gradient  (19). 

When  n  is  near  1/2  the  exchange  resonance  coincides  with  both  integer  and 
half  integer  resonances.  An  examination  of  orbit  behavior  at  this  triple 
coincidence  shows  that,  as  a  practical  matter,  the  amplitude  dependent  frequency 
shift  in  the  betatron  oscillation  due  to  mode  coupling  is  not  sufficient  to 
stabilize  the  t  •  1  integer  resonance  (though  presumably,  as  in  the  case  of 
accelerators  not  employing  toroidal  fields,  higher  order  resonances  will  be 

O 

subject  to  non-linear  stabilization  ).  This  fact  makes  it  advisable  to  allow,  in 
the  design  of  an  experiment,  for  acceleration  with  constant  or  nearly  constant 
ratio  3gQ/3zo  thereby  holding  the  tunes  approximately  fixed  in  time. 
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Appendix 

In  this  Appendix  we  discuss  the  solution  to  Equations  (14)  and  (15)  in  the 


text.  Writing 


Af  -  afexp(i<j>f)  Ag=>  asexp(i<j>s) 


(Al) 


where  af  s  and  (j>,  are  real  we  find  in  a  straightforward  way  an  equation 
*•  >**  t  ) s 


r  1  2 

for  P  =  2  as  : 


i  p'2  +  V(p)  =  0 


(A2) 


where 


v(P)  =  4rP3  -  4 Dp 2  +\c2 


r  -  ri  r2 


C  =  r.  aca  cos(4£-  26  )  =  constant 
1  f  s  Yf  rs 


and  the  constant  D  is  defined  in  the  text,  Eq.  (18).  The  other  quantities  are 
given  in  terms  of  p  by 


1/2 

ag  -  (2P)i/Z  (A3) 

-  C/(2p )  (A4) 

Af  -  (C  -  ip" )  exp(2i$g).  (A5) 


These  expressions  hold  for  p  >  0.  If  p  »  0  at  0  ■  0,  say,  then,  from  (14) 

and  (15)  A  =  0  and  A,  remains  fixed  for  all  950.  The  modal  frequency  shifts 
s  f 

are  given  directly  by  (A4)  for  the  slow  mode  and  may  be  obtained  from  (A5)  for 
the  fast  mode. 

It  may  be  shown  that  V(p)  has  one  negative  and  two  positive  roots.  Denoting 
these  by  p^  ^  3  with  p^  >  pj  >  0  >  p^  we  find  that  the  exchange  period  (period  of 
p(9))  is  given  by 


1/2 


K(m) 


r(P]_  "  p3>' 


(A6) 


where  m  =  (p^  -  P2^^pi  ~  P3)  atl(*  w^ere  K  is  the  usual  complete  elliptic 
integral. 

The  special  cases  C  =  0  and  T  =  0  lead  to  motion  of  infinite  period  and 


p  =  constant,  respectively. 
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10.  20.  50.  40.  50 

MAJOR  PERIODS 

Fig.  2:  |Af | ,  | |  vs.  major  periods  (A  major  period  is  a  change  of  6  by 

2ir.)  n  *  1/2;  n ^  =  5/8.  Initial  (0  =  0)  values  of 
-3  * 

A^=  5.7735  x  10  (1  -  i)  and  Ag  *  2A^  correspond  to  those  of  a 

particle  initially  at  x  **  y  «  .02  with  zero  transverse  velocity, 
same  initial  values  are  used  in  all  subsequent  figures. 
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z  vs.  r  -  rQ;  n  -  1/2;  n2  =0;  5B  =  0.5  gauss.  Particle  leaves 
area  after  ~  10  major  periods. 
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Fig.  13:  z  vs.  r  -  tq;  n  *»  1/2;  n2  =  -4;  6B  =  0.5  gauss.  Particle  remain 
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MAJOR  PERIDDS 


Fig.  14:  |Aj|,  |Ag|  vs.  major  periods;  n  =  1/2;  n 2  =* 

e2  -  5  x  10“3. 
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1/2;  ^  =*  -10;  6B  =  0.5  gauss.  Particle  remains  in 
major  periods. 
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Abstract 


Current  thresholds  for  longitudinal  and  transverse  instabilities  are 

calculated  for  beams  of  specified  dimensions  in  conventional  (B  =  0)  and 

0 

modified  0)  betatrons,  using  simple  models  for  the  longitudinal  and 

transverse  impedances.  Self  field  effects  of  the  beam  are  included  and  lead 
to  a  novel,  competitive  effect  between  the  stabilization  mechanism  and 
instability  growth.  This  competition  results  in  a  multi-valuedness  in  the 
limiting  current  vs  beam  energy  spread  plot,  even  for  conventional 
betatrons.  Accessibility  of  the  various  limiting  current  levels  appears  to 
depend  upon  the  rate  of  beam  injection.  The  stabilizing  effects  of  betatron 
oscillations  are  discussed  and  written  as  the  sum  of  three,  physically 
interpretable  contributions  to  an  effective  energy  spread.  We  find  that  the 
presence  of  a  strong  toroidal  field  can  significantly  improve  the  current 
carrying  capacity  of  the  accelerator. 


I.  Introduction' 

The  addition  of  a  toroidal  magnetic  field  to  a  conventional  betatron  has 

been  shown  theoretically^  **  to  increase  the  equilibrium  current  which  may  be 

2 

confined  by  a  factor  1/2  (B  /B  )  ,  for  large  values  of  B./B  .  In  conventional 

O  Z  0  z 

accelerators  without  solenoidal  focusing,  however,  beam  stability  considerations 

6— ft 

place  the  actual  limit  on  beam  current.  °  Therefore  it  becomes  important  to 

analyze  the  stability  conditions  and  associated  limiting  currents  for  a  given 

beam  equilibrium  in  the  presence  of  a  toroidal  field.  In  this  paper  we  present 

such  an  analysis  for  both  longitudinal  and  transverse  modes. 

A  device  in  which  a  toroidal  magnetic  field  is  superimposed  on  the  usual 

weak  focusing  betatron  field  has  come  to  be  called  a  "Modified  Betatron”.  See 

Figure  1.  A  stability  analysis  of  this  accelerator  necessarily  must  include  the 

strong  self  and  induced  (wall  image)  fields  of  the  electron  beam.  It  is 

primarily  the  inclusion  of  these  fields  which  distinguishes  this  work  from  the 

o 

stability  analysis  performed  for  the  so-called  plasma  betatron  in  which  self 
fields  are  much  less  important.  These  self  field  effects,  however,  will  be  seen 
to  have  a  dramatic  effect  on  the  current  versus  energy  spread  scaling;  namely  we 
predict  the  existence  of  more  than  one  stable  value  of  current  for  a  given  beam 
energy  spread.  This  somewhat  surprising  result  will  be  discussed  later.  The 
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central  result  of  this  work,  however,  is  that  significantly  more  current  may  be 
carried  by  a  beam  in  a  modified  betatron  configuration  than  in  a  conventional 
betatron  of  the  same  dimensions,  assuming  equal  beam  sizes  and  energy  spreads. 

In  reaching  this  conclusion  we  have  included  not  only  the  effects  of  self  fields, 
but  also  the  (stabilizing)  effects  of  betatron  oscillations  and  the 
(destabilizing)  effects  of  short  wavelength  enhancements  to  the  longitudinal  and 
transverse  impedances  due  to  chamber  resonances.  Below  we  discuss  the  dispersion 
relation  for  arbitrary  toroidal  fields  and  currents,  describe  our  model  for  the 
impedances,  and  present  analytical  and  numerical  results  from  the  dispersion 

relation. 


II.  Discussion 

A  dispersion  relationship  for  both  the  longitudinal  and  transverse  modes  in 
a  modified  betatron  accelerator  configuration  has  been  derived.*®  Included  in 
the  derivation  are:  beam  self  field  effects,  induced  field  effects  arising  from 
wall  image  charges  and  currents  as  well  as  finite  chamber  wall  conductivity 
effects.  Toroidal  corrections  to  the  equilibrium  beam  self  fields  and  chamber 
wall  image  fields  have  been  neglected.  The  longitudinal  and  transverse 
impedances,  which  characterize  the  beam  environment,  are  incorporated  in  a 
phenomenological  way  in  the  short  wavelength  limit.  The  dispersion  relation, 


therefore,  treats  disturbances  of  all  wavelengths,  including  wavelengths  much 

longer  or  shorter  than  the . chamber  minor  radius. 

The  short  wavelength  model  for  the  impedances  contains  effects  associated 

with  propagating  chamber  modes.  These  effects  can  significantly  affect  the 

instability  growth  rates.  With  the  inclusion  of  the  short  wavelength 

contributions  to  the  impedances  a  realistic  estimate  for  the  current  limitations 

due  to  the  various  instabilities  can  be  obtained  for  the  modified  betatron 

configuration  and  compared  to  those  of  a  conventional  betatron.  To  perform  a 
meaningful  comparison  we  will  choose  identical  parameters  for  the  two  types  of 
betatrons,  i.e.,  same  geometry,  injection  energy,  field  index,  etc.  The  only 
difference  of  course  will  be. that  the  modified  betatron  configuration  will 
include  a  toroidal  magnetic  field. 

The  dispersion  relation*®  for  the  longitudinal  and  transverse  modes  of  a 
cold  beam  may  be  written 


where  (w  -  £.toc)/uc,  u  is  the  complex  mode  frequency,  £.  =  1,2, 3,... is  the 

longitudinal  (toroidal)  harmonic  mode  number,  o>  «  ft  /y  is  the  electron  rotation 


frequency,  ftz  -  lelB^/ra^c  is  the  non-relativistic  cyclotron  frequency,  y  = 

22—1/2  ~  2  2 
(1  -  /c  )  ,  vq  is  the  longitudinal  (toroidal)  beam  velocity,  *»  -  2i£ 

(v/y)(Z B/iZo),  v  is  Budker's  parameter(v  =  I[A]/17  x  103),  Z (  =  Z ^(u)  is  the 

total  effective  longitudinal  impedance,  ZQ  =  Art/c  (ZQ  =  377  ft  in  MKS 
~  2  3  2 

units),  =  1/2  -  2( v/y  )(ro/rfe)  “  2i(v/Y)(rQ  Z  /Zo),  rQ  is  the  major  electron 
beam  radius,  r^  is  the  minor  electron  beam  radius,  Z^  =  Z^Cw)  is  the  total 
effective  transverse  impedance,  b  =>  B  /B  and  B  is  the  toroidal  magnetic 

v  Z  0 

field.  In  (1)  finite  amplitude  betatron  oscillations  were  neglected,  the 

external  field  Index  was  taken  to  be  1/2  and  the  electron  beam  was  assumed  to  be 

mono-energetic,  highly  relativistic  (vq»  c)  and  circular  in  cross  section. 

III.  Approximate  Representation  of  Impedances 

In  our  model  the  longitudinal  impedance  Z|(u)  is  taken  to  consist  of  three 

terms 


z  =  Z  +  z  +  z 
I  I , s  1,0  t  ,r 


(2) 


The  first  term  in  (2)  is  the  long  wavelength  space  charge  shielding  contribution 
associated  with  a  smooth  infinitely  conducting  chamber  and  is  given  by 
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(2a) 


Z|,'s“  1  Zo(i/2*w  )(1-+  2  *n.a/rb>» 


where  y  is  the  relativistic  factor  corresponding  to  the  wave  phase  velocity  and 
w 

2 

a  is  the  minor  radius  of  the  toroidal  chamber.  Due  to  the  l/y  factor  this 

w 

impedance  term  is  typically  quite  small;  it  may  even  change  sign. 

The  part  of  the  longitudinal  impedance  due  to  the  resistive  nature  of  the 


chamber  wall,  in  the  long  wavelength  regime,  is 


Z.  -  £  Z  (1  -  i)  fi/2a. 
1,0  o 


“1/2 

where  6  =  c  (2Trolwl)  is  the  skin  depth  and  o  is  the  wall  conductivity.  It 
has  been  assumed  in  (2b)  that  the  skin  depth  is  small  compared  to  the  thickness 


of  the  chamber  wall. 

Finally,  the  last  term  in  (2)  represents  the  resonance  contribution  to  the 

total  longitudinal  impedance  and  arises  from  the  fact  that  the  chamber  can 

support  propagating  waves.  To  obtain  the  exact  form  for  Z  would  require  a 

i ,  r 

rather  involved  analysis  of  the  beam-chamber  structure  and  is  beyond  the  goals  of 
the  present  paper.  We  will,  therefore,  simply  represent  this  contribution  by  the 


phenomenological  expression 
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(2c) 


f-  ■ 

zl. 

I-’  .*■ 

K:  • 


,  w/u  -  i  Q(1  ■  (u/u_)  .)  % 
Z.  =  R  —  ( - = - - - - 

,,T  °  Wr  (w/u>  )2  +  QZ<1  -  (w/(o  )  V 


where  RQ  defines  the  chamber  shunt  impedance,  u>r  a  2.4  c/a  is  the  cutoff 
frequency  of  the  lowest  order  chamber  mode,  and  Q  is  the  quality  factor 
associated  with  the  chamber.  The  chamber  shunt  impedance,  RQ,  can  be  estimated 
by  noting  that  near  a  resonant  frequency  tu  *=  tur  the  longitudinal  impedance  is 
roughly  equal  to  the  free  space  impedance 


(/3  -  *)  * 


1/3 


(3) 


% 


1/3 

It  follows,  therefore,  that  R^  «  Z  i.  where  l  is  the  toroidal  mode  number 

o  o  r  r 

associated  with  the  resonant  frequency,  i.e.,  =  u>r/wc. 

The  expression  for  the  total  transverse  impedance  Z^(w)  is  also  written  as 


the  sum  of  three  contributions 


4 
I 


Z,=Z,  +  z,  +  z  . 

X  XjS  1)0  X|T 


(A) 


The  long  wavelength  space  charge  part  of  the  impedance,  for  a  smooth  infinitely 


Sr 

p  V ' 


conducting  chamber,  is 
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Zl,s  ■  1  roV1/rb2  -  <4»> 

The  second  term  is  the  long  wavelength  contribution  due  to  the  resistive  nature 
of  the  wall  and  is  given  by 


Z,  -  r  Z  (1  -  i)  6/a  . 
1,0  o  o 


Chamber  resonances  contribute  to  the  transverse  impedance  a  part  which  we  will 


represent  by  the  form1 


ZjL  r  =2  Z»  r’ 

l,r  £a* 


where  Z  is  defined  in  (2c). 
i  ,r 


In  the  present  work  we  will  not  be  concerned  with  resistive  wall  effects 


which  lead,  in  conventional  acclerators,  to  well  known  longitudinal  and 


transverse  instabilities0  having,  however,  comparatively  slow  growth  rates. 


Equations  (2b)  and  (4b)  are  included  here  only  for  reference  purposes. 


We  remark  that  with  the  impedances  as  defined  above  the  dispersion  relation 

Is  virtually  independent  of  the  beam  minor  radius,  r^,  except  for  the  weak 
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logarithmic  dependence  in  (2a).  This. is  reasonable  if  we  think  of  the  dynamics 
involved  in  the  various  instabilities,  that  is,  if  we  consider  that  it  is  the 
beam  centroid  which  moves  transversely,  even  in  the  "longitudinal"  or  negative 
mass  instability  in  which  beam  bunching  occurs  as  the  beam  centroid  moves  in  or 
out  radially.  Motion  of  the  beam  centroid  is  affected  by  the  externally  applied 
fields,  including  those  due  to  wall  image  charges  and  currents;  these  fields, 
unlike  beam  self  fields  which  are  carried  along  transversely  by  the  moving  beam,, 
do  not  depend  on  the  minor  radius  of  the  beam. 

IV.  Stability  Condition  and  Limiting  Current 

To  obtain  the  limiting  current,  based  on  stability  requirements,  for  the 

modified  and  conventional  betatron  a  stability  criterion  is  needed.  If  the 

distribution  in  particle  rotation  frequencies  is  Lorentzian  in  shape  the 

criterion  for  stability  is  simply 

r  <  JtlASlI,  (5) 

where  T  is  the  growth  rate  in  the  absence  of  a  frequency  spread  and  Afi  is  the 
half  width  of  the  frequency  spread  on  the  beam.  It  should  be  noted  that  the 


large  tails  associated  with  a  Lorentzian  distribution  make  the  criterion  in  (5) 

somewhat  less  stringent  than  would  be  the  case  if  a  more  realistic  choice  of 

frequency  distribution  were  used.  However,  the  use  of  a  more  realistic 

distribution  of  particle  frequencies  would  result  in  a  considerably  more  involved 

stability  criterion.  Since  we  are  interested  here  mainly  in  making  a  comparison 

of  the  limiting  currents  of  the  modified  and  conventional  betatrons,  consistent 

use  of  a  Lorentzian  for  both  devices  should  serve  our  purpose. 

Both  an  intrinsic  longitudinal  energy  spread  on  the  beam  electrons  as  well 

as  finite  amplitude  betatron  oscillations  will  produce  a  spread  in  rotation 

frequencies.  Solving  the  particle  orbit  equations  correct  to  second  order  in  the 

betatron  oscillation  amplitude,  with  self-field  effects  and  intrinsic  energy 

spread  effects  included,  the  average  spread  in  the  beam's  rotation  frequency  is 

found  to  be  given  by 

Iasi!  -  UnlAE  +  I ah! B,  (6) 


where 

i 


I  Aftl  ajj  -  j  wc  lal(AE/E) 


(6a) 
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(6b) 


IaqL  -  u  <r  A2r  r)bZ/2-n  -  3/2  j, 

B  C  D  O  S 

-1  -2 

with  a  =  (1/2  “  n  )  -  y  *  The  two  terms  on  the  right  hand  side  of  (6)  are  due 

s 

to  an  intrinsic  energy  spread  and  to  finite  amplitude  betatron  oscillations, 

respectively.  In  (6a)  the  fractional  intrinsic  energy  spread  is  denoted 

2 

by  A E/E  where  E  =  ymoC  is  the  beam  particle  energy,  (y  »  1)  and 
3  2 

n  =  2(v/y  )(r  /r  )  is  the  field  index  associated  with  the  beam's  self  fields. 

S  O  D 

In  obtaining  (6a, b)  we  assumed  a  circular  beam  cross  section  and  an  external 

field  index  of  1/2.  Note  that  it  is  here,  in  the  relation  between 

AE/E  and  AJ2,  that  self  fields  play  an  important  role. 

By  utilizing  the  beam  envelope  equation,  the  frequency  spread  term  Ia^L  cai 

be  expressed  in  a  more  Illuminating  form.  The  condition  for  a  matched  beam,  i.e 

12 

non-oscillating  minor  beam  radius,  is 

X  »2  "  "s-  l'2  +  <‘oS.>2'<rb2*)2>  (7> 

where  e  is  the  normalized  transverse  beam  emittance  as  measured  in  the  Larmor 

n 

frame.  Substituting  (7)  into  (6b)  gives 


In  (8),  the  first  term  in  brackets  is  the  familiar  longitudinal  energy  spread  due 

to  emittance,  the  second  term  is  a  toroidal  correction  to  the  first,  while  the 

last  term  is  the  energy  spread  associated  with  the  electrostatic  potential  drop 

across  the  beam.  Both  contributions  to  the  total  frequency  spread, 

|AnlAE  and  I  Aft  I B,  are  proportional  to  the  various  energy  spreads.  The  two 

different  proportionality  factors,  lal  in  the  case  of  } A£2 1 and  y  in  the  case 

of  I  Aft  I B,  arise  because  the  intrinsic  particle  energy  spread  produces  a 

rotational  frequency  spread  primarily  by  changing  the  particle's  radial  position 

whereas  the  various  energy  spreads  contributing  to  lAftlg  merely  result  in  a 

longitudinal  velocity  spread.  Hence  the  various  contributions  to  the  longitudinal 

energy  spread  contribute  differently  to  the  frequency  spread. 

The  desired  stability  criterion  is  obtained  by  substituting  (6)  into  (5)  and 

becomes 

r  <  i«c  [~  lal  ( AE/E)  +  Uen/rb)2/2  -  f  (r^r^V  +  v/2y|/Y2]  (9) 

Given  the  intrinsic  energy  spread,  beam  radius  and  emittance,  the  criterion  in 


(9)  implies  a  limiting  beam  current  which  if  exceeded  will  result  in  instability. 
As  a  simple  illustration  we  will  first  consider  the  negative  mass  mode  in  a  low 


current  conventional  betatron,  i.e.,  B„  =  0,  n  «  1/2.  From  (1),  the  dispersion 

0  s 

~  2  2  2 

relation  for  the  negative  mass  instability  is  Am  =  -  to  /u. 

X  U  X 

«  4i£ Z(v/y)(Z  u(w)/i.Zo),  where  we  have  assumed  that  lArn^l  «  Im^l  »  1//2. 
Approximating  Z^w)  by  Z((Jl(i)c)  we  find  that  the  growth  rate  is 

T  =  2£u)c(v/Y/^^  |  (Z  ^(£0^))  |  *^/(£Zq)^^.  Using  (9)  and  neglecting  the  betatron 
oscillations  we  recover  the  well  known  negative  mass  stability  condition 


Y  “ 


o _  rt E^2 

^T)T  eJ  • 


In  obtaining  (10)  we  have  assumed  that  the  real  and  imaginary  parts  of  Z|  are 
approximately  equal. 

Next  we  consider  the  full  dispersion  relation  (1)  for  arbitrary  b.  Here  and 
below  we  will  continue  to  neglect  the  effect  of  betatron  oscillations  on  the 


stability  of  the  beam.  The  effect  is  generally  small  and,  due  to  the  bz 
dependence  in  (6b),  it  favors  the  modified  betatron.  Therefore,  neglect  of  the 


betatron  oscillations  is  conservative  when  comparing  the  modified  and 


conventional  betatrons. 
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If  we  may  continue  to  approximate  Z „  (w)  by  Z  •  (Zw  )  and  if  we 

i,r  l,r  c 

set  6=0  then  (1)  is  in  general  a  sixth  order  (fourth  order,  if  b  =  0) 
polynomial.  By  neglecting  resistive  wall  effects  we  are  assuming  that  these  are 
negligible  for  the  negative  mass  branch  of  the  dispersion  relation  in  which  we 
will  be  interested  here.  Of  the  six  roots,  two  pair  correspond  to  transverse 
modes  (a  pair  being  a  forward  and  backward  wave,  essentially)  and  one  pair  to  the 
longitudinal  mode. 

For  a  given  set  of  parameters,  we  have  found  the  roots  numerically,  for  a 
large  range  of  Z;  a  maximum  value  of  T/Z  is  then  found  and  substituted  into  the 
stability  criterion  from  which  a  value  of  AE/E  is  computed.  An  important  point 
to  note  is  that  the  toroidal  field  has  more  of  a  stabilizing  effect  on  the 
high  Z  modes  than  on  the  low  Z  modes;  consequently,  as  b  is  increased  the 
behavior  of  the  impedance  for  small  Z  tends  to  determine  the  maximum  r /Z  and 
therefore  the  limiting  current. 

Typical  results  are  illustrated  in  Figures  2  and  3.  Here  we  plot  the  beam 
current  vs  the  Lorentzian  full  width  energy  spread  required  for  stable  motion  for 
b  =  0,  5,  10,  20,  rQ/a  =  6,  a/r^  =  3,  RQ/Z0  =  4,  and  Q  =  10.  Figures  2  and  3  are 
plotted  for  y  =  3  and  6  respectively.  The  dashed  lines  are  plots  of 
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For  consistency  we  must  be  to  the  left  of  these  lines.  This  restriction  may  be 
understood  by  considering  the  displacement  of  a  particle  from  the  center  of  the 

beam: 


6r 

r 

o 


6E/E 

1 


+  betatron  oscillations 


where  6E  is  the  difference  between  the  particular  particle  energy  and  the  average 
energy  of  the  particles  in  the  beam.  It  follows  that 


rb  .  | AE/E  , 
r  -  l-2n  1 


(12) 


o  s 

where  AE  is  the  full  width  of  the  distribution. 

The  effects  of  the  self  fields,  as  represented  by  n  In  (6a)  and  (11)  are 

s 

immediately  evident  in  the  plots  of  Figures  2  and  3.  If  ng  were  zero  then  the 
dashed  boundary  lines  would  be  a  single  vertical  line  and  the  curves  would  be 
monotonic.  As  it  is,  the  effect  of  the  self  fields  is  basically  traceable  to  the 
increasing  (as  ng  increases  toward  1/2)  then  decreasing  (as  ng  increases  beyond 
1/2)  factor  lal*  The  multi-valuedness  of  the  curves  may  then  be  understood  as 


follows  (Refer  to  Figure  4.):  For  very  small  currents  (Branch  I,  in  Figure  4) 
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the  cold  beam  growth  rate  is  small  and  an  increasing  function  of  current.  The 

self  field  index,  ng,  is  negligible  compared  to  1/2  and  the  energy  spread 

required  for  stability  is  an  increasing  function  of  current.  This  is  the  regime 

in  which  virtually  all  conventional  accelerators  operate.  As  the  current  is 

further  increased,  however,  a  second  branch  becomes  accessible,  shown  as  Branch 

II,  in  Figure  4:  While  the  cold  beam  growth  rate  continues  to  increase  with 

increasing  current,  the  intrinsic  energy  spread  stabilization  mechanism  becomes 

more  effective  due  to  the  increase  in  | ot I  until,  near  n  <  1/2  a  very  small 

energy  spread  results  in  a  large  spread  in  angular  velocity;  the  instability  is 

therefore  easily  stabilized.  More  simply  said,  we  have  stability  on  the  low 

current  branch  (Branch  I)  because  the  growth  rate  is  small  and  on  the  high 

current  branch  (Branch  II)  because  the  stability  mechanism  is  strong. 

There  is  a  third  branch  which  appears  in  the  example  of  Figure  2  above 

ns  =  1/2  (I  >  316  A)  for  b  =  10  and  20  and  is  illustrated  as  Branch  III  in  Figure 

4.  This  region  is  accessible  only  in  the  modified  betatron  since  we  are 

constrained  in  the  conventional  betatron  by  the  equilibrium  condition  n  <  1/2. 

s 

For  b  =  5  in  the  example  shown  in  Figure  2  the  stable  points  fall  to  the  right  of 
the  dashed  line  and  so  are  not  shown.  As  the  current  is  increased  beyound  316  A 
the  growth  rate  increases  and  the  stability  mechanism  becomes  less  effective 
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as  la!  decreases.  Consequently  as  the  current  is  increased  beyond  this  point  the 
required  energy  spread  increases  monotonically .  This  is  illustrated  by  the  third 
branch  in  the  figure. 

Finally  we  comment  on  the  accessibility  of  the  various  branches  available 
for  small  energy  spreads.  If  beam  injection  proceeds  slowly,  over  many  growth 
times,  say,  it  appears  that  only  the  lowest  branch  is  accessible;  attempting  to 
add  more  current  will  drive  the  system  unstable.  If  however  current  can  be 
introduced  into  the  accelerator  more  rapidly,  the  higher  current  branches  may 
become  accessible.  Only  a  carefully  designed  experiment  can  test  this 
speculation.  For  b  **  5  and  10  for  the  parameters  of  Figure  2  typical  growth 
times  are  ~  3  particle  circulation  periods  so  that  high  current  injection  on  this 
time  scale  is  a  practical  experimental  possibility.  The  third  branch  in  Figure  2 
is  clearly  the  most  promising  for  very  high  current  operation. 


V.  Conclusions 

We  have  shown  that  the  addition  of  a  toroidal  magnetic  field  to  a 


conventional  betatron  may  significantly  improve  the  current  carrying  capacity  of 


the  betatron  by  controlling  the  collective  instabilities  which  limit  the  current. 


The  calculation  has  included  self  field  effects  and  a  simple,  though  realistic 


model  for  the  longitudinal  and  transverse  impedances.  The  stabilizing  effects  of 
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betatron  oscillations,  which  have  been  shown  to  include  the  effects  of  emittance, 
toroidal  geometry,  and  energy  shear  due  to  the  electrostatic  potential  drop 
across  the  beam,  become  stronger  as  the  toroidal  field  is  increased,  given  a 
fixed  beam  radius.  Inclusion  of  self  field  effects  In  the  stability  criterion 
(6)  has  been  shown  to  lead  to  a  multi-valuedness  in  the  current  vs  energy  spread 
plot  which  has  been  interpreted  as  the  result  of  the  competition  between  the 
growth  and  stabilization  mechanisms.  Accessibility  of  the  high  current  branches 
may  depend  on  the  duration  of  the  injection  process. 


—3 
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Figure  Captions 


Figure  1.  Modified  Betatron  configuration. 

Figure  2.  Instability  threshold  current  vs  energy  spread  for  b  =  0,  5,  and 
10.  The  beam  and  chamber  parameters  are  rQ/ a  =  6,  a/r^  = 

3,  Y  =  3,  R0/Z0  =  4,  Q  =  10. 

Figure  3.  Same  as  Fig.  2  except  y  =  6. 

Figure  4.  Sample  plot  of  limiting  current  vs  energy  spread,  illustrating  the 
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By  adding  an  {  =2  stellarator  field  to  a  betatron  accelerator,  a  new  configuration  is  ob¬ 
tained  which  is  capable  of  accelerating  multikiloamp  beams  and  which  will  tolerate  a  large 
(more  than  50%)  mismatch  between  the  particle  energy  and  the  vertical  magnetic  field.  The 
additional  field  is  a  twisted  quadrupole  which  acts  as  a  strong-focusing  system.  Hits  de¬ 
vice  has  been  analyzed  both  analytically  and  numerically. 

PACS  numbers:  52.75.Di,  29.20.Fj 


Conventional  betatrons1-*  are  current  limited  at 
injection.  Recently,  efforts  have  been  made  to  ex¬ 
tend  the  current-carrying  capability  of  the  beta¬ 
tron.  For  example,  the  plasma  betatron3  employs 
a  toroidal  magnetic  field  in  the  direction  of  the 
particle  orbit  to  contain  the  plasma.  Current  In¬ 
terest  is  focused  on  high- current  nonneutral 
electron  acceleration  in  modified  betatrons.4'® 

By  adding  a  stellarator  field  to  a  cyclic  accel¬ 
erator,  a  strong-focusing  system7  is  obtained 
which  can  sustain  high  currents  and  large  mis¬ 
match  between  particle  energy  and  vertical  field. 
The  energy  bandwidth  relaxes  the  design  require¬ 
ments  for  the  injector  and  the  magnetic  field  sys¬ 
tem.-  Unlike  fixed-field  alternating-gradient  beta¬ 
trons,®  the  stellar ator-betatron  (or  stellatron)  in¬ 
cludes  a  strong  toroidal  field  to  confine  very  high 
currents.  Figure  1  shows  a  sketch  of  the  stel¬ 
latron  configuration. 

We  have  quantitatively  studied  the  stellatron 
configuration.  Our  studies  have  consisted  of  nu¬ 
merical  and  single-particle  orbit  calculations, 
as  well  as  analytical  linearized  orbit  theory,  in¬ 
cluding  the  beam  self-fields. 

We  may  study  the  behavior  of  an  intense  elec¬ 
tron  beam  in  the  stellatron  quantitatively  by  con¬ 
sidering  small  departures  from  a  “reference 
orbit,”  a  circle  located  at  the  null  point  in  the 
quadrupole  field,  at  r  =  r0,  z  =  0.  Here  and  below 
we  use  a  cylindrical  (r,6,z)  coordinate  system 
with  origin  at  the  center  of  the  torus’s  major 
cross  section.  Quantities  evaluated  at  the  ref¬ 
erence  orbit  will  carry  a  subscript  0  below;  de¬ 
partures  from  this  orbit  will  carry  a  subscript  1. 


The  twisted  quadrupole  field,  of  period  2 jt/to, 
then  is  written  as 


Br  s  kBl(-r1sinmO+zl  cos m0)t 

BM  =  kBt{rlcosmO+zlshvnO)t  Be^Beoj 


(1) 


Where  k,BttB0o  are  constants,  and  the  betatron 
field  is 

Bre-nBj0Zi/r0,  B,sBt0(l-n(r1/r0)],  (2) 

where  B ^  is  the  vertical  field  at  the  reference 
orbit  and  n  Is  the  usual  field  index. 

We  consider  the  motion  of  an  electron  located 
within  a  beam  whose  center  is  located  at  r  =r0 
+Ar,  z=Az;  the  electron’s  position  is  r=r0 
+Ar  +  &r  =  r0+r,,  z=Az  +  6z£z,.  Using  a  cylin¬ 
drical  approximation  for  the  beam  self-fields,  we 
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find 


n  +  fi«,U-«  +  Mcosm0)rl-.|^V^6r+^Arj  =  -  {ntoMOa slnme)zl  +  neoz1),  (3a) 

1  w  ^  ^  \  t 

3'i  +  ««)*(”-MCO8w0)z1---^rf62+-^-Azj  =  -(pn,0asln mO)rt-Qeplt  (3b) 


fil=-(i'l/y0)n*o» 

where  toMO=eBt0/my0c,  y0  Is  the  relativistic  factor 
for  the  reference  orbit,  o>*a  =  4 irn#a/my0J  n0  is 
the  beam  number  density,  r*  is  the  minor  radius 
of  the  beam,  a  is  the  minor  radius  of  the  (perfect¬ 
ly  conducting)  chamber,  and  p  =fer0B*/B,0. 

By  performing  a  formal  ensemble  average  of 
(3a)-(3c)  one  may  find  equations  governing  single¬ 
particle  motion  and  that  of  the  beam  centroid. 
Details  will  be  published  elsewhere.  By  changing 
the  independent  variable  from  t  to  0,  and  making 
the  transformation  £  =  (r1  +  iZ1)/r0=tp  exp(im0/2)t 
one  obtains  an  equation  for  ip  which  may  be  solved 
in  the  special  case,  n  -  { ,  with  the  following  re¬ 
sults. 

Particle  motion  is  oscillatory  (under  certain 
conditions;  see  below)  about  a  center  located  at 

*L  _ _  (4) 

r„  R>  +  pa(»i2  +  m6-fi>)*1  »  '  ' 

where  ft,  =£  -n„  n,  =ioba/2y0aQMOaf  b=Beo/BM, 

A  —Yi/P<?Yot  and  angular  brackets  denote  an  en¬ 
semble  average.  There  are  five  characterisitc 
oscillation  frequencies,  and  (m/2±vt 
where 

vta*n+\ma±(fimat  ii2)l,x  (5) 

with  n  =ft*  +56*,  m  =  m  +  b.  These  frequencies 
are  real  when  the  system  is  located  within  the 
regions  of  the  plane  of  Fig.  2  marked  "stable.” 

We  remark  that  for  low- current  beams  (n0—  0) 
the  stability  condition  reduces  to 

|£m2  +  m&  -  l|  >|2p|.  (6) 

The  "most”  stable  configuration  results  when  the 
field  lines  are  twisted  clockwise  ( m  >0)  when 
viewed  in  the  direction  of  Be8,  i.e.,  in  the  same 
sense  as  electron  gyration  about  B0. 

Similarly,  the  motion  of  the  beam  center  is  it¬ 
self  oscillatory  about  a  center  located  at 

—  - - -  (7) 

r0  +  pa(ma  +  mb -ft*)"1  *  '  1 

where  ft*  =5  -  (r*  a/aa)nt ,  with  characteristic  fre¬ 
quencies  as  in  (5),  under  the  replacement  n* 

-  (rb*/aa)n,. 

Two  important  features  of  the  solution  are  worth 


(3c) 


'pointing  out.  First,  stable  motion  is  possible 
throughout  an  injection- acceleration  cycle.  This 
has  been  checked  for  many  possible  time  histo¬ 
ries.  A  typical  trajectory  in  the  stability  plane  is 
shown  in  Fig.  2.  The  unstable  region  on  the  left 
of  the  diagram  would  not  be  entered  in  this  case 
even  if  the  acceleration  were  continued;  “ u ” 
never  changes  sign  in  this  case. 

The  second  important  feature  of  the  solution 
pertains  to  the  energy  bandwidth  of  this  machine. 
We  note  that  the  radial  shift  (7)  of  the  orbit  of  a 
mismatched  beam  is,  as  expected,  much  smaller 
than  that  in  a  weak  focusing  (p  =  0)  device,  (p 
can  easily  exceed  100-200  in  designs  we  have 
considered.)  The  stellatron’s  large  energy  band¬ 
width  has  very  helpful  consequences  for  injector 
and  magnetic  field  design  tolerances. 

The  introduction  of  fixed  toroidal  and  helical 
fields  to  the  betatron  causes  the  betatron  wave¬ 
lengths  to  depend  on  energy,  resulting  in  reso¬ 
nant  instabilities  driven  by  field  errors  during 
acceleration.  If  the  toroidal  field  is  sufficiently 
large,  the  betatron  wavelengths  will  be  insensitive 
to  beam  current.  Such  instabilities  may  be  avoid¬ 
ed  by  holding  all  the  fields  in  constant  ratio  dur¬ 
ing  acceleration.  Alternatively,  the  effect  of  the 


FIG.  2.  Stellatron  stability  plane  (n  ).  The  dotted 
line  is  the  trajectory  of  an  experiment  with  /  =10  kA, 
£0(  =  5  kG,  e,  s  2fi/mb  =1,  m  =20,  r*  =  l  m,  while  £,* 
is  raised  from  118  to  1700  G,  corresponding  to  an  in¬ 
crease  in  energy  from  ~  3.5  to  50  MeV. 
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instabilities  may  be  minimized  if  the  energy  gain 
per  revolution  is  large  enough  to  pass  rapidly 
through  each  resonance. 

A  single-particle  code,  which  integrates  the 
relativistic  equations  of  motion  for  an  electron  in 
an  applied  magnetic  field,  has  been  utilized  to 
study  certain  nonlinear  aspects  of  the  stellatron 
configuration.  Unlike  the  analytical  analysis  of 
the  preceding  paragraphs,  this  analysis  does  not 
employ  a  paraxial  approximation  for  the  electron 
motion  and  does  not  use  an  expansion  in  the  parti¬ 
cle  displacement  from  a  reference  orbit.  Also, 
the  applied  field  in  this  analysis  includes  toroidal 
corrections  to  first  order  in  the  inverse  aspect 
ratio. 

The  total  magnetic  field  utilized  by  the  code 
may  be  expressed  as  B  =  Bt  +B, ,  where  Bfc  is 
the  conventional  betatron  field,  given  by  Eq.  (2), 
and  B,  is  the  stellarator  field,  given  by  B,  =V4>4 , 
in  terms  of  the  magnetic  scalar  potential,  4>, , 
which  may  be  expressed  as  $<<0>+4><a>,  where 

**<0,(P»^>s)”®eo{s  +  (ei/af^i(*)s*nW(<P  -  <*$)]}» 

Here,  x=fap,  a  =  2v/L,  L  is  the  helix  pitch 
length,  and  /,  represents  the  modified  Bessel 
function.  The  coordinates  (p,<p,s)  form  a  local 
cylindrical  system  centered  on  the  minor  axis, 
where  s  =Ro0  is  distance  measured  along  the 
minor  axis  for  toroidal  angle  6,  and  (p,9>)  are 
polar  coordinates  in  the  plane  transverse  to  the 
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FIG.  3.  Stellatron  single-particle  bandwidth.  The 
bandwidth  Au/u{,  where  u  -yP,  is  plotted  against  c, 
a  2ft /mb.  The  accelerator  is  matched  at  y"7  with  B, 0 
•=118  G,  and  >B0o  =  l  kG. 


minor  axis  at  s.  The  toroidal  correction,  <ptU) , 
is  given  to  first  order  in  the  inverse  aspect  ratio 
by  Danllkln.®  All  calculations  have  been  per¬ 
formed  for  1  =  2.  The  variables  p  and  m,  which 
describe  the  helical  field  in  the  previous  analyti¬ 
cal  analysis,  are  given  by  p  =  emb/ 2  and  m 
=  2aR0for  1  =  2. 

This  model  has  been  utilized  to  investigate  the 
single-particle  bandwidth  of  the  stellatron.  As 
e,  increases,  the  allowed  mismatch  in  the  stel¬ 
latron  becomes  too  large  to  be  correctly  modeled 
by  the  linearized  theory.  Figure  3  shows  the  re¬ 
sults  from  both  models  for  bandwidth  versus  e, . 
These  calculations  assume  a  torus  having  a  1-m 


FIG.  4.  Single-particle  orbits,  (a)  Without  the 
helical  field  components  (ej  =  0),  A«/#,=  ±3%j  (b)  stel¬ 
latron  orbit  with  c,  *=i  ,  Att/u0*=  ±  50%. 
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major  radius  and  a  10- cm  minor  radius.  A  test 
electron  is  launched  tangent  to  the  minor  axis 
with  relativistic  momentum  u=  y0,  which  differs 
from  the  matched  momentum,  u0>  by  varying 
amounts,  An.  The  figure  shows  A n/n0  against  e, , 
and  represents  the  maximum  |An/n0|  for  which 
the  test  orbit  remained  confined  in  the  device. 
Mismatch  in  excess  of  50%  can  be  tolerated  for 
these  parameters. 

Figure  4(b)  shows  typical  stellatron  orbits,  pro¬ 
jected  on  the  minor  cross  section,  for  Beo  =  5 
kG,  e,  for  ±  50%  mismatch.  The  betatron  field 
is  again  118  G  with  n  =5.  Without  the  helical  field 
'contribution,  Fig.  4(a)  shows  that  as  little  as 
±  3%  mismatch  is  not  tolerabie. 

The  superposition  of  twisted  quadrupole,  toroi¬ 
dal,  and  conventional  betatron  magnetic  fields  ap¬ 
peals  to  offer  significant  practical  advantages 
for  the  confinement  and  acceleration  of  large 
electron  currents  (tens  of  kiloamperes)  to  moder¬ 
ate  energies  (hundreds  of  megaelectronvolts). 
Foremost  among  these  advantages  is  the  greatly 
improved  energy  bandwidth  over  that  of  a  weak- 
focusing  device.  The  large  bandwidth  of  the  stel¬ 
latron  relaxes  the  requirements  for  monoener- 
getic  injection,  for  a  uniform  (within  a  few  per¬ 
cent)  magnetic  field  configuration,  and  for  a  rigid 
mechanical  design.  Injection  should  not  be  any 
more  difficult  than  for  other  high-current  ac¬ 
celerator  concepts,  and  is  facilitated  by  the  ex¬ 
ternally  applied  rotational  transform  of  the  stel- 
lerator  field.  The  orbits  should  remain  stable 


from  injection  up  to  the  highest  energies  achiev¬ 
able  by  conventional  inductive  acceleration. 
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ABSTRACT 


The  combination  of  a  bumpy  torus  field  and  a  conven¬ 
tional  betatron  field  leads  to  an  interesting  strongly- 
focused,  high-current  accelerator  configuration.  We  discuss 
the  single  particle  orbital  stability  question  and  show  that 
the  strong-focusing  in  this  accelerator  can  easily  lead  to 
greater  than  20%  bandwidth  in  allowed  mismatch  between  the 
vertical  magnetic  field  and  the  average  beam  particle  energy. 


Conventional  betatrons  are  current-limited  due  to 

the  space  charge  at  injection.  Until  recently,  the  approach 

to  multi-kiloampere  betatrons  has  concentrated  exclusively 

on  overcoming  this  space  charge  limit.  High  current  conven- 

2 

tional  betatrons  employ  high-energy  injectors,  while  modi- 
tied  betatrons  '  employ  a  toroidal  magnetic  field  to  prevent 
space  charge  blow-up.  In  both  of  these  cases  however,  a  mis¬ 
match  between  the  injection  energy  and  vertical  field  of  a 
few  percent  will  cause  the  beam  to  hit  the  wall.  The  maximum 
allowed  error  in  the  vertical  field  is  typically  on  the  order 
of  a  few  gauss.  Recently  it  was  shown5  that  the  combination 
of  an  %  =  2  stellarator  and  betatron  field  results  in  a  strong 
focusing  high-current  betatron  or,  "stellatron, "  with  a  large 
energy  bandwidth.  Such  a  configuration  has  the  advantages  of 
relaxing  the  vertical  field  and  injector  tolerances.  In  addi¬ 
tion,  the  strong  focusing  introduces  a  threshold  for  the 
negative  mass  instability,  so  that  this  instability  does  not 
operate  at  injection.  In  this  note  we  report  analytical  and 
numerical  results  on  the  bandwidth  and  stability  of  an  alter¬ 
native  strong-focusing  scheme,  namely,  a  combination  "bumpy 
torus"  and  betatron  field,  corresponding  to  the  i  =  0  stellatron 
The  bumpy- torus  betatron  field  consists  of  a  super¬ 
position  of  an  Z  =  0  stellarator  field  and  the  field  of  a 
conventional  betatron.  Near  the  minor  axis  at  r  =  rQ,  z  =  0, 
this  field  has  the  form 


Be  =  Beo 


(■*5 


cos  m0 


Bz  =  Bzo(1  -  nx^  +  2  6Be  s;*-n  m8 


where  x  =  (r-rQ)/ro,y  =  z/rQ,  0  is  the  azimuthal  angle, 

n  is  the  vertical  field  index,  and  m  is  the  number  of  bumpy- 

torus  field  periods  around  the  torus.  B  ,  BQ  ,  and  6B_ 

zo  0o  6 

are  constants. 

In  the  paraxial  approximation,  the  equation  for  parti¬ 
cle  motion  is,  for  n  = 


+  JL  \2  - 

d<j>2  m2 


4n  +  b2  (1  +  e  cos2(j))  2  ip 
s 


=  jL  e  2m 


( 2 4>  +  e  sin  2<J>) 


m2  Po 


where  <J>  =  md/2,  ip  =  (x+iy)exp 


[(i/2)/ 


b (1  -  e  cos  m0 )  d0 


b  =  B0O/BZO»  e  =  CBfi/Bflo»  P„  Is  tlie  momentum  of  a  particle 


0'~0o 


which  would  circulate  on  the  minor  axis,  Sp  is  the  "momen¬ 
tum  error,"  and  the  self-fields  are  included  via  n  e  w2/ 

s  b 

(2Yo  Ko)  Where  wb'  ^zo  are  ^eam  plasma  frequency  and 

the  vertical-field  cyclotron  frequency,  respectively. 

Eq. (2)  is  a  Hill  equation,  which  has  characteristic 
bands  of  stability,  as  shown  in  Figure  1.  The  shaded  regions 
in  the  figure  are  unstable  portions  of  the  plane,  e  vs.  b/m, 
for  the  case  ng  =  30  and  m  =  30.  The  intersections  of  the 
unstable  regions  with  the  abscissa  are  given  by 

(b2  +  2  -  4n  )/m2  =  q2,  where  q  =  0,  1,  2,  .... 
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As  we  increase  B  to  accelerate  the  electrons,  we 

z 

typically  would  not  wish  to  increase  Bq .  The  result  is  that 
the  operating  point  of  the  accelerator  will  tend  to  move  from 
right  to  left  in  Figure  1.  Consequently,  the  accelerator 
should  be  run  in  the  left-most  stable  band  to  avoid  crossing 
unstable  bands.  These  considerations  require  m>b  at  injection 
and  force  the  use  of  a  large  number  of  field  periods  in  the 
design  of  the  strong-focusing  system. 

The  left-most  unstable  band,  corresponding  to  q  =  0, 

is  due  to  beam  space-charge  and  rapidly  disappears  during 

acceleration  since  the  self-field  index,  n  ,  is  proportional 

s 

to  y-3,  where  y  is  the  relativistic  factor.  The  left-most 
stable  band,  therefore,  becomes  broader  during  acceleration. 

In  conventional  betatrons,  resonances  are  avoided  by 
increasing  y  and  the  accelerating  field  in  synchronism.  The 
introduction  of  non- synchronous  fields,  as  in  the  modified 
betatron  and  the  stellatron,  makes  the  betatron  wavelengths 
energy-dependent,  which  can  lead  to  the  crossing  of  resonant 
instabilities  driven  by  field  errors  during  acceleration.  As 
in  all  strong-focusing  devices,  the  occurence  of  orbital 
resonances  plays  an  important  role  in  the  operation  of  the 
bumpy- torus  betatron.  Using  the  Floquet  solutions  to  (2) 
it  is  possible  to  obtain  a  condition  for  the  integer  reso¬ 
nances,  when  space-charge  effects  may  be  neglected: 


where  i|>  ( 4> )  is  the  solution  to  (2)  with  6P  =  0  satisfying 
4>j(o)  =  1,  ’  (o)  =  0  where  n  is  an  integer,  the  Fourier 

component  number  of  the  dipole  field  error.  Equation  (3) 
provides  the  basis  for  numerical  calculation  of  contours  in 
the  stability  plane  where  (3)  is  satisfied  for  a  given  n; 
an  example  is  given  in  Figure  2. 

If  all  the  fields  cannot  be  made  synchronous  with  the 
particle  energy,  the  effect  of  resonant  instabilities  could 
be  minimized  by  making  the  energy  gain  per  pass  large. 

Finally,  we  consider  the  question  of  containment  of 
beams  whose  average  momentum  is  not  matched  to  the  vertical 
betatron  field,  i.e.  the  question  of  the  momentum  compaction 
of  this  configuration.  Figure  3  shows  the  allowed  mismatch, 
Ap/Pq,  plotted  against  e  =  6B0/Bg  for  B0o  =  2kG,  Bzq  =  118G, 
n  =  ^-,  rQ  =  100cm  and  m  =  30.  This  plot  is  generated  numeri¬ 
cally  by  launching  particles  on  the  minor  axis  along  the 
toroidal  direction  with  various  amounts  of  mismatch.  The 
figure  shows  the  largest  mismatch  for  which  the  calculated 
orbits  are  contained  in  a  10cm  minor  radius  chamber.  Con¬ 
tainment  of  particles  with  a  mismatch  of  +20%  is  obtained 
for  e^O. 2. 

In  conclusion,  we  find  the  spatially  alternating  trans¬ 
verse  magnetic  field  associated  with  a  bumpy-torus  provides 
an  alternating  field  gradient  on  the  minor  axis,  which  leads 
to  a  potentially  interesting  strongly-focused  accelerator 
configuration.  This  new  accelerator  is  seen  to  have  a  region 
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of  stable  orbits,  and  to  have  a  significant  bandwidth  in 
allowed  mismatch  between  the  vertical  magnetic  field  and  the 
particle  momentum. 

This  work  was  supported  by  the  Naval  Research  Labora¬ 
tory.  We  wish  to  acknowledge  discussions  with  members  of 
the  Advanced  Accelerator  Project  at  NRL. 


References 

1.  See,  e.g.  D.W.  Kerst,  Handbuch  der  Physik,  XLIV,  13  (1959) . 

2.  E.  Lee  et  al,  Proc.  1983  Particle  Accelerator  Conference 
(to  be  publTshed) . 

3.  P.  Sprangle  and  C-A.  Kapetanakos,  J.  Appl.  Phys.  49. 

1  (1978) . 

4.  N.  Rostoker,  Comm.  Plasma  Phys.  (5,  91  (1980)  . 

5.  C.  Roberson,  A.  Mondelli,  and  D.  Chernin,  Phys.  Rev.  Lett. 
50,  507  (1983). 


39.7 


LIST  OF  FIGURES 


Figure  1: 


Figure  2: 


Figure  3: 


Stability  plane  for  bumpy-torus  betatron,  for 
the  case  n  =  m  =  30.  The  shaded  regions  are 
unstable  for  particle  motion. 

Stability  plane  for  bumpy-torus  betatron,  with 

the  single  particle  resonance  lines  n  =  0,  5, 

10,  15,  20,  25,  indicated  for  the  case  n  =  0, 

s 

m  =  30. 

Single  particle  bandwidth.  Data  points  indicate 
the  maximum  value  of  momentum  mismatch  tolerated 
by  the  device  vs.  the  bump  size,  e,  for  particles 
initialized  on  the  minor  axis,  for  the  specific 
case  Bzo  =  118G,  BQo  =  2kG,  rQ  =  100cm,  m  =  30. 
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Program  LBE 


This  program  uses  the  IMSL  subroutine  DVERK  to  integrate  the 
linearized  equations  of  motion  for  either  the  beam  centroid 
or  for  individual  particles  within  a  beam  in  a  modified 
betatron. 
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particle  equations 
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Chamber  major,  minor  radii,  field  index 
2 

Current  (kA)  ,  beam  minor  radius 

Dimensionless  correction  quantities  added 
to  &E  and  JlB,  respectively. 

(Logarithms  appearing  in  toroidal  correction 
terms . ) 

(JlE  =  Jln(a/rb)+  AE;  Jig  =  Jln(a/rb)+  1  +  AB) 

.  3 

Dimensionless  time  step  ,  total  number  ^ 
of  steps,  output  file  writing  frequency  , 
integration  error  parameter 

Initial,  final  vertical  field  values; 
acceleration  time^ 

Initial,  final  toroidal  field  values; 
time  for  change  of  Bfl^ 


Notes: 


All  quantities  in  cgs  units  unless  otherwise  specified. 

2 

If  I  is  entered  as  a  negative  number,  the  density  used 
in  calculating  ng  is  ( }  1 1/17)  /tsx^2xq,  as  for  positive  I, 
but  all  toroidal  corrections  to  coefficients  are  set 

to  zero . 

3 

Time  step  =  Dt/ (eBQ . /my • c) . 

oi  i 

4 

e.g.  If  KPT  =  2,  output  is  written  (to  file  FOR003.DAT) 

every  other  time  step. 
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A  linear  variation  in  time  is  used. 
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APPENDIX  BB 


Program  BTRAK 


This  program  uses  the  IMSL  subroutine  DVERK  to  integrate  the 
single  particle  equations  of  motion  in  various  accelerator 
field  configurations. 


ro 

Bzo'  n'  n2 '  n3 

b4>o 

KSTEL 

Model 


Input  Data  File1:  FOR025.DAT 
:  Major  radius 

:  Vertical  magnetic  field  and  field 

indices  at  r^ 

:  Toroidal  field  at  r 

o 

:  (jZf/l)  =  (no,yes)  stellarator  field  present 

:  (if  KSTEL  =  1,  for  stellarator  field3) 

4 

:  Field  error  model  :  0  =  No  field  error 


1  =  Cylindrical  limit 

(r£C°) 

2  =  (r,z)  cylindrical 

coordinates 


SB,  m,  H  : 

SB,  a,  XrQ,  m,  NZ  : 
Dx,  KTOT,  TOL,  MXOUT: 


Model  =  1 

Model  =  2  / 

.  5 

Dimensionless  step  size  ,  total  number  of 
steps,  integrator  error  control,  max 
number  of  points  written  to  output  file. 


Sr,  6z,  <p 

:  Initial  values 

(Sr,max'  (6z)it,ax 

!  If  l5rl>Srmax  1 

tion,  program 

Y.  fr.  £z 

:  ^initial'  fr  = 

max 


Notes: 


All  quantities  in  cgs  units 
2fiz  =  Bzo(1_nx  +  |n2x?  +  |<n-n2)y2  +  n3xV6 

-  |-(n+n2+n3)  xy2 ) 
x  =  (r-r0)/r0;  y  »  2/r0 
3 

Br  =  Bzoli  [-x  siranij)  +  y  cosm<J>] 

Bz  =  Ix  cosm<j>  +  y  sinmcj>] 


Model  1: 


B  = 
r 


B*  - 


B  = 
z 


6BIm[(k5r  +  ik6z)  £-1eiln<l>] 
SB^-Rel  (k6r  +  ik<$z) 
6BRe[(k6r  +  ik6z)  £“1eim<i>] 


Br  " 


B4>  " 


2-6Bk6r  sinm<j> 
6B  cosmij) 


Bz  = 


^-SBkiSz  sinm<J> 


where  k  =  m/r 
Model  2: 


Br  - 


6B 


■j^r-  [a  +  2£x]  3  (Xz)  sin  mcj) 


if  Z  /  0 


if  5L  «  0 


B<f>  =  IF-1111  I1  +  ax  +  &x2]  js  (Az)cos  m<j> 
y  o 

Bz  =  5B[1  +  ax  +  $x2  ]  z '  (Xz)  sin  m<|) 


where  x  =  (r-r  ) /r 
o  '  o 

3  =  |-(m2  -  a  -  (XrQ)  2) 

2(Xz)  =  sinh(Xz)  NZ  =  1 
2 ( X z )  =  cosh (Xz)  NZ  ^  1 

5dt  =  t*dT  where  t*  =  2Tr/max(ft  ,o  ). 
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APPENDIX  CC 


Program  TUNES 

This  program  checks  for  resonances  of  the  form 

n+v+  +  n_v_  =  P 

where  v+  are  the  tunes  of  the  fast  and  slow  modes  and  n_^,  P 
are  integers.  The  program  checks  those  n_^_  satisfying 
In.  |  +  In  |  <  3  and  for  P<P _ _  where  P__  is  input. 

Input  Data  File1:  FORj355.DAT 
I FLAG  :  (o,l)  =  (beam,  particle)  motion 

rQ,  a,  r^  :  Chamber  major,  minor  radii,  beam  minor  radius 
I  :  Current  (kA) 

AE,  AB  :  Toroidal  corrections:  =  fcn£t/rb)+  AE; 

iQ  =  £n  ^i/r^)  +  1  +  AB 
2 

Bzi,  B zf,  n:  Initial,  final  vertical  field  values  ,  field  index 

.  .  2 
B0i'  B0f  :  Initial,  final  toroidal  field  values  . 

P„,„„  :  Maximum  value  of  P  to  be  checked, 

max 


^All  quantities  in  cgs  units,  unless  otherwise  specified 
2 

Linear  variation  in  time  used.  (T:  0-1) 


Notes : 


Interactively  Generated  Plots 
(Uses  JAYCOR  plotting  routine) 


Response  to  plot  option?  request: 


No  plot  (Stop) 


v+,  v_  plane,  showing 
resonance  lines 


v+  vs  time 


v_  vs  time 
y  vs  time 


Output  File:  FOR057.DAT 


List  of  resonances  crossed,  with  values  of 
v+,  v_,  T,  y,  Bz,  and  BQ  at  crossing. 
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